Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,362,529 T→C 15.3% intergenic (‑176/+22) yjdQ ← / ← pheU pseudogene, P4‑like integrase remnant/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,362,5290TC15.3% 390.1 / 62.1 170intergenic (‑176/+22)yjdQ/pheUpseudogene, P4‑like integrase remnant/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base T (56/88);  new base C (18/8);  total (74/96)
Fisher's exact test for biased strand distribution p-value = 5.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-02

Only 100 of 171 total aligned reads displayed.

GCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAA  >  NC_000913/4362432‑4362628
                                                                                                 |                                                                                                   
gCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGACAAACGGACCTCCGTgg                                                                                                   >  2:3101551/1‑100 (MQ=255)
gCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTgg                                                                                                   >  2:1851981/1‑100 (MQ=255)
  tGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAg                                                                                                 <  2:1190609/100‑1 (MQ=255)
   ggTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAgg                                                                                                >  1:229772/1‑100 (MQ=255)
   ccgatcttgGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAgg                                                                                                <  1:3150367/93‑1 (MQ=255)
      gATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTcc                                                                                             >  2:1475262/1‑100 (MQ=255)
        tgtgGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGt                                                                                           <  1:3260203/100‑1 (MQ=255)
         gtgGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGtt                                                                                          >  2:2555108/1‑100 (MQ=255)
          tgGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAggagatcgg                                                                                         >  2:3150346/1‑93 (MQ=255)
           ggggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACTGACCTCCACGGAGGTCCGtttt                                                                                        <  2:2289918/100‑2 (MQ=255)
           ggggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTta                                                                                        <  1:2964077/100‑1 (MQ=255)
           ggggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTta                                                                                        <  1:3232257/100‑1 (MQ=255)
           ggggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTta                                                                                        >  2:1435410/1‑100 (MQ=255)
           ggggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTta                                                                                        >  2:3747817/1‑100 (MQ=255)
           ggggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTta                                                                                        >  2:3599801/1‑100 (MQ=255)
           ggggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTta                                                                                        <  2:1903950/100‑1 (MQ=255)
           ggggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGtttt                                                                                        <  1:3677664/100‑2 (MQ=255)
             ggAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTtata                                                                                      <  1:1161122/100‑1 (MQ=255)
                gaATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGa                                                                                   >  2:669999/1‑100 (MQ=255)
                gaATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGtttttcgtt                                                                                   <  1:3582296/100‑7 (MQ=38)
                  aTCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAAt                                                                                 >  1:2550700/1‑100 (MQ=255)
                  aTCTGCTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAACTAGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTacatac                                                                                   >  1:1175825/1‑96 (MQ=255)
                   tCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCACGGAGGTCCGtttttcgtttca                                                                                <  2:591673/100‑10 (MQ=37)
                     tGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTg                                                                              >  1:1523516/1‑100 (MQ=255)
                       gTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGt                                                                            <  2:1549532/100‑1 (MQ=255)
                       gTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGt                                                                            >  1:1853494/1‑100 (MQ=255)
                            gTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCg                                                                       >  1:237938/1‑100 (MQ=255)
                                 gTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTc                                                                  >  1:204997/1‑100 (MQ=255)
                                   aGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTTAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCgg                                                                <  2:2845168/100‑1 (MQ=255)
                                   aGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCgg                                                                <  2:229772/100‑1 (MQ=255)
                                   aGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCgg                                                                <  1:1056406/100‑1 (MQ=255)
                                     aaTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGaa                                                              <  1:2653947/100‑1 (MQ=255)
                                     aaTTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGaa                                                              >  1:2131747/1‑100 (MQ=255)
                                          aTTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGaa                                                         <  2:697074/100‑1 (MQ=255)
                                           tttGGAGGCAGGACGCTTAAAGCGTGGCGTACTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCACGGACTCGTAATCGa                                                          <  2:1133950/98‑1 (MQ=255)
                                           tttGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAGCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAAc                                                        <  2:1320158/100‑1 (MQ=255)
                                           tttGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAAc                                                        <  1:796411/100‑1 (MQ=255)
                                            ttGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAAcc                                                       <  2:1258777/100‑1 (MQ=255)
                                             tGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCa                                                      <  1:1243801/100‑1 (MQ=255)
                                               gAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAg                                                    <  1:2468630/100‑1 (MQ=255)
                                               gAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAg                                                    <  1:848691/100‑1 (MQ=255)
                                               gAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCGCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCaan                                                    <  2:2723978/100‑2 (MQ=255)
                                                aGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAgg                                                   <  1:2315141/100‑1 (MQ=255)
                                                   cAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGaca                                                <  1:1548865/100‑1 (MQ=255)
                                                    aGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGacac                                               <  2:2673958/100‑1 (MQ=255)
                                                      aaCGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACgg                                             <  2:2391304/100‑1 (MQ=255)
                                                       aCGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACggg                                            <  2:2255347/100‑1 (MQ=255)
                                                       aCGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTCGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACggg                                            <  2:2517557/100‑1 (MQ=255)
                                                          cTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGAt                                         >  2:3453198/1‑100 (MQ=255)
                                                          cTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGAt                                         >  1:1200670/1‑100 (MQ=255)
                                                            tagatggTGGCGTCCTGGAGCGAAGAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACggg                                            <  2:3565925/89‑1 (MQ=37)
                                                              aaTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTc                                     <  1:650087/100‑1 (MQ=255)
                                                              aaTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTc                                     <  1:3866702/100‑1 (MQ=255)
                                                              aaTCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTc                                     <  1:1239076/100‑1 (MQ=255)
                                                                tCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCaa                                   <  1:718200/100‑1 (MQ=255)
                                                                 cGTGGCGTCCTGAAACGAAAGACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAAt                                  <  2:3900999/100‑1 (MQ=255)
                                                                 cGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCaa                                   <  1:2905219/99‑1 (MQ=255)
                                                                  gTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATc                                 >  1:3448436/1‑100 (MQ=255)
                                                                  gTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATc                                 >  2:1620356/1‑100 (MQ=255)
                                                                      cGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCt                             <  1:1780686/100‑1 (MQ=255)
                                                                         ccTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGct                          >  1:2352447/1‑100 (MQ=255)
                                                                          cTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGctc                         >  1:2856021/1‑100 (MQ=255)
                                                                           tGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGctct                        <  2:2988087/100‑1 (MQ=255)
                                                                           tGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGctct                        >  2:98643/1‑100 (MQ=255)
                                                                            gAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTa                       >  1:1430019/1‑100 (MQ=255)
                                                                             aaaCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTAc                      <  1:2739438/100‑1 (MQ=255)
                                                                             aaaCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTAc                      >  1:1114971/1‑100 (MQ=255)
                                                                              aaCGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTAcc                     <  2:3671731/100‑1 (MQ=255)
                                                                               catataaaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCg                    >  1:1291814/6‑100 (MQ=18)
                                                                                 gAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGAc                  <  1:1793528/100‑1 (MQ=255)
                                                                                 gAAAAAAGGACATCCGTGTAGGTCCGTTTATATGAATTTGGGTCCCGGGCTCGGAATCGATCCAAGGACACGGGGTTTGGCAAGCCCCGGcgttac                      >  2:727303/1‑91 (MQ=255)
                                                                                  ataaaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACt                 >  2:3243684/3‑100 (MQ=255)
                                                                                  aaaaaCGGACCTCCGTGGAGGTCCGTTTATACGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACt                 <  2:2338114/100‑1 (MQ=255)
                                                                                   taaaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTg                >  2:2625116/2‑100 (MQ=255)
                                                                                   taaaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTg                >  1:2797048/2‑100 (MQ=255)
                                                                                   aatgcGGCCCTCCGTTGAGCTCTCTTTATATGTATATGTCGGGCGCACTCCGACTCCAACCTAGCACACGGCGATTTTCAATCCCCTGCTTTAc                      <  2:3030091/90‑1 (MQ=2)
                                                                                   aaaaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTg                <  1:1989798/100‑1 (MQ=255)
                                                                                    aaaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACTCGGGGATTTTCAATCCCCTGCTCTACCGACTGa               <  2:2365400/100‑1 (MQ=255)
                                                                                    aaaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTGTCAATCCCCTGCTCTACCCAGTGa               >  2:2852488/1‑100 (MQ=255)
                                                                                    aaaCGGACCTCCACGGAGGTCTGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa               >  2:934999/1‑100 (MQ=255)
                                                                                    aaaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGa               >  1:2248826/1‑100 (MQ=255)
                                                                                    aaaCGGACCTCCACGGAGGGCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCTCTGCTCTACCGACTGa               >  1:581781/1‑100 (MQ=18)
                                                                                     aaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAg              >  2:3040241/1‑100 (MQ=255)
                                                                                     aaCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAg              <  1:1609402/100‑1 (MQ=255)
                                                                                     aaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAg              >  2:1091948/1‑100 (MQ=255)
                                                                                     aaCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAg              >  2:1472431/1‑100 (MQ=255)
                                                                                      aCGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGc             <  1:82576/100‑1 (MQ=255)
                                                                                      aCGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACAAGGGGATTTTCAATCCCCTGCTCTATCGAATGAg              >  2:2484642/1‑99 (MQ=18)
                                                                                       cGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCt            >  1:1099746/1‑100 (MQ=255)
                                                                                       cGGACCTCCACGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCt            >  2:1109657/1‑100 (MQ=255)
                                                                                           ccTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATcc        >  2:420020/1‑100 (MQ=255)
                                                                                           ccTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATcc        >  1:870482/1‑100 (MQ=255)
                                                                                            cTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCg       <  1:3453221/100‑1 (MQ=255)
                                                                                            cTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCg       >  1:3064404/1‑100 (MQ=255)
                                                                                             tCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCgg      <  1:3599829/100‑1 (MQ=255)
                                                                                              ccGTGGATGTCTGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTCTCAATCCCCTGCTCTACCGACTGAGCTATCCggg     >  2:2077210/1‑100 (MQ=255)
                                                                                              ccGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGCATCGAACCAAGGACACGGGGATTTTCAATGCCCTGCTCTACCGACTGAGCTATCCggg     >  2:874220/1‑100 (MQ=255)
                                                                                              ccGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCggg     <  1:3217349/100‑1 (MQ=255)
                                                                                               cGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGc    <  1:3677739/100‑1 (MQ=255)
                                                                                                 tGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCaa  >  2:919992/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GCTGGTGATGTGGGGAGAATCTGGTTGAGTTCGGTAGAATTGATTTGGAGGCAGAACGCTTAAATCGTGGCGTCCTGAAACGAAAAACGGACCTCCGTGGAGGTCCGTTTATATGAATTTGGTGCCCGGACTCGGAATCGAACCAAGGACACGGGGATTTTCAATCCCCTGCTCTACCGACTGAGCTATCCGGGCAA  >  NC_000913/4362432‑4362628

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: