Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,989,323 (A)8→7 100% coding (20/483 nt) yqeJ → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,989,3160A.80.1% 537.5 / 179.0 206coding (13/483 nt)yqeJuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (30/11);  new base . (76/89);  total (106/100)
Fisher's exact test for biased strand distribution p-value = 2.71e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

Only 100 of 225 total aligned reads displayed.

CGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGC  >  NC_000913/2989217‑2989414
                                                                                                   |                                                                                                  
cGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATtata                                                                                                    <  2:1936382/100‑1 (MQ=255)
cGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATtata                                                                                                    <  2:3388340/100‑1 (MQ=255)
 ggggATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATaa                                                                                                   <  2:1741731/100‑1 (MQ=255)
   ggATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATaaaa                                                                                                 <  1:1062975/100‑1 (MQ=255)
      tgatTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATaaaaaaa                                                                                              >  1:149404/1‑100 (MQ=255)
       gatTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATaaaaaaac                                                                                             <  1:935093/100‑2 (MQ=255)
        atTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATaaaaaaact                                                                                            <  1:3368935/100‑3 (MQ=255)
        atTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAAc                                                                                            >  1:904412/1‑100 (MQ=255)
        atTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAAc                                                                                            >  1:494447/1‑100 (MQ=255)
           gTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACttat                                                                                        <  1:3097758/100‑1 (MQ=255)
            ttATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACttatt                                                                                       >  2:1837828/1‑100 (MQ=255)
             tATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACttatta                                                                                      >  1:2996723/1‑100 (MQ=255)
                 gCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTAttattta                                                                                  >  1:2389979/1‑100 (MQ=255)
                     gCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTAtttt                                                                              >  1:1288404/1‑100 (MQ=255)
                     gCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTAtttt                                                                              >  1:2853743/1‑100 (MQ=255)
                      cAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTAtttt                                                                              >  2:1315185/1‑100 (MQ=255)
                      cAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTa                                                                             >  1:3517557/1‑100 (MQ=255)
                      cAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTa                                                                             >  1:2092194/1‑100 (MQ=255)
                      cAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTAACTAT‑AAAAAAACTTATTATTTATTTTa                                                                             >  2:3488373/1‑100 (MQ=255)
                        aTAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAg                                                                            <  2:1070922/100‑1 (MQ=255)
                        aTAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGt                                                                           >  1:72505/1‑100 (MQ=255)
                        aTAAAAACTATAACCAACATGAAGTGGCAAACTGGTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGt                                                                           >  1:3439238/1‑100 (MQ=255)
                         tAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGt                                                                           >  1:1457923/1‑100 (MQ=255)
                         tAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGtt                                                                          <  2:2401836/100‑1 (MQ=255)
                          aaaaaCTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACAGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGtt                                                                          >  2:2869109/1‑100 (MQ=255)
                          aaaaaCTATAACCAACATGAAGTGGCAAACTGTTTAACCCTATCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGttt                                                                         >  2:2270804/1‑100 (MQ=255)
                            aaaCTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGttttt                                                                       <  1:2510055/100‑1 (MQ=255)
                            aaaCTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGttttt                                                                       >  1:3780132/1‑100 (MQ=255)
                             aaCTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTa                                                                      <  2:2956680/100‑1 (MQ=255)
                              aCTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTAt                                                                     <  2:2116432/100‑1 (MQ=255)
                               cTATAACCAACATGAAGTGGTAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATc                                                                    <  1:2220359/100‑1 (MQ=255)
                                tataACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCa                                                                   <  2:757411/100‑1 (MQ=255)
                                  taACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGt                                                                 >  2:3920594/1‑100 (MQ=255)
                                    aCCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTgg                                                               >  1:364236/1‑100 (MQ=255)
                                     ccAACATGAAGTGGCAAACTGTTTAACCCTTTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTg                                                                <  2:395880/98‑1 (MQ=255)
                                     ccAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTgg                                                               >  2:453627/1‑100 (MQ=255)
                                     ccAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTggg                                                              <  2:2987250/100‑1 (MQ=255)
                                       aaCATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGt                                                             >  1:1304064/1‑100 (MQ=255)
                                       aaCATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGtt                                                            <  1:3578756/100‑1 (MQ=255)
                                          aTGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTAt                                                         <  2:4115086/100‑1 (MQ=255)
                                            gAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATtc                                                       >  2:313772/1‑100 (MQ=255)
                                            gAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATtc                                                       <  1:2662850/100‑1 (MQ=255)
                                            gAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATtc                                                       <  2:3608347/100‑1 (MQ=255)
                                            gAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATtc                                                       <  1:2317750/100‑1 (MQ=255)
                                            gAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGAATAAAAAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATtc                                                       <  2:1201926/100‑1 (MQ=255)
                                                tGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTctct                                                    <  1:611280/100‑1 (MQ=255)
                                                tGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCtt                                                   <  1:3157142/100‑1 (MQ=255)
                                                 ggCAACCTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGCTTGATTAT‑AAAAAAACTTATTATTTATCTTAGTTTTTATCAGTGGGTTGATTCGCTTc                                                  >  2:3618536/1‑100 (MQ=255)
                                                 ggCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCtt                                                   >  1:2795796/1‑100 (MQ=255)
                                                 ggCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCtt                                                   >  1:4046557/1‑99 (MQ=255)
                                                 ggCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTc                                                  >  1:3422201/1‑100 (MQ=255)
                                                  gCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTc                                                  <  2:3198213/99‑1 (MQ=255)
                                                  gCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCa                                                 >  1:4133763/1‑100 (MQ=255)
                                                  gCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCa                                                 >  2:636743/1‑100 (MQ=255)
                                                   cAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCAc                                                <  2:3179730/100‑1 (MQ=255)
                                                    aaaCTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACt                                               >  1:1153171/1‑100 (MQ=255)
                                                       cTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTc                                            <  1:4168375/100‑1 (MQ=255)
                                                       cTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTc                                            <  2:2250305/100‑1 (MQ=255)
                                                         gcttAACCCTGTCTTACCGTGAGGTTAATAGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCt                                           >  2:311264/3‑100 (MQ=255)
                                                         gTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTa                                          >  2:566326/1‑100 (MQ=255)
                                                         gTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTa                                          >  1:2191828/1‑100 (MQ=255)
                                                            tAACCCTGTCTTACCGTGATTTTAATCGTAATTAAATATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTAtt                                        <  2:1132721/100‑1 (MQ=255)
                                                            tAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTc                                       >  1:570133/1‑100 (MQ=255)
                                                            tAACCCTGTCTTACAGTTAGGTTAATCGTGATTTATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTc                                       <  2:3573802/100‑1 (MQ=255)
                                                              aCCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCtt                                     >  2:3303380/1‑100 (MQ=255)
                                                              aCCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCtt                                     <  2:2886752/100‑1 (MQ=255)
                                                              aCCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCtt                                     <  2:4014282/100‑1 (MQ=255)
                                                               cccTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCtt                                     >  2:2461816/1‑99 (MQ=255)
                                                               cccTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTa                                    <  2:1621655/100‑1 (MQ=255)
                                                               cccTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTa                                    <  1:1992887/100‑1 (MQ=255)
                                                                ccTGTCTTACCGTTTGGTTAATCGTGATTGATTAT‑GAAAGACCTTATGATTTATTGTTGTGTTTATCAGGGGGTTTATTCTCTTCACTGTCTATTCTTac                                   <  2:1366988/100‑1 (MQ=255)
                                                                ccTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTac                                   <  1:1637802/100‑1 (MQ=255)
                                                                     cTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTg                               <  1:2607449/100‑1 (MQ=255)
                                                                     cTTACCGAGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTATTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGc                              >  1:2916686/1‑100 (MQ=255)
                                                                      ttattgTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCa                             <  2:2637942/95‑1 (MQ=255)
                                                                       tACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTacac                                 <  2:518377/95‑1 (MQ=255)
                                                                        aCCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGATTTATTCTCTTCACTGTCTATTCTTACACTGCAg                            >  2:2707707/1‑100 (MQ=255)
                                                                        aCCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTTTTACACTGgaga                           >  2:3662953/1‑100 (MQ=255)
                                                                              aGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAAt                     <  2:1178231/100‑1 (MQ=255)
                                                                               ggTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAAt                     >  1:3566358/1‑100 (MQ=255)
                                                                                gtgaaTCGGCATTGAATAT‑AAAAATACTTACTATTTATTTAAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCCTACACTGCAGAAAAAATGa                   <  2:829349/97‑1 (MQ=255)
                                                                                 ttAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCATTGGATTTATTCTCTTTACATTGTATTTTTATAATGTATAAAatat                     >  2:655258/1‑97 (MQ=255)
                                                                                 ttAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTTCACTGCAGAAAAAATg                    <  2:1220873/98‑1 (MQ=255)
                                                                                 ttAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGat                  <  2:1085235/100‑1 (MQ=255)
                                                                                 ttAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGa                   >  1:1534946/1‑99 (MQ=255)
                                                                                 ttAATCGTGATTGATTAT‑AAAAAAACTTATTAATTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGat                  >  1:3958439/1‑100 (MQ=255)
                                                                                  tAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTAAACTGCAGAAAAAATGat                  >  2:414745/1‑100 (MQ=255)
                                                                                  tAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGata                 <  1:791486/100‑1 (MQ=255)
                                                                                    atagTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCATTGTCTATTCTTACACTGCAGAAAAAATGatat                <  2:216283/96‑1 (MQ=255)
                                                                                    aTCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGatata               >  1:400263/1‑100 (MQ=255)
                                                                                     tCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGatatat              >  2:1446967/1‑100 (MQ=255)
                                                                                       gTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATaa            >  1:1427830/1‑100 (MQ=255)
                                                                                       gTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGAGTTTATTGACTTCATTTTCTATTCTTACATTGCAGAAAAAATGatatct              >  2:2662059/1‑96 (MQ=255)
                                                                                         gattgattATAAAAAAAACTTATTATTTATTTTAGTTATTATCAGTGGGTTTATTCTATTTATTGACTATTTTTATATTGAATAAAAAATGATATATaa            >  2:3294345/1‑99 (MQ=255)
                                                                                         gattgattAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACg          <  2:1573156/100‑1 (MQ=255)
                                                                                           ttgattAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGaa        <  1:1815676/100‑1 (MQ=255)
                                                                                             gattAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTATATTCTTACACTGCAGAAAAAATGATATATAACGAAAc      >  1:346386/1‑100 (MQ=255)
                                                                                             gattAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGAGATTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATaa            >  2:2513904/1‑94 (MQ=255)
                                                                                              attAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACa     <  1:328329/100‑1 (MQ=255)
                                                                                                tataAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACAt    >  1:2959795/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGC  >  NC_000913/2989217‑2989414

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