Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 22.6% intergenic (+47/‑149) 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,5810AC22.6% 203.9 / 34.3 142intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (46/62);  new base C (2/32);  total (48/94)
Fisher's exact test for biased strand distribution p-value = 4.58e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730
                                                                                                                                                     |                                                                                                                                                      
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                       >  2:2617888/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                       >  2:2614896/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGACTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                       >  2:4011043/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGAGATATCCCAAACATGTCTGGCCGGACGACACGGCAAATACTGCACCGACGAGACGGACGAAAAAGAATTCGTAAGGTTGTTATACTCCCACTCCCTGTGGGAGAGGGCCTGGGAGAGGGCaa                                                                                                                                                       >  1:3701325/1‑149 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCCCCCCGGTCGGGCCGGACGACCCGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                     <  2:6098997/150‑1 (MQ=38)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                     <  2:6627968/150‑1 (MQ=255)
     tggTGCAAAAATATATGAATGGGCGTTATCCCAAACTTGTGTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCTTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGTGGGCATCAgc                                                                                                                                                   <  2:3990575/149‑1 (MQ=37)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgc                                                                                                                                                 >  1:1100840/1‑150 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                 >  1:2872497/1‑150 (MQ=255)
            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                            >  1:5104648/1‑150 (MQ=255)
             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgg                                                                                                                                           >  1:6518762/1‑149 (MQ=255)
               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGttg                                                                                                                                         >  2:732680/1‑149 (MQ=255)
                gtgaTGAACGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtca                                                                                                                                        <  2:1038553/148‑1 (MQ=255)
                gagaTGAATGGGCGTTATCCCAAACAGGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTTAGGTTGGTTTTGTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtca                                                                                                                                        <  2:4448866/150‑1 (MQ=38)
                gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtca                                                                                                                                        <  2:3142659/150‑1 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGtggcccc                                                                                                                                     >  1:241871/1‑144 (MQ=255)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTGTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCCGCGCGCCCGtgggccc                                                                                                                                     >  2:6369418/1‑144 (MQ=38)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcccc                                                                                                                                    >  2:6397875/1‑143 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccac                                                                                                                                   >  2:4242859/1‑148 (MQ=255)
                      aaaGGGCGTTATCCCAAACATGTCTGGCCGGATGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccaca                                                                                                                                  >  2:4924802/1‑150 (MQ=255)
                       atgggCGTTATCCCAAACATGTCTGGCCGTACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcac                                                                                                                                 <  2:5638487/148‑1 (MQ=255)
                         gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccaccccc                                                                                                                               >  2:6153194/1‑150 (MQ=255)
                            cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                            <  2:2812880/150‑1 (MQ=255)
                              ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCaaaac                                                                                                                          >  2:3431119/1‑148 (MQ=255)
                                 tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCGAAcccac                                                                                                                       >  1:6779710/1‑148 (MQ=255)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCCCCCCAACCcact                                                                                                                      >  2:4911979/1‑150 (MQ=38)
                                         ttttCTGGCCGGACGACCGGGCAAATCCTGCACCTACGCTAGGTACGAAAATTTTTTGGTATGTTTGTTTTTCTCCCTGTCCCTGTGGGTGTGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                               <  2:2209790/147‑1 (MQ=255)
                                         aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCGCTCCCGCaa                                                                                                               >  2:2804900/1‑150 (MQ=38)
                                          ggTCGGGCCGGACGACCCGGCAACTCCGGCACCGCCGCGACGGCCGAAAAAGTATTGGTACGGTGGGTTTTCTCCCTCGCCGTGTGGGAGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                              <  2:3813977/149‑1 (MQ=255)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                        <  2:5051218/150‑1 (MQ=255)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCAAACACCCACCCACaaacggcag                                                                                                        >  1:1295927/1‑143 (MQ=25)
                                                gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCGAACCCGCTCCCCCCCGGGagag                                                                                                        >  2:4297485/1‑150 (MQ=37)
                                                cccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                        <  2:5613542/149‑1 (MQ=255)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGGGCCCCCGCCCCCAAACCCCCACCCCCACGGGCGCggg                                                                                                      >  1:6681851/1‑150 (MQ=25)
                                                  cGGACGACCCGGAAAACACTCCACCGACGCGACGGACTAAAATGTATTCTTAAGGTTGTTTCTCTCCCTCTCCCTGGGGTAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                      <  1:4964440/150‑1 (MQ=255)
                                                    gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCACCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                    <  2:3765842/150‑1 (MQ=255)
                                                    gacgacCCGGAAAATACTGCACCGACGCGAGGGACTAAAATTTATTCGTATGTTTGTTTTACTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                    <  1:1524928/150‑1 (MQ=255)
                                                     ccgccccGGCAAATACGCCCCCCCCGCGACGCCCGACACAGTACTCGTAAGGTGGGTTGTCGCCCCCCCCCGGGGGGCGCGGGCCGGGGTGCGGGCCTCAGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                   <  2:6398420/146‑1 (MQ=255)
                                                          ccgcaaaaTCCTGCACCGCCGCGACGGACTAAAAAGTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGGGGGCGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCCATc                                                                                              <  1:4381672/145‑1 (MQ=255)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCAAACCCACACCCACCAGGGCGAggggggccatc                                                                                              >  2:3113986/1‑143 (MQ=25)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGGACCCCCACCCCCACCCACGCCCGACAGGGGGGGGGGGCCgcgc                                                                                              >  1:912170/1‑147 (MQ=21)
                                                           cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGAACCGCACCCCAACCCTCAACCGCCAGGGGGGGGGGACCCAtag                                                                                             >  1:6177423/1‑148 (MQ=25)
                                                                     tGCACCGACGCGATGGACGAAAAAGTATACGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCGCCCAAACCCAAAAACACTCCCACAAGTGGGGGGGGACCCCTCGAGccccaaat                                                                                   >  2:368846/1‑143 (MQ=16)
                                                                     tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCCCCCCAACCCCCGCCCGCAAGGGTGAGGGGGCCCAGGCAGCACAAAtt                                                                                   >  2:556607/1‑150 (MQ=25)
                                                                       ccccGACGCGCGGGCCGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCGGTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                 <  1:4176319/148‑1 (MQ=255)
                                                                        cccGACGCGCGGGACTAAAATGTTTCGGTAGGGTGGGTTTTCCCCCTCCCCGGGGGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                <  2:1225123/149‑1 (MQ=255)
                                                                         ccgccgCGCCGGCCGAAAATTTATTCGTATGGTTGTTTTTCTCCCTCGCCCTGTGGGAGCGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                               <  2:2708679/150‑1 (MQ=255)
                                                                         ccgccgCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTCCCCCCCCTCCCGGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                               <  2:2708689/150‑1 (MQ=255)
                                                                         ccGACGCGCGGGACTAAAATTTATTGGTTGGGTTGGTTTTCTCCCTCCCCCGGGGGGTGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                               <  2:5281446/150‑1 (MQ=255)
                                                                         ccGACGCGCCGGAGGAAACAGTTGTCGTAGGGTTGTTTTTCTCCCGCTCCCTGTGGTTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                               <  2:3447454/150‑1 (MQ=255)
                                                                         ccGACGCGAGGGACGAAAAATTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                               <  2:2708831/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACTAAAATTTTTTCGTATGGTTGGTTTTCGCCCTCTCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                               <  1:6621754/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                               <  2:1092879/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCGCCCTCTCCCGGTGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                               <  2:907964/150‑1 (MQ=255)
                                                                                      cgAAAACGTATCGGTCCGGTTGGTTTTCTCCCTCGCCGTGTGGGAGCGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCCCCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                                  <  1:3102623/150‑1 (MQ=38)
                                                                                       gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCAACCCCCAAACCACACCCCAAAGGGAGAGGGGACCCAACGAGCGCCAAATTTGATAGaaataataaac                                                                 >  2:1119815/1‑143 (MQ=16)
                                                                                         aaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCAACCACACCCCCACCCACACCCACCCCGGCGGGGGGCCCCAGCGAGCCCAAAAATTGATAGCTAtcttcatagg                                                               >  1:6134513/1‑145 (MQ=16)
                                                                                           atgTATTCGTACGGTTGGTTTCCCCCCTCCCCCGGGGGGTGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                             <  1:3637765/148‑1 (MQ=255)
                                                                                             tttttCTTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                           <  1:3431119/147‑1 (MQ=38)
                                                                                             gTATTCGTAAGGTTGGTTTTCCCCCGCGCCCGGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                           <  1:4813732/150‑1 (MQ=255)
                                                                                              tATTCGTATGGTTGGTTTTCTCCCACTCCCGGGGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                          <  1:1393055/150‑1 (MQ=255)
                                                                                               aTTCGTACGTTTGTTTTTCCCCCCCTCCCTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                         <  1:5298138/150‑1 (MQ=255)
                                                                                                  cGTACGGTGTGTTTTCACCGTCTCCCTGTGTTTGCGGGCCGGGTTGTGGGCCTCAGCGCGCTCGTTCACCCTCACCCTATCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTacc                                                      <  2:2925239/150‑3 (MQ=255)
                                                                                                    tAAGGTTGGTTTTCTCCCTCTCCGTGTGGGTGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGTGGTCCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                    <  2:2057805/150‑1 (MQ=255)
                                                                                                    tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGAACCACACCCAAAACCAAAACCACCAGGGGGAGGGGGCCGAAAGAGCACACAATTTTAGAGATATCATATTCGGTGATGTAACag                                                    >  2:6753391/1‑150 (MQ=12)
                                                                                                     tcggttggttTTCCCCCTCTCCGTGTGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                   <  1:5773294/148‑1 (MQ=255)
                                                                                                     ctggttggttTTCACCCTCGCCGTGTGGGCGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                   <  1:5773315/148‑1 (MQ=255)
                                                                                                        gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                <  2:6710808/150‑1 (MQ=255)
                                                                                                         ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGACATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                               <  1:6449257/150‑1 (MQ=255)
                                                                                                          tggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                              <  2:3359819/150‑1 (MQ=255)
                                                                                                           ggttTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                             <  1:3191513/150‑1 (MQ=255)
                                                                                                             ttgcccccctctccgtGGGGGTGTGGGCGGGGTTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGa                                          <  1:2516105/144‑1 (MQ=255)
                                                                                                                tctccctctcccTGTGGGAGAGGGCCGGGGTGAGGGCATCGGCGCGCACGTGCACCCCCACCCCAAACCACGACCCCCAGGGAGAGGGGGCCGCTCGCGCACAAAATTTGGGTGATATGGTATTCGGTCTTGTAACGGAAGAACAGagaa                                       >  2:1568859/1‑150 (MQ=16)
                                                                                                                  tccctctccctGGGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                     <  1:5473923/150‑1 (MQ=255)
                                                                                                                    cctctccctGTGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  1:5026851/150‑1 (MQ=255)
                                                                                                                    cctctcccaCAGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                   <  2:4983338/71‑1 (MQ=25)
                                                                                                                      tctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCTCTCCCGCAAGGGCGAGGGGCCCGATCGAGCACAAAATTTGGGAGATATATTATTATGTCTTGTAAAAGAAGAACAGAAAAACCGGc                                 >  2:1602668/1‑150 (MQ=25)
                                                                                                                        tccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACAACCACAAGGGGGGGGGGCCCCATCGAGCACAAAATTTGGGAGATATATTATTAAGTATTGTTACACAAGAAAAGAAAAACGGGcgt                               >  2:3835433/1‑148 (MQ=16)
                                                                                                                        tccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCAAAACCCCACCCCCCAGGGCGGGGGGCCCCCGCGAGCCCAAATTTTGGGTGATAAGTTCTTATGTGGTGTAACAGAAGAACACGaaaaagggcgt                                >  2:3268650/1‑143 (MQ=16)
                                                                                                                        tccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGCCCCCCACCCGAAACCACACCCACCACGGGGGGGGGGGCGAGCGAGCACAAATTTTGTGTGATGTATTATTATGTCTGGTAACAGAAGAACAGAAAAAAGGGcgt                               >  2:3364224/1‑148 (MQ=16)
                                                                                                                            tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACTCCCACAAGGGAGAGGGGGCCCATCGAGCACAAATTTTTAGAGATAAATTCTTATGTGTTTTAACAGAATAACAGAAAATCGGGCtttttc                           >  2:3988397/1‑148 (MQ=255)
                                                                                                                              tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCGCCCGCACCCGAAACCACGCCCACAAGGGGGGGGGGACCCATAGAGCCCAAAATTTGTGGTGTATGATAATGTGTGGTGTGACACAAGAAAAGAAAAACGGGCGTTTgggc                         >  2:474988/1‑150 (MQ=255)
                                                                                                                              tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                         <  2:6626478/150‑1 (MQ=255)
                                                                                                                               gggtgtgggCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                        <  1:1249053/144‑1 (MQ=255)
                                                                                                                                ggtgtgggCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:4138803/145‑1 (MQ=255)
                                                                                                                                ggAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  2:4213785/150‑1 (MQ=255)
                                                                                                                                    ggggCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                   <  1:4242859/149‑1 (MQ=255)
                                                                                                                                        cgggggggCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg               <  2:5712369/148‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  1:2593539/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  1:6570135/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:2022569/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGGCATCCGCGCGCACGGTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  2:6074169/150‑1 (MQ=255)
                                                                                                                                          ggggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGCTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg             <  1:4699597/150‑1 (MQ=255)
                                                                                                                                            ggTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  2:5753495/150‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  2:5753335/150‑1 (MQ=255)
                                                                                                                                             gTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag          <  2:6270103/150‑1 (MQ=255)
                                                                                                                                                  ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  1:2033530/150‑1 (MQ=255)
                                                                                                                                                  ggCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  2:1999485/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  2:726673/150‑1 (MQ=255)
                                                                                                                                                     aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  2:434100/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  NC_000913/164432‑164730

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: