Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,595 T→G 11.0% intergenic (+61/‑135) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,5950TG11.0% 369.2 / 15.1 133intergenic (+61/‑135)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (14/101);  new base G (15/2);  total (29/104)
Fisher's exact test for biased strand distribution p-value = 3.86e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 142 total aligned reads displayed.

AAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATG  >  NC_000913/164446‑164742
                                                                                                                                                      |                                                                                                                                                   
acgataTGAAGGGGGTTTATCCCAAACATGTCTGGCGGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAATTTTTTCTTATGTTTGTTTTTCTCCCTCTCCCTGTGGTTTTGGGCGGGGTTTTGGGCATCAGCGCGCACGtt                                                                                                                                                     <  2:2343872/145‑1 (MQ=255)
   agaTGAATGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAATGTATTCTTAAGTTTGTTTTTCTCCCTCTCCCTTTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcac                                                                                                                                                  <  2:3040885/150‑1 (MQ=38)
   agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAATTTATTCGTAAGGTTGGTTTTCGCCCTCTCCCTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcac                                                                                                                                                  <  2:3370062/150‑1 (MQ=255)
     aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaccc                                                                                                                                                >  2:2643346/1‑150 (MQ=255)
       gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcgccccc                                                                                                                                              >  2:1906669/1‑148 (MQ=255)
        aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                             <  2:3540528/150‑1 (MQ=255)
              cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                                       <  2:3908165/150‑1 (MQ=255)
               gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCGCCCACACCCAAAc                                                                                                                                      >  1:331383/1‑150 (MQ=38)
               gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGGCCCCCACCCCaaa                                                                                                                                      >  1:2676431/1‑149 (MQ=38)
                ttATCCAAAACTTGTCTGGCGGGACGACCCGGAAAATCCTGCACCGACGCGAGGGACGAAAATTTTTTCGTACGGTTGTTTTTCTCCCCCCCCGGGGGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                                     <  2:5023776/150‑1 (MQ=255)
                    cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                                 >  1:4799987/1‑150 (MQ=255)
                     ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCctcta                                                                                                                                >  2:491047/1‑149 (MQ=255)
                     ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCACAcccacac                                                                                                                                >  2:5325887/1‑146 (MQ=38)
                      cAAACCTGGCTGGCCGGACGACCCGGCAAATACTGCCCCGACGCGCCGGACGAAAATGTATGGGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                               <  2:3441073/150‑1 (MQ=37)
                      cAAACATGTCTGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                               <  2:4872175/150‑1 (MQ=255)
                      cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                               <  2:3441600/150‑1 (MQ=255)
                       aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGCCCCCCACCCCAACCCACAccc                                                                                                                              >  1:358372/1‑150 (MQ=37)
                         aCATGTCTGGCGGGACGACCCGGAAATTACTGCCCCGACGCGACGGACTAAAATGTATTCTTACGGTTGGTTTTCTCCCCCTCCCTGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTc                                                                                                                            <  2:5657769/150‑1 (MQ=255)
                           tttagtGGCCGGCCCCCCGGGCAAATCCTGCCCCGACGCCCCGTACGAAAAAGTTTTCTTATGTTGGGTTTTCCCCCCCGCCCGGGGGGTGTGGGCGGGGTTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                          <  2:5050563/145‑1 (MQ=255)
                           aTGTCTGGCCGGACGCCCCGGCAAATACTGCCCCGACGCGACGGACGACAAAGTAGTCGTCCGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                          <  2:4887186/150‑1 (MQ=37)
                            tgtggggccggccGCCCCGGAAAATACTGCCCCGCGGCGCCGTACTAAAATTTTTTGGTTTGGTTGTTTTTCACCCACTCCGTGGGGGTGTGGGGGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                         <  1:5851939/145‑1 (MQ=255)
                            tGTCTGGCGGGACGACCCGGCAAATCCGGCACCGCCGCGACGGACTAAAATTTTTTCGTTTGGTTGGTTTTCCCCCTCTCCGTGTGGGTTTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                         <  2:1711122/150‑1 (MQ=255)
                            tGTCGGGCCGGCCGACCCGGAAAATCCTGCCCCGACGCGACGGACTAAAAAGTATTCTTACGGTTGGTTTTCACCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                         <  2:332692/150‑1 (MQ=255)
                            ggtgtGGCCGGACGACCGGGAAAATACTGCACCGCGGCGACGGACTAAAATGTATTCTTATGTTTGGTTTTCTCCCTCTCCGTGTGGTAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                         <  2:1711019/146‑1 (MQ=255)
                                tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACAACCACAAGGgag                                                                                                                     >  2:4504978/1‑150 (MQ=38)
                                tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACAACCAAACCCACTCCCACAAggggg                                                                                                                     >  2:4505283/1‑148 (MQ=38)
                                 ggCGGGACGCCCCGCCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGTTTGGTTTTCGCCCTCACCGGGGGGGCGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                    <  2:4156635/150‑1 (MQ=255)
                                 ggCCGGACGACCCGGCAATTACTGCACGGACGCGACGGACGAAAAAGTTTTGGTAAGGTTGGTTTCCTCCCTCGCCCGGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                                    <  2:4814256/150‑1 (MQ=255)
                                  gCCGGACGACCCGGCAAATACTGCACCGACGCGACGTACGAAATAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                   <  2:4864509/150‑1 (MQ=255)
                                  gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGCCCCGCACCCCAAACCACGCCCGCAAggggggg                                                                                                                   >  1:1478369/1‑146 (MQ=37)
                                  gCCGGACCCCCCGCAAAATCCTGCACCCACGCGACGGACTAAAAATTTTTCGTATGTTTGGTTTTCCCCCTCTCCCTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                   <  2:5647137/150‑1 (MQ=255)
                                  cccGGACCCCCCGGAAAATACTGCACCGCCGCGACGCACGAAAAATTATTGTTACGGTTGTTTTTCGCCCACCCCCGGTGGGCGAGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                   <  2:2216801/149‑1 (MQ=255)
                                   cgggACTACCGGGAAAATACTGCACGTCCGCTAGGTACTAAAATTTTTTGTTATGGTTGGTTTTCCCCCTCTCCGTGTGGGTTTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                                  <  1:2888469/148‑1 (MQ=255)
                                   ccGGACCACCGGGAAAATACTGCACCCACGCGAGGGACCAAAATGTTTGCGTACGGTGGGTTTTCACCCACCCCCTGTGGGAGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                                  <  2:1597878/150‑1 (MQ=255)
                                    cgtacgacCCCTAAAATATTCCACCTACGCGCCTCACTAAAATTTTTTCTTATGTTTGTTTTTGCCCCTCTCCTTGTGGTTGTGGGCGGGGGTTTGGGCATCAGCGCGCACGTTCACACTCACCCTAACACTCTCACTCAAGGGAGAggg                                                                                                                 <  2:5054183/147‑1 (MQ=255)
                                     ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCACCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCACACCCAAAACCACTCCCACAAGGGAGAgggg                                                                                                                >  2:1498538/1‑150 (MQ=38)
                                     cggcccccccGCAAATAACTGCCCCGAGGCGAGGGACTAAAATTTTTTCGTACGGTTGGTTTACTCCCACTCCGTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                                 <  1:3978622/144‑1 (MQ=255)
                                      gacgccCCGGAAAATACGGCCCCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCGCGCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                               <  2:729266/150‑1 (MQ=38)
                                      gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCCCCCCAACCCGCGCCCGCCACGGGGAggggg                                                                                                               >  2:5470533/1‑149 (MQ=37)
                                      gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCGAACCCACACCCGCACGGGGGAggggg                                                                                                               >  2:5470541/1‑149 (MQ=37)
                                      gaccacCCGCCAAATACTGCCCCGCCGCGACGGACGAAACAGTATGCGTAAGGTGGGTTTCCTCCCCCCCCCTGTGGGAGAGGGCCGGGTTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGTGAGAGGGGa                                                                                                               <  2:729221/150‑1 (MQ=255)
                                           cccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAAACCACAACCGCAAAGGGGAGGGGACCGaa                                                                                                          >  1:3510512/1‑149 (MQ=25)
                                            ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAAACCAAAACCACAAGGGGGAGGGGccccaac                                                                                                         >  2:3833282/1‑146 (MQ=25)
                                             cGGCAAATACTGCACCGACGCGACGGACGACAAAGTATTCGTCAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCGCACCCCAACCCGCGCCCCCCCGGGCGCGGGGAccccgcg                                                                                                        >  2:2491752/1‑145 (MQ=21)
                                             cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCGCACCCAAACCCACACCCACAACGGAGCGGGGGCCGAACg                                                                                                        >  1:5740558/1‑150 (MQ=25)
                                               gCAAATACTGCCCCGACGCGACGGACCAAAAAGTTTTCTTAAGGTGGGTTTTCACCCTCTCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAg                                                                                                      <  1:4356040/150‑1 (MQ=255)
                                                aaaaTCCTCCACCGACGCGACGGACTAAAATTTATTCTTCTGTTTGGTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                                     <  1:2872994/149‑1 (MQ=255)
                                                 aaaTACTTCACCGAGGCGCCTGACTAAAATTTTTTCTTATGTTTGTTTTTCTCCATACCCCTGTGGGAGTGGGCCGGGTTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                                    <  2:200355/150‑1 (MQ=255)
                                                       cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCGAAACCACACCCTCCCGGGAGCGGGGGCCGGTCGCGCCCCAAt                                                                                               >  1:5880434/1‑150 (MQ=25)
                                                         gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACACCCAAAACCAAAACCACAAAGGGGAGGGGGCCCCACGAGcacacaaaa                                                                                             >  2:4702842/1‑145 (MQ=21)
                                                          cACCGCCGCGACGGACGAAAATGTTTTCTTATGGTTGGTTTTCGCCCTCTCCGGGTGGGAGTGGGCGGGGGGGTGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                            <  1:1362674/150‑1 (MQ=255)
                                                          cACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                            <  1:853925/150‑1 (MQ=255)
                                                            ccgccgCGACGCACTAAAAAGTTTTCGTACGGTGGGTTTACGCCCTCCCCCTGGGGGTGCGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                          <  2:2246253/150‑1 (MQ=255)
                                                            ccGACGCGACGTACGAAAAAGTATGCGTACGTTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                          <  1:3729605/150‑1 (MQ=38)
                                                            ccGACGCGACGGACTAAAATGTTTTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                          <  2:3654856/150‑1 (MQ=38)
                                                            ccGACGCGACGGACGAAAAATTATTCTTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGTGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                          <  2:1459364/150‑1 (MQ=38)
                                                            ccGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                          <  2:4451034/150‑1 (MQ=255)
                                                             ggccgcgacgGACCAAATTGTAGTCGTACGGTTGGTTTTCCCCCTCTCCGTGGGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                                         <  1:4757748/147‑1 (MQ=255)
                                                                      ggacgAAAATGTATTCGTATGTTTGGTTTCCGCCCTCTCCCTGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGTGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                                <  2:4235759/150‑1 (MQ=38)
                                                                          gAACAAGTATTCGTACGGGTGGTTTTCTCCCTCGCCCTGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                                            <  2:1967817/150‑1 (MQ=255)
                                                                          gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCCAACCAACACACACAAAGGGAAGGGGGCCCAGCCAGCCCAACATTTTAGAGATAtcaaataa                                                                            >  1:4189009/1‑144 (MQ=12)
                                                                              aattATTCGGTAGGTTGGTTTGCCCCCTCACCCTGGGGTTGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                                        <  2:4111279/147‑1 (MQ=255)
                                                                                gTATTGGTTGGGTTGGTGTTCTCCCTCGCCGGGGGGGTGTGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                      <  2:609880/150‑1 (MQ=255)
                                                                                gTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                      <  2:3576614/150‑1 (MQ=255)
                                                                                    tggTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCTCCCTCACGCTAACCCTCTCCCTGACGGGAGAGGGGCCCGATGGAGCACAACTTTTGAGAGATCTCGTCTTCTGTGTTGTaa                                                                  <  2:3976378/148‑1 (MQ=25)
                                                                                      gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCTGGGGTGAGGGCATCAGCGCGCACGTTCAACCGCACCCCAAACCGCGACCCCAAGGGAGAGGGGGCCGCGCGCGCCCAAAATTTGTGAGCTATCGTCTTCTGTCGTGTGACa                                                                >  2:3495373/1‑150 (MQ=16)
                                                                                         cgtttggttTCCCCCCCCCCCCGGGGGGAGTGGGCCGGGGTGCGGGCATCCGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                             <  1:1633985/147‑1 (MQ=255)
                                                                                         cggtgggttTTCGCCCGCGCCGTGTGGGCGCGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                             <  2:1704955/149‑1 (MQ=255)
                                                                                           gttgtttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                           <  1:4505283/150‑1 (MQ=255)
                                                                                           gttggttTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                           <  1:4505439/150‑1 (MQ=255)
                                                                                             tggttTACTCCCCCCCCCTGGGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                         <  1:1868395/150‑1 (MQ=255)
                                                                                              ggttTTCGCCATCTCCGGGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                        <  1:2760826/150‑1 (MQ=255)
                                                                                                   tctccctctcccGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                                   <  2:1446265/150‑1 (MQ=255)
                                                                                                     cccctccccttGTGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                 <  1:512413/149‑1 (MQ=255)
                                                                                                     cccctcaccccGGGGGAGCGGCCCGGGGCGCGGCCATCAGCGCGCACCTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTGTGAGAGAGATCTTCCTCTGTCTTGTAACAGAAGAACAGAAAATc                                                 <  2:243296/149‑1 (MQ=255)
                                                                                                      ccctctccctGTGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTATTGTAACAGAAGAACAGAAAATCg                                                <  2:3707100/150‑1 (MQ=255)
                                                                                                       cctctccgtGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                               <  1:3277824/150‑1 (MQ=255)
                                                                                                        ctcgccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                              <  2:3592649/150‑1 (MQ=255)
                                                                                                        ctccccctGTGGGAGAGGCCCGGGGTGAGGGCATCAGCGCGCACGTACACCCTCACCCTAACCCACTCCCTCAAGGGAGAGGGGACCGAACGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                              <  2:5482463/150‑1 (MQ=255)
                                                                                                             cctGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                                         <  2:5817471/150‑1 (MQ=255)
                                                                                                                 tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACACCCGCAAGGGGGTGGGGGCCCATCGAGCACAAAATTTTAGATATATATTAATATGTCTTTTAACAGAAGAAAAGAAAAAAGGGCGTTTtcgc                                     >  2:3007351/1‑150 (MQ=255)
                                                                                                                 ggggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                     <  1:1386496/149‑1 (MQ=255)
                                                                                                                   ggAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                   <  2:2375443/150‑1 (MQ=255)
                                                                                                                   ggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTCACCCAAAACCAAGACCAAAAGGGGGGGGGGGCCCGGCGCGGACAAAATTTGGGAGAAATGTTGTTGTGTGTTTTTACAGAAGAACAGAAAATCGGGCgttttggccg                                   >  2:3368834/1‑144 (MQ=255)
                                                                                                                        gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                              <  2:4376385/150‑1 (MQ=255)
                                                                                                                         ggCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                             <  1:1906669/150‑1 (MQ=255)
                                                                                                                           cgggggggTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACATAAAATCGGGCTTTTGCGCCTGAATATTg                           <  1:734269/148‑1 (MQ=255)
                                                                                                                            cGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                          <  1:3324446/150‑1 (MQ=255)
                                                                                                                             ggggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                         <  1:5550300/150‑1 (MQ=255)
                                                                                                                             ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                         <  1:1950152/150‑1 (MQ=255)
                                                                                                                             ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                         <  1:5995372/150‑1 (MQ=255)
                                                                                                                                   aGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa                   <  1:1364337/150‑1 (MQ=255)
                                                                                                                                   aGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCACTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa                   <  2:4799987/150‑1 (MQ=255)
                                                                                                                                    gggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCGGAATATTGCGGAGaaaa                  <  2:2431823/150‑1 (MQ=255)
                                                                                                                                    gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa                  <  2:4623265/150‑1 (MQ=255)
                                                                                                                                      gcctCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg                <  2:2565258/147‑1 (MQ=255)
                                                                                                                                        aTCAGCGCGCACGTTCACCCACACCCAAAACCACAACCAAAAGGGAGAGGGGACCCATCGAGCACAAAATTTGAGAGAAATAATATTATGTCTTGTAAAAGAAGAACAGAAAAAAGGGCGTTTTCGCCGGAATAATTCGGAGAAAAAACa              >  2:3498292/1‑150 (MQ=255)
                                                                                                                                         tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAACGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCAt             <  2:5250969/150‑1 (MQ=255)
                                                                                                                                              gcgcACGTTCACCCTCACCCTAACCCTCTCCCTCTAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg        <  1:3694863/150‑1 (MQ=255)
                                                                                                                                                    gTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTATTGTAACACAAGAACAGAAAATCGGGCTTTTGCGCCAGAATATTGCGGAGAAAAAGCATGGCCGGGAATg  >  2:4903144/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
AAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATG  >  NC_000913/164446‑164742

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