Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,581 A→C 19.5% 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,5810AC19.5% 208.5 / 25.7 128intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (37/64);  new base C (4/23);  total (41/87)
Fisher's exact test for biased strand distribution p-value = 3.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729
                                                                                                                                                     |                                                                                                                                                    
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAATAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                      >  2:3208393/1‑150 (MQ=255)
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                      >  2:2148546/1‑150 (MQ=255)
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                      >  2:5240606/1‑150 (MQ=255)
ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGAACCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGCTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCAGGGGTGCGGGCa                                                                                                                                                      >  2:5762826/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     >  1:5795045/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     >  2:4069562/1‑150 (MQ=255)
 cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     >  2:4069496/1‑150 (MQ=255)
   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   >  2:1939067/1‑150 (MQ=255)
       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAgcgc                                                                                                                                               >  1:4583478/1‑150 (MQ=255)
           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                           >  1:5566750/1‑150 (MQ=255)
               agagATGAAAGGGCGATATACCAAACATGTCTGGCCGGACGACCCGTCAAAGACTGCACCGACGCGACGGACGTAAAAGTATACGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGCGAGGGCCGGTGTGAGTGCAGCAGCGCGCACGGtc                                                                                                                                       >  1:2404596/1‑150 (MQ=37)
                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccc                                                                                                                                   >  2:747474/1‑150 (MQ=255)
                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTggaccca                                                                                                                                  >  1:3819868/1‑149 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctc                                                                                                                                 >  1:4191711/1‑150 (MQ=255)
                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccgc                                                                                                                                 >  2:4845379/1‑148 (MQ=255)
                      aacgggCGTTATCCCAAACAGGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAATGTATTCTTATGGTTGTTTTTCGCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                <  2:3312883/147‑1 (MQ=38)
                      aaaGGGCTTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                <  2:3312890/150‑1 (MQ=255)
                        aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacccacccc                                                                                                                              >  1:4362331/1‑145 (MQ=255)
                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaccccaaacca                                                                                                                            >  2:2782739/1‑143 (MQ=255)
                          ggCGTTAGCCCCACCCGGCCTGGCCGGCCGCCCCGGCAAATCCTGCCCCGACGCGACGGACTAAAAAGTATTCGTAGGGTGGGTTTTCGCCCTCTCCCGGTGGGCGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTcaccctcaccct                                                                                                                            <  2:4656628/150‑1 (MQ=255)
                              ttATCCCAACCTTGTGGGGCCGGACGACCGGGAAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTACGGTTGGTTTTCTCCCTGGCCGTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                        <  2:4218656/150‑1 (MQ=255)
                              ttATCCCAACCATGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                        <  2:3689041/150‑1 (MQ=255)
                                 tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACACCCCAAccccc                                                                                                                     >  1:1209411/1‑148 (MQ=38)
                                  cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCACCCGCACCCCAAcccgca                                                                                                                    >  2:3862555/1‑147 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGTCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCCCACCCCAAaccacaa                                                                                                                   >  1:4089375/1‑146 (MQ=38)
                                   ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCACCCCAAaccacac                                                                                                                   >  1:5790161/1‑146 (MQ=38)
                                    caacccggTCGGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCGGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                  <  2:5289973/143‑1 (MQ=255)
                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCACAcccaccac                                                                                                                  >  2:3940269/1‑145 (MQ=38)
                                    cAAACATGGCTGGCCGGAGGACCGGGAAAATACGGCACCGAGGCGAGGGAGTAAAATGTTTGCGTTTGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTcc                                                                                                                  <  2:4191711/150‑1 (MQ=25)
                                     aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCCCACCCAAAccccacacc                                                                                                                 >  1:4131386/1‑144 (MQ=37)
                                        cATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCGCAACCAAAACCACTCCcaca                                                                                                              >  1:1939662/1‑150 (MQ=37)
                                              tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCGAACCCACTCCCACCAgggtg                                                                                                        >  2:2420776/1‑148 (MQ=38)
                                                gCCGGACGACCCGCAAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTAAGGTTGGTTTTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  2:396764/150‑1 (MQ=255)
                                                gCCGGACCCCCCGGCAAATACTGCACCGACGCGAGGCACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCCCCGTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                      <  2:5085468/150‑1 (MQ=38)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCACCCCCCAAACCCCGCCCCCAAGGGCGAggg                                                                                                    >  1:862568/1‑150 (MQ=37)
                                                  cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                    <  2:3529648/150‑1 (MQ=255)
                                                    gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCTCGCACGTTGGCCCGCACCCAAAACCACGCCCACCAGGGAggggggg                                                                                                  >  1:146409/1‑149 (MQ=37)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGGGCGCCCACCCCCAAACCCACACCCACCACGGCGGGGGGGCCgcgc                                                                                            >  1:4198110/1‑147 (MQ=21)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   <  2:3573016/95‑1 (MQ=255)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                   >  1:3573016/1‑95 (MQ=255)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCAAAACCACACCCACAAGGGGGGGGGGGcccaag                                                                                            >  1:887148/1‑146 (MQ=25)
                                                          ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCACACCCAAAACCAAACCCACAAggggggggggggcgggc                                                                                            >  1:518888/1‑143 (MQ=21)
                                                            gtaaaaTACTGCACCTACGCGTCGGACTAAAATCTATTCTTTAGTTTTTTTTTCTCCATATCCCTTTGGTTGTGGGCGGGGGTTTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:2404596/147‑1 (MQ=255)
                                                            ggaaaaTCCTGCCCCGCCGCGACGGCCGAAAATGTAGTGGTACGTTTGGTTTGCCCCCCCCCCGGGGGGTTGTGGGCCGGGGTGGGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:2072678/147‑1 (MQ=255)
                                                              aaaaaACTGCACCGCCGCGACGCACTAAAACCTAGTCTAATGGTTGGTTTCCACCCTCCCCGGGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                        <  1:2875702/149‑1 (MQ=255)
                                                                      gCACCGACGCGACGGACGAAAATGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                <  2:3433062/150‑1 (MQ=255)
                                                                       cACCGCCGGGCCGTACGAAATTGTTTTCTTATGTTTGGTTTCCCCCCTCTCCGGGGGGGTGCGGGCGGGGTGGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  1:4460303/150‑1 (MQ=255)
                                                                       cACCGACGCGACGGACGAAAAATTATTCTTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  2:1916106/150‑1 (MQ=255)
                                                                       cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                               <  2:1916362/150‑1 (MQ=255)
                                                                        aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCGCCCACACCCAAAACCACACCCACAAGGGCGCGGGGGCCGCACGAGCACAAAATTTg                                                                              >  1:2379124/1‑150 (MQ=25)
                                                                        aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACACACACCCAAAAACAAAACCAAAAGGGAGAGTGGGCCCAACGAACACAAAatttt                                                                              >  2:553910/1‑149 (MQ=18)
                                                                         ccgccgCGCCGGACGAAAAAGTATTCGTATGGTTGGTTTTCGCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                             <  2:2963718/150‑1 (MQ=255)
                                                                         ccGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCTAAACCACACCCACAAGGGGGAGGGGCCCGATAGCGCACaaaaattgt                                                                             >  2:2862913/1‑144 (MQ=25)
                                                                               cgcggtacgAAAACTTATGCTTAGGGTTGTTGTCCTCCCACGCCGGGTGGGTGTGGGCGGGGGTGCGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTCGAGAGata                                                                       <  2:3914976/144‑1 (MQ=255)
                                                                                 acgtactAAAAAGTATTCTTACGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                     <  1:1334586/150‑1 (MQ=38)
                                                                                 acggactAAAATGTATTTGTATGGTTGGTTTTCCCCCTCCCCCTGTGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                     <  1:114993/150‑1 (MQ=255)
                                                                                          aaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                            <  2:4381737/150‑1 (MQ=255)
                                                                                           aaGTTTTCTTACGGTTGGTTTTCTCCCTCGCCCTGTGGGCGCGGGCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTc                                                           <  1:1810135/150‑1 (MQ=255)
                                                                                             gTATTCGTAAGGTTGGTTTTCGCCCTCGCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                         <  1:1385112/150‑1 (MQ=255)
                                                                                              ttttCGTTGGGTGGGTTTCCCCCCCCGCCCTGTGGGCGAGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                        <  1:2070262/148‑1 (MQ=255)
                                                                                               aTTCCTAAGTTGGGTTTTCGCCCTCTCCGTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                       <  2:5584959/150‑1 (MQ=38)
                                                                                                  cgcacggttgtttTACACCCCCCCCCTGGGGGAGGGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                    <  1:495068/145‑1 (MQ=255)
                                                                                                    tatggttggttTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                  <  1:3508384/147‑1 (MQ=255)
                                                                                                    tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCCCCCCCACCCCAACCCGCTCCCGAAACGGCGAGGGGCCCCAACGCGCACAAAAGGTTGGAGTTATATTCTTCTGTCTTGTAACag                                                  >  2:3907231/1‑150 (MQ=16)
                                                                                                      cggttggttTTCTCCCACTCCCTGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                <  1:4556712/149‑1 (MQ=255)
                                                                                                      aGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                <  1:447783/150‑1 (MQ=255)
                                                                                                       gtttgtttTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                               <  1:1929166/148‑1 (MQ=255)
                                                                                                       ggttggttTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                               <  2:2872899/150‑1 (MQ=255)
                                                                                                                   ccctctccctGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGTGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                   <  1:4845379/150‑1 (MQ=255)
                                                                                                                   ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                   <  2:3218282/150‑1 (MQ=255)
                                                                                                                    cctctccgtGTGGGCGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                  <  1:31702/150‑1 (MQ=255)
                                                                                                                         cccggggggAGTGGGCCGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                             <  2:236705/144‑1 (MQ=255)
                                                                                                                          cctGGGGGTGGGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                            <  1:747474/150‑1 (MQ=255)
                                                                                                                           gtgtgGGAGTGGCCCGGGGTGAGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                           <  2:5498591/149‑1 (MQ=255)
                                                                                                                           ctGTGGTAGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                           <  2:1939662/150‑1 (MQ=255)
                                                                                                                              ggggAGTGGGCCGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                        <  1:4408209/149‑1 (MQ=255)
                                                                                                                               ggtagcgggCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                       <  1:1799261/144‑1 (MQ=255)
                                                                                                                                ggAGAGGGCCGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTGGAGAGATATCCTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                      <  2:1255586/150‑1 (MQ=255)
                                                                                                                                 gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCAACCGCACCCAAACCCGAACCCGCCACGGAGGGGGGACCGATCGAGCACAAATTTTTTGAGATATCATCTTCTGTCGTGTACAAGAAGAACAGCAAAACCGGCTGTTGCGccag                     >  2:5310312/1‑148 (MQ=255)
                                                                                                                                  agagGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                    <  1:1864830/150‑1 (MQ=255)
                                                                                                                                     gggCGGGGGTGCGGCCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                 <  1:1013270/150‑1 (MQ=255)
                                                                                                                                      ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACACCCCCCAGGGGGAGGGGGCCGATAGAGCCCAAAATTTGTGAGATAAATTATTATGTATTTTAACAGAACCACAGAAAAAAGGGCTTTTTGGCCCGAaaat                >  2:3055707/1‑150 (MQ=255)
                                                                                                                                        ccGGGGTGAGGGCATCAGCGCGCACGTTGACCCACAACCAAACCCACAACCACAACGGAGAGGGGACCCATCGAGCACAAAATTTGTTAGTTATATTATTAATTATTGTAAAAGAAGAAAAAAAAAAAGGGCTTTTTAGCCGGAATAATg              >  2:4775643/1‑150 (MQ=255)
                                                                                                                                        ccGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCAAAACCACTCCCACAACGGGGGGGGGGCCCCTAGCGCCCAAAATTTTTGATATATAATATTATTTGTTTTTAAAGAAGAAAACAAAAAACGGCGTTTTAGCCGGAATATTg              >  2:594105/1‑150 (MQ=255)
                                                                                                                                         gggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  1:185420/149‑1 (MQ=255)
                                                                                                                                         cGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  1:1778392/150‑1 (MQ=255)
                                                                                                                                         cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  2:3159654/150‑1 (MQ=255)
                                                                                                                                          ggggggAGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:2646136/150‑1 (MQ=255)
                                                                                                                                          ggggTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  2:4535029/150‑1 (MQ=255)
                                                                                                                                          ggggTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  2:5318152/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATGGTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  2:4534936/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCCTCAGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:1786061/150‑1 (MQ=255)
                                                                                                                                          ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg            <  1:3479189/150‑1 (MQ=255)
                                                                                                                                           gggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTTTAACAGAAGAACAGAAAATCGGTCTTTTGCGCCTGAATATTGCgg           <  2:3641622/150‑1 (MQ=255)
                                                                                                                                           gggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg           <  2:1560555/150‑1 (MQ=255)
                                                                                                                                            ggggAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga          <  2:1609059/147‑1 (MQ=255)
                                                                                                                                            ggTGAGGGCCTCAGCGCGCCCGTTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga          <  2:890129/150‑1 (MQ=255)
                                                                                                                                                 gggCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa     <  2:5270309/150‑1 (MQ=255)
                                                                                                                                                  ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  2:3719479/150‑1 (MQ=255)
                                                                                                                                                  ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa    <  2:5795045/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  NC_000913/164432‑164729

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: