Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,520,645 C→T E247K (GAA→AAA)  SAUSA300_RS07380 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,520,6450CT98.6% 191.6 / ‑3.8 70coding (739/1245 nt)SAUSA300_RS07380hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (2/67);  minor base . (0/1);  total (2/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CAAACATCCAGTTGTATACATTAACTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTCATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTGAT  >  NC_007793/1520553‑1520741
                                                                                              |                                                                                                
gatctaTCCCGTTGTATACATTAACTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTtttttattttca                                                                                           <  1:5097711/95‑1 (MQ=255)
  ctccaccAGTTGTATACATTAACTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAAT‑TTTTttttattttcat                                                                                          <  1:3461444/95‑1 (MQ=255)
        atgTTGTATACATTAACTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGa                                                                                   >  1:4660922/3‑100 (MQ=255)
           ttGTATACATTAACTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTg                                                                                <  1:2042496/100‑1 (MQ=255)
            tGTATACATTAACTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTgg                                                                               <  1:1633190/100‑2 (MQ=255)
                  cATTAACTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATAc                                                                         <  1:832116/100‑1 (MQ=255)
                  cATTAACTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATAc                                                                         <  1:2788619/100‑1 (MQ=255)
                   aTTAACTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACa                                                                        <  1:3311575/100‑1 (MQ=255)
                   aTTAACTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACa                                                                        <  1:1726189/100‑1 (MQ=255)
                      aaCTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTGCCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGt                                                                     <  1:2903547/100‑1 (MQ=255)
                      aaCTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGt                                                                     <  1:135063/100‑1 (MQ=255)
                      aaCTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGt                                                                     <  1:3298808/100‑1 (MQ=255)
                      aaCTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGt                                                                     <  1:2062812/100‑1 (MQ=255)
                      aaCTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGt                                                                     <  1:5117092/100‑1 (MQ=255)
                      aaCTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGt                                                                     <  1:237844/100‑1 (MQ=255)
                       aCTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAAT‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGt                                                                     <  1:826803/100‑1 (MQ=255)
                       aCTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTa                                                                    <  1:4709891/100‑1 (MQ=255)
                        ctctTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAATTTTTTTTTTATTTTCATATAATGATTGTAAATACAAGt                                                                     <  1:157741/100‑1 (MQ=255)
                          ctTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGcc                                                                 <  1:1499343/100‑1 (MQ=255)
                           tttATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCt                                                                <  1:2732878/100‑1 (MQ=255)
                           tttATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCt                                                                <  1:1711404/100‑1 (MQ=255)
                            ttATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCtt                                                               <  1:2227531/100‑1 (MQ=255)
                            ttATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCtt                                                               <  1:3706229/100‑1 (MQ=255)
                            ttATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCtt                                                               <  1:4732258/100‑1 (MQ=255)
                             tataAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCtgg                                                              <  1:3322469/100‑3 (MQ=255)
                             tataAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTg                                                              <  1:4104387/100‑1 (MQ=255)
                              ataAAATTGACTTAGTCGTAATCCTACTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCttag                                                             <  1:4155/100‑3 (MQ=255)
                                      gACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATc                                                     <  1:3381795/100‑1 (MQ=255)
                                      gACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATc                                                     <  1:422734/100‑1 (MQ=255)
                                      gACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATc                                                     <  1:5307818/100‑1 (MQ=255)
                                      gACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATc                                                     <  1:3516832/100‑1 (MQ=255)
                                      gACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TGTTTTTTATTTTCATATAATGATTGAAAATACAAGTAGCCTTGAATATaagc                                                     <  1:3805525/100‑3 (MQ=255)
                                       aCTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCa                                                    <  1:2340362/100‑1 (MQ=255)
                                       aCTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCa                                                    <  1:3838770/100‑1 (MQ=255)
                                        cTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAAtatg                                                   <  1:5499601/100‑4 (MQ=255)
                                         ttAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAgg                                                  <  1:845994/100‑1 (MQ=255)
                                           aGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGag                                                <  1:3430125/100‑2 (MQ=255)
                                           aGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGAATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAggcg                                                <  1:5673344/100‑3 (MQ=255)
                                            gTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGCAta                                               <  1:2265131/100‑2 (MQ=255)
                                              cGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAgcacagt                                             <  1:1654671/100‑7 (MQ=255)
                                               gTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCttt                                            <  1:2728183/100‑1 (MQ=255)
                                               gTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCttt                                            <  1:5217161/100‑1 (MQ=255)
                                               gTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCttt                                            <  1:2373409/100‑1 (MQ=255)
                                               gTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCttt                                            <  1:1153208/100‑1 (MQ=255)
                                               gTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCttt                                            <  1:1450823/100‑1 (MQ=255)
                                               gTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCttt                                            <  1:2798421/100‑1 (MQ=255)
                                               gTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCttt                                            <  1:1706940/100‑1 (MQ=255)
                                               gTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCttt                                            <  1:5709438/100‑1 (MQ=255)
                                               gTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCttt                                            <  1:3958930/100‑1 (MQ=255)
                                                           tACCCGTTTCAATTGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGt                                <  1:1825305/100‑1 (MQ=255)
                                                           tACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGt                                <  1:4563833/100‑1 (MQ=255)
                                                           tACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGt                                <  1:4270717/100‑1 (MQ=255)
                                                           tACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGt                                <  1:4807674/100‑1 (MQ=255)
                                                           tACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGt                                <  1:4852863/100‑1 (MQ=255)
                                                           tACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGt                                <  1:1054604/100‑1 (MQ=255)
                                                           tACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGt                                <  1:1491987/100‑1 (MQ=255)
                                                           tACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGt                                <  1:2724870/100‑1 (MQ=255)
                                                           ggcccGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGt                                <  1:3475165/98‑1 (MQ=255)
                                                             ttccgatctaTGGCATCTGGATAA‑‑TTTTTTTTATTTTCCTATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTat                              <  1:939566/91‑2 (MQ=255)
                                                             tccgatctaaTGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAg                              <  1:4055534/92‑1 (MQ=255)
                                                               cGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATGATCAGGATCTTTAGAAAGTAATGTAGTc                            <  1:5486754/100‑1 (MQ=255)
                                                               cGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTc                            <  1:5631855/100‑1 (MQ=255)
                                                               cGTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTc                            <  1:5376395/100‑1 (MQ=255)
                                                                gTTTCAATGGCATCTGGATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCa                           <  1:3937630/100‑1 (MQ=255)
                                                                             cTGGATAA‑‑TTTTTTT‑ATTTTTATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGct             <  1:2263721/100‑1 (MQ=255)
                                                                               ggATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGctt            >  1:5115900/1‑100 (MQ=255)
                                                                               ggATAA‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGctt            <  1:1658879/100‑1 (MQ=255)
                                                                                   ta‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTg        <  1:3949171/99‑1 (MQ=255)
                                                                                   aa‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTg        <  1:4606668/100‑1 (MQ=255)
                                                                                   aa‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTg        <  1:27238/100‑1 (MQ=255)
                                                                                   aa‑‑TTTTTTTTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTg        <  1:5639671/100‑1 (MQ=255)
                                                                                     acttttttttATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTg        <  1:4650432/98‑1 (MQ=255)
                                                                                       ttttttttATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAg      <  1:851579/100‑1 (MQ=255)
                                                                                       attttttTATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAg      <  1:68620/99‑1 (MQ=255)
                                                                                         ttttttATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTg    <  1:3849435/100‑1 (MQ=255)
                                                                                         ttttttATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTg    <  1:1001165/100‑1 (MQ=255)
                                                                                         ttttttATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTg    <  1:4372405/100‑1 (MQ=255)
                                                                                         ttttttATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTg    <  1:2219693/100‑1 (MQ=255)
                                                                                         ttttttATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTg    <  1:387294/100‑1 (MQ=255)
                                                                                         ttttttATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTg    <  1:2833820/100‑1 (MQ=255)
                                                                                         ttttttATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTg    <  1:3625359/100‑1 (MQ=255)
                                                                                          tttttattttcaTATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTgt   <  1:468741/100‑2 (MQ=255)
                                                                                          tttttattttcaTATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTgg   <  1:1220615/100‑2 (MQ=255)
                                                                                           ttttattttcatATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTgga  <  1:5061352/100‑3 (MQ=255)
                                                                                           ttttattttcatATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTGag  <  1:3702201/100‑2 (MQ=255)
                                                                                           ttttattttcatATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTGAt  <  1:4265341/100‑1 (MQ=255)
                                                                                              |                                                                                                
CAAACATCCAGTTGTATACATTAACTCTTTATAAAATTGACTTAGTCGTAATCCTTCTTTACCCGTTTCAATGGCATCTGGATAA‑‑TTTTTTTCATTTTCATATAATGATTGTAAATACAAGTAGCCTTGAATATAATCAGGATCTTTAGAAAGTAATGTAGTCATTATTTTAATTGCTTCTTGAGTGAT  >  NC_007793/1520553‑1520741

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑