Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,613 T→A 11.9% intergenic (+79/‑117) 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,6130TA11.9% 857.1 / 26.6 295intergenic (+79/‑117)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (42/203);  new base A (35/2);  total (90/205)
Fisher's exact test for biased strand distribution p-value = 7.55e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

CCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGC  >  NC_000913/164467‑164762
                                                                                                                                                  |                                                                                                                                                     
 caccctggTCGGGCGGGCCGCCCCGGCAAATACGGCCCCGCCGCGACGGACGAAATATTTTTGGTACGGTGGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                                                   <  2:5325098/143‑1 (MQ=25)
     ctggTCGGGCCGGACCCCCCGGAAATCCCGGCCCCGACGCGCGGGCCGAAAAAGTTTTGGTTGGGTGGGTTTTCTCCCCCCCCCGGTGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                                                               <  2:2028699/147‑1 (MQ=255)
      aTGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                              <  2:8145478/150‑1 (MQ=255)
       tGTCGGGCCGGACGCCCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                             <  2:8215960/150‑1 (MQ=255)
        gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCTCACCCTAAACCTCTCCCTCaatg                                                                                                                                            >  2:8321813/1‑148 (MQ=255)
        gTCTGGCCGGACGACCCGGAAAATACTGCACCGCCGCGACGGACGAAAATTTATTCTGTCGGTTGGTTTTCTCCCTCTCCGTGTGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                            <  2:4380714/150‑1 (MQ=255)
         tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTTCGCCCACACCCACACCCACACCCGCAAggg                                                                                                                                           >  2:7901970/1‑150 (MQ=37)
              ccGGACGCCCCGCCAAATACTGCACCGCCGCGACGGACTAAAAAGTATTCGTATGGTTGTTTTCCTCCCCCTCCCGGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                      <  2:2382922/150‑1 (MQ=255)
               cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGGCCCGCACCCAAACCCACACCCGCAAGgggggggg                                                                                                                                     >  1:4448471/1‑150 (MQ=25)
                   cgacCCGGCAAATACTGCACCGACGCTCCGGACTACAAATTATTCGTATGGTTGGTTTTCTCCCTATCCGGGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                                                 <  2:5692465/150‑1 (MQ=255)
                   cgacCCGGCAAATACGGCCCCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCGCCCCCTCCCGGTGGGTGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                                                 <  2:10717633/150‑1 (MQ=255)
                     acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCCCACCCGCACCCACACGGGGGCGGGGGcccg                                                                                                                               >  1:8600067/1‑148 (MQ=25)
                       ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCGAACCCACACCCACAAGGGGGCGGGGGCCGCTc                                                                                                                             >  1:10193711/1‑150 (MQ=37)
                       ccGGCAAAAACTGCACCGACGCGAGGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                                             <  2:6451638/150‑1 (MQ=255)
                        cGGCAAATCCGGCCCCGTCGCGAGGTACTAAAATTTTTTCTTATGGTTGGTTTTCCCCCTCTCCGTGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                                                            <  1:9294079/150‑1 (MQ=255)
                        cGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCTCTCCCCCAACGGGGGGGGGGCCgcgcg                                                                                                                            >  2:7364037/1‑146 (MQ=37)
                         ggCAACTACTGCACCGCCGCGACGGACTAAAATGTAGTCTTACGGTTGGTTTTCTCCCCCGCCCGGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                                                           <  2:9786042/150‑1 (MQ=255)
                          gCAAATACTGCACCGACGCGACGGACGAAAAAGTGTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCACAACCCCACCCCCACCGGGGAGGGGGCCgggaggg                                                                                                                          >  1:2340931/1‑144 (MQ=21)
                                aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAACCCACACCCCCACGGGGGAGGGGCCCGCTCGTGCACaaa                                                                                                                    >  1:7103696/1‑150 (MQ=25)
                                 cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCGAACCCTCGACCCCCAGGGCGAGGGGGCCGATCGAGCACACAt                                                                                                                   >  2:7283855/1‑150 (MQ=37)
                                   gCATCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGAGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                                                 <  2:8061770/150‑1 (MQ=255)
                                    cACCGACGCGACTTACTAAAATTTTTTCTTATGTTTGTTTTTCGCCCCCTCCGTGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                <  1:1123953/150‑1 (MQ=255)
                                    cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCGCACCCGAACCCACTCCCACAAGGGCGAGGGGGCCGAGCGAGcacacaattg                                                                                                                >  2:7299431/1‑144 (MQ=37)
                                      ccgccgCGAGGGTCGAAAATTTTTTGGTTGGGTTGGGTTCCCCCCTCGCCGGGTGGGTGGGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                              <  1:11040002/150‑1 (MQ=255)
                                      ccgccgCGACGGACTAAAATGTATGCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                              <  1:10866627/150‑1 (MQ=255)
                                      ccGACGCGCCGGACTAAAAAGTTTTCGTCCGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                              <  1:9889806/150‑1 (MQ=255)
                                      ccGACGCGACGGCCGAAAATTTTGTCGTAAGGTTGGTTTTCCCCCTCTCCGTGTGGGTGTGGGCGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                              <  1:10735636/150‑1 (MQ=255)
                                      ccGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                              <  2:11093008/150‑1 (MQ=255)
                                      ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCAAACCCAAAACCACACACACAAGGGGGAGGGGGCCGAAAGAGCCCAAAatttta                                                                                                              >  1:9098699/1‑148 (MQ=21)
                                           gcgacgTACTAAAATGTTTTCTTTTGGTTGTTTGTCTCCCTCTCCGTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                                                         <  1:4287811/150‑1 (MQ=255)
                                           gcgacgGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                                                         <  1:9973580/150‑1 (MQ=255)
                                             gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCACACCCCAACCCCCACCCGCACGGGGGGGGGGCCCGGGCCCGCACAAAATTTGTGCGatat                                                                                                       >  2:10064247/1‑150 (MQ=21)
                                               gggactAAAATTTATTCTTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                                                     <  2:8706149/149‑1 (MQ=255)
                                                ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCCCACCCAAAACCGCAACCACCAGGGAGGGGGGGCCCGGCGAGCACAAAATTTTAGTGGTAtatt                                                                                                    >  1:8162743/1‑150 (MQ=16)
                                                    gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAACCCACACCCACACCGGAGAGGGGACCCAGCGACCACCAAAATTTAGAGATAtcttattc                                                                                                >  2:4516856/1‑150 (MQ=21)
                                                     aaaaaGTATTCGTACGGTTGGTTTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                                                               <  2:3725573/150‑1 (MQ=255)
                                                       aatttATTCGTTTGGTTGTTTTTCGCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                                                             <  1:1637537/146‑1 (MQ=255)
                                                          gTATGCGTAAGGTTGGTTTACTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                          <  2:4318428/150‑1 (MQ=255)
                                                               cGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                                                     <  2:7111709/150‑1 (MQ=255)
                                                                gTACGGTTGGTTTTCTCCCTCGCCCTGTGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                                    <  2:21645/150‑1 (MQ=255)
                                                                  atggttggttTTCCCCCGCTCCGTGGGGGTGGGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                                                  <  1:4068517/148‑1 (MQ=255)
                                                                       tgttttaCTCCCTCGCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                                             <  1:6248069/147‑1 (MQ=39)
                                                                       tggttTTCTCCCTCTCCCTGGGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                                             <  1:2811416/150‑1 (MQ=255)
                                                                       gggttGTCTCCCTCCCCCTGTGGGTGCGGGCCGGGGTGAGCGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                                             <  1:8967225/149‑1 (MQ=255)
                                                                         gtgttCGCCCTCCCCCGGGGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                                                           <  2:3451114/147‑1 (MQ=255)
                                                                             tctccctctcccTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                                                       <  2:1355471/150‑1 (MQ=255)
                                                                                ccctctccctGTGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAACTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                                                    <  2:6974794/150‑1 (MQ=255)
                                                                                  ctctccctGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                                                  <  1:1817178/150‑1 (MQ=255)
                                                                                   tctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCTCCCGCACCCACACCCTCTCCCACAAGGGGGTGGGGCCCGGTCGAGCACAAATTTTTAGAGATATAATATTAGGTCTTGTAACAGAAGGACAGTAAATCGGGc                                                                 >  2:2680202/1‑150 (MQ=21)
                                                                                      ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCCAAACCGCACCCGCCAGGGAGAGGGGACCCCACCAGCACAAAATTTGTGAGATAAATTATTATGTGGTGTAAAAGAAGTACAGAAAAAAGGGCttt                                                              >  2:10982265/1‑150 (MQ=16)
                                                                                           ggggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGCTCGCGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCGTTTgcgc                                                         <  2:9542644/149‑1 (MQ=255)
                                                                                             ggtgtgggCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  2:3217462/145‑1 (MQ=255)
                                                                                             ggtgcgggCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  2:4871928/145‑1 (MQ=255)
                                                                                             ggAGTGGGCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                                                       <  1:7649937/150‑1 (MQ=255)
                                                                                              gagtGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTg                                                      <  2:284490/150‑1 (MQ=255)
                                                                                               cgtgggCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                                                     <  1:9563551/147‑1 (MQ=255)
                                                                                               agagGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCACACCCAAACCCACTCCCACAAGGGAGAGGGGGCCCCTCGACCCCAAAATTTGTGATATATATTCTTAATTATTTTTACAGAAAAAAAAAAAAAAGGGCTTTTTCGCCAGa                                                     >  2:2597859/1‑150 (MQ=255)
                                                                                                   ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                                                 <  2:8885283/150‑1 (MQ=255)
                                                                                                      cGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                              <  2:7633023/150‑1 (MQ=255)
                                                                                                      cGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                              <  1:7535004/150‑1 (MQ=255)
                                                                                                      cGGGGTGAGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                              <  2:8997430/150‑1 (MQ=255)
                                                                                                      cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                              <  1:10113578/150‑1 (MQ=255)
                                                                                                       ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  1:5923874/150‑1 (MQ=255)
                                                                                                       ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  1:6684250/150‑1 (MQ=255)
                                                                                                       ggggTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  1:4676592/150‑1 (MQ=255)
                                                                                                           tGAGGGCCTCAGCGCGCACGTTCACCCACACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga                                         <  1:5645625/150‑1 (MQ=255)
                                                                                                             cgggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACACATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa                                       <  2:9398732/149‑1 (MQ=255)
                                                                                                                       cgcgcACGTTCCCCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGcc                             <  1:10554525/150‑1 (MQ=255)
                                                                                                                         cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                           <  1:3539675/150‑1 (MQ=255)
                                                                                                                           cACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                         <  1:588296/150‑1 (MQ=255)
                                                                                                                            aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGaa                        <  1:4112731/150‑1 (MQ=255)
                                                                                                                             cGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt                       <  1:1363014/150‑1 (MQ=255)
                                                                                                                              gTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATg                      <  2:8856565/150‑1 (MQ=255)
                                                                                                                               ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGa                     <  1:10006682/150‑1 (MQ=255)
                                                                                                                                  accctcaccctAACCCTATACCTCAAGGGAGAGGGGACCGATAGAGCCCAAAATTTTATAGATATCATATTAGGTAGTGTTACAGAAGAACAGAAAATCGGGCTTTTGCGCCGGAAAATTGCGGGGAAAAAGCCTGGCCGGGAAAGGCcg                  >  2:2729530/1‑150 (MQ=255)
                                                                                                                                    cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcgcg                >  2:10382215/1‑150 (MQ=255)
                                                                                                                                    cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTGGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACcgcg                >  2:7429694/1‑150 (MQ=255)
                                                                                                                                     ctcaccctAACCCTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               <  1:3443564/150‑1 (MQ=255)
                                                                                                                                     ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               >  2:6773450/1‑150 (MQ=255)
                                                                                                                                     ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               <  1:10443954/150‑1 (MQ=255)
                                                                                                                                     ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               <  2:9194379/150‑1 (MQ=255)
                                                                                                                                     cccaccctACCCCTCTCCCTCACGGGAGAGGGGACCGAGCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAACTCGGGCTTTTGCGGCTGAATATTGCGGAGAATAAGCATGGCCGGGAATGACCGCGa               <  1:7299431/148‑1 (MQ=255)
                                                                                                                                       caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc             <  1:6894168/150‑1 (MQ=255)
                                                                                                                                           ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAt         <  1:5376180/150‑1 (MQ=255)
                                                                                                                                            tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt        <  1:5032503/150‑1 (MQ=255)
                                                                                                                                            tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtc        <  1:7904878/150‑2 (MQ=255)
                                                                                                                                            gtccccTCTCCCTTGAGGGAGAGGGTACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt        >  2:4426366/4‑150 (MQ=255)
                                                                                                                                             aaCCCTCTCCCTCTAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGGGAGCCAAttt       <  2:6623919/150‑2 (MQ=255)
                                                                                                                                               cccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCAATTGGa     >  2:10120200/1‑150 (MQ=255)
                                                                                                                                               cccTCTCCCTCAAGGGAGAGGGGACCGCTCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     <  2:4781742/150‑1 (MQ=255)
                                                                                                                                                ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCAATTGGAc    >  2:10999697/1‑150 (MQ=255)
                                                                                                                                                ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGCCCGCGAGCCCATTGGAc    >  2:6475710/1‑150 (MQ=255)
                                                                                                                                                ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    >  2:10695302/1‑150 (MQ=255)
                                                                                                                                                ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTGTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    <  1:7414529/150‑1 (MQ=255)
                                                                                                                                                  tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAACATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAGTGGACGc  >  1:6585537/1‑150 (MQ=255)
                                                                                                                                                  tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc  >  1:3187002/1‑150 (MQ=255)
                                                                                                                                                  tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc  >  1:6839672/1‑150 (MQ=255)
                                                                                                                                                  tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc  >  2:5955855/1‑150 (MQ=255)
                                                                                                                                                  tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc  >  1:6944160/1‑150 (MQ=255)
                                                                                                                                                  tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc  >  1:9731620/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGC  >  NC_000913/164467‑164762

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: