Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,613 T→A 9.8% 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,6130TA 9.8% 714.9 / 18.7 239intergenic (+79/‑117)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (35/170);  new base A (22/2);  total (66/173)
Fisher's exact test for biased strand distribution p-value = 2.61e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

CCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGC  >  NC_000913/164466‑164762
                                                                                                                                                    |                                                                                                                                                     
 ccAAACTTGTCGGGCCGGACGACCCGTCAAATCCTGCACCGACGCGACGGACGAAAATGTTTTGGTATGGTTGGTTTTCTCCCTCTCCCGGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                                                    <  2:6090762/150‑1 (MQ=37)
   aaccctGTCGGGCCGGACGCCCCGGCAAATACTGCACCGACGCGACGGACGAAATATTATCGGAATGGTGGGTTTGCTCCCTCCCCCGGGGGTTGGGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                                                  <  2:5041571/145‑1 (MQ=255)
   aaaCATGTCTGGCCGGACGACCCGCCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGACGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTccc                                                                                                                                                  <  2:4585300/150‑1 (MQ=255)
       aTGTCTGGCCGGACCACCCGGAAAATACGGCCCCGACGCGACGGACTAAAATTTTTTGGTATGGTTGTTTTTCTCCCACCCCGTGGGGGATTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                              <  2:3175311/150‑1 (MQ=255)
        tggcTGGCCGGACGACCCGGCAAAAACTGCACCGACGCGACGGACGAAAAAGTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                                                             <  2:2008261/147‑1 (MQ=255)
         gTCTGGCCGGACGACCCGGCAAAAACTGCACCGACGCGACGTACCAAAAATTATTCTTAAGGTTGGTTTTCCCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                                                            <  2:1789872/150‑1 (MQ=38)
             ggCCGCCCGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATGCGTACGGTTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGgaga                                                                                                                                        <  1:6103928/150‑1 (MQ=38)
              gCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCAAACCCTCAAGGGagag                                                                                                                                       >  2:7287283/1‑150 (MQ=255)
              gCCGGACGACCCGCCAAATCCGGCACCGACGCGACGGACGAAAAAGTATGCGTATGGTGGGTTTTCTCCCTCACCCTGTGGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                       <  1:8389784/150‑1 (MQ=255)
               ccGGACGCCCCGGCAAATACTGCACCGACGCGCCGGACGAAAACGTATTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                      <  2:5463411/150‑1 (MQ=38)
                cGGACGCCCGGTAAAATACTGCACCGCCGCGACGTACTAAAATTTTTTCGTATGGTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAggg                                                                                                                                     <  1:2853805/150‑1 (MQ=255)
                     gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCGCTCCCACCAGGGAGAGGGGCCCg                                                                                                                                >  2:4881150/1‑150 (MQ=38)
                         ccgCAAATTACTGCCCCGAGGCTCCGGACTAAATTTTTTTGTTATGGTTGGTTTTCCCCCTCCCCGTGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                                                             <  1:6570261/148‑1 (MQ=255)
                            cAAATCCTGCCCCGACGCGCCGGCCTAAAAAGTATTCGTAAGGTTGGTTTTCGCCCTCCCCCTGGGGGTGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                                                         <  2:6810620/150‑1 (MQ=255)
                                    tGCACCGACGCTACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                                                  <  2:4762457/150‑1 (MQ=255)
                                    tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCCAACCCACACCCACAAGGGCGAGGGGGCCGCAGGAGCCCaaaaa                                                                                                                  >  1:4913708/1‑148 (MQ=25)
                                    ggccccgccgCGACGCACTAAAATTTATTCTTTTGGTTGTTTTTCCCCCCCCCCGGGTGGGCGCGGGCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                                                  <  1:797503/146‑1 (MQ=255)
                                      cACCGCCGCGCCGGACTAAAATTTTTTCTTTTGTTTGGTTTTCGCCCGGGCGGGGGGGGTGGGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                <  2:2862715/150‑1 (MQ=255)
                                          gacgcgacgGACTAAAAAGTTTTGGTTAGGTTGGTTTTCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                                                            <  1:6016581/150‑1 (MQ=38)
                                                  ggacgAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                                                    <  1:7698925/150‑1 (MQ=255)
                                                  ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGCCCCCACCCAAAACCACACCCCCAACGGGGGGGGGGCCGAACGCGCCCAAAATTTTTGAGATAtatt                                                                                                    >  1:974418/1‑150 (MQ=18)
                                                  ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCGCCCCAAACCACAACCACAAAGGAGAGGGGGCCGCGAGAGCACCAAATTTgtgaggtttaat                                                                                                    >  1:3780794/1‑143 (MQ=16)
                                                  ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTCCACCCGCACCCAAAAACACAACCAAAAAGGGGGGGGGGCCCATCGAGCCCAAAATTTTATAGATAtatt                                                                                                    >  1:4267577/1‑150 (MQ=16)
                                                        aaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCGCACCCACAACGGCGAGGGGGCCCCTACAGCACAAATTATGAGAGAAATCTTATTATg                                                                                              >  2:4527098/1‑150 (MQ=21)
                                                            gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                          <  2:6982084/150‑1 (MQ=255)
                                                             tagccGTAAGGTTGGTTTTCGCCCGCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                                         <  2:4490259/146‑1 (MQ=255)
                                                             tATTGGTATGGTGGGTTTACGCCCTCTCCCTGTGGTTGCGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                                                         <  2:1761942/150‑1 (MQ=255)
                                                                  gTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                                    <  2:7210112/150‑1 (MQ=255)
                                                                  gTATGGTTGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                                    <  2:7420385/150‑1 (MQ=255)
                                                                  gTAAGGTTGGTTTTCTCCCTCTCCCCGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                                                    <  2:4934525/150‑1 (MQ=255)
                                                                    atggttggttTTCCCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                                                  <  2:615431/148‑1 (MQ=255)
                                                                     agtttggttTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                                                 <  1:4441059/147‑1 (MQ=255)
                                                                        ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAACCGCACCCAAAACCACACCCCCCAGGGGGTGGGGCCCGAGCGAGCCCAAAAGTTGTGAGATAAGTTATTGTGTGTTGTTAGAgaagaa                                                                              >  2:850550/1‑150 (MQ=12)
                                                                        ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACTACCACAAGCGAGGGGGGGCCCAGCGCGCGCAAATTTTGTGAGGTATATTATTTTGTGGTGTTAAagaacta                                                                              >  1:6907370/1‑147 (MQ=12)
                                                                         tggttTCCTCCCCCTCCCTGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                                             <  1:5874214/150‑1 (MQ=255)
                                                                          gtttttCTCCCTAACCCTGTGGGAGTGGCCGGGGGGGCGGGCCTCAGCGCGCACGTTCACACTCACCCTAACCCACTCCCTCAAGGGAGAGGGGACCGATAGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                            <  1:6485182/148‑1 (MQ=255)
                                                                            ttttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCCAACCCACTCCCACAAGGGGGGGGGGACCGCAAGAGCCCAAATTTTGTGGGGTATCTTCATCTGTCTTGTTACAGAAGAACAGa                                                                          >  2:2118352/1‑150 (MQ=21)
                                                                                 tccctctccctGTGGGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                                                     <  1:6336314/150‑1 (MQ=255)
                                                                                  ccctctccctGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                                                    <  1:8343241/150‑1 (MQ=255)
                                                                                  ccctctcccgGGGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                                                    <  1:729024/150‑1 (MQ=255)
                                                                                        ccgtgtgggtgtgggCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                                                              <  1:6870721/147‑1 (MQ=255)
                                                                                        ccctGTGGGCGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                                                              <  1:3516318/150‑1 (MQ=255)
                                                                                          cggggggAGCGGGCCGGGGTGCGGCCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                                                            <  1:2315442/146‑1 (MQ=255)
                                                                                             tgGGAGAGGGCCGGGGTGCGCGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                         <  2:5523194/150‑1 (MQ=255)
                                                                                             ggggAGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                         <  2:8194598/149‑1 (MQ=255)
                                                                                             ggggAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                                                         <  1:2630080/149‑1 (MQ=255)
                                                                                                     ggCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCCCCCCCCAACCCACACCCACAACGGACAGGGGCCCGAACCAGCCCCAAATTTGAGATATAAAATAATATGTCTTGTAAAACAACAAAACAAAAACCGGCGTTTTCGccagtaaat                                                 >  2:1282281/1‑143 (MQ=255)
                                                                                                     ggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCCCAACCGCCAGGGGGCGGGGGACCCTCGAGCCCAAAAAGTGTGTGATATATTAGTGTGTCGTGGTAAAGAAGGAAAGCCAAAGGGGCTTTTTCGccggtagat                                                 >  2:5304429/1‑143 (MQ=255)
                                                                                                       ccGGGGTGCGGGCATCAGCGCGCACAGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg                                               <  1:1853170/150‑1 (MQ=255)
                                                                                                        cGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc                                              <  1:3274301/150‑1 (MQ=255)
                                                                                                         ggggTGAGGGCCTCAGAGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  1:629969/150‑1 (MQ=255)
                                                                                                         ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  1:4734530/150‑1 (MQ=255)
                                                                                                         ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  2:4663731/150‑1 (MQ=255)
                                                                                                         ggggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg                                             <  1:6123261/150‑1 (MQ=255)
                                                                                                           ggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga                                           <  2:5246802/150‑1 (MQ=255)
                                                                                                                 ggCCTCCGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                     <  1:2433577/150‑1 (MQ=255)
                                                                                                                   cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc                                   <  2:419167/150‑1 (MQ=255)
                                                                                                                    ctCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa                                  <  2:6534914/149‑1 (MQ=255)
                                                                                                                      cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATg                                <  2:7532947/150‑1 (MQ=255)
                                                                                                                          gcgcCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg                            <  1:7135257/150‑1 (MQ=255)
                                                                                                                          gcgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCg                            <  2:6436580/150‑1 (MQ=255)
                                                                                                                           cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                           <  2:8246550/150‑1 (MQ=255)
                                                                                                                           cgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCgg                           <  2:4388591/150‑1 (MQ=255)
                                                                                                                             catgTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                         <  1:1892250/147‑1 (MQ=255)
                                                                                                                             cACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                         <  2:5724611/150‑1 (MQ=255)
                                                                                                                              gctgatgccctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGaa                        >  2:1862559/8‑150 (MQ=255)
                                                                                                                               cGTTCACCCTCACCCTAACCCTCTCACTCAAGTGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTTCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt                       <  1:557666/150‑1 (MQ=255)
                                                                                                                               cGTTCACCCTCACCCTAACCCTCCCCCTCAAGGGAGAGGCGACCGATCGAGCACAACATTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAACAACGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAAt                       <  2:7761854/150‑1 (MQ=255)
                                                                                                                                    accctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                  <  2:4982242/150‑1 (MQ=255)
                                                                                                                                      cctcaccctAACCCTCTCCCTCAAGGGAGGGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCAGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcgcg                >  1:5097567/1‑150 (MQ=255)
                                                                                                                                      cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGGATATTGCGGAGAAAAAGCATGGCCGGGTATGACcgcg                >  1:1362472/1‑150 (MQ=255)
                                                                                                                                       ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACCGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               <  1:4082007/150‑1 (MQ=255)
                                                                                                                                       ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcgcgg               >  1:7881377/1‑149 (MQ=255)
                                                                                                                                       ctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa               <  2:3679711/150‑1 (MQ=255)
                                                                                                                                        tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg              <  2:7362107/150‑1 (MQ=255)
                                                                                                                                        tcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAg              >  1:2543115/1‑150 (MQ=255)
                                                                                                                                         caccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGc             <  2:3598091/150‑1 (MQ=255)
                                                                                                                                            cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCaa          >  2:5978989/1‑150 (MQ=255)
                                                                                                                                            cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa          <  1:245186/150‑1 (MQ=255)
                                                                                                                                            cctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa          >  1:7031386/1‑150 (MQ=255)
                                                                                                                                              tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt        >  2:677715/1‑150 (MQ=255)
                                                                                                                                              tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt        >  1:6694310/1‑150 (MQ=255)
                                                                                                                                              tAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt        >  2:2437010/1‑150 (MQ=255)
                                                                                                                                              gaaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAAtt        <  1:1734868/149‑1 (MQ=255)
                                                                                                                                               tccccTCTCCCTTGAGGGAGAGGGTACCGATCGAGCACAAAGTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       >  2:3993192/3‑150 (MQ=255)
                                                                                                                                               tccccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       >  2:7944323/3‑150 (MQ=255)
                                                                                                                                               aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       <  2:4267578/150‑1 (MQ=255)
                                                                                                                                               aaCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTg       <  2:3780794/150‑1 (MQ=255)
                                                                                                                                                aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGAGATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgg      <  1:3994341/150‑1 (MQ=255)
                                                                                                                                                 cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     <  1:2024353/150‑1 (MQ=255)
                                                                                                                                                 cccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa     <  1:1302563/150‑1 (MQ=255)
                                                                                                                                                  ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    <  2:2979721/150‑1 (MQ=255)
                                                                                                                                                  ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    >  2:7981085/1‑150 (MQ=255)
                                                                                                                                                  ccTCTCCCGTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGAc    >  1:3139035/1‑150 (MQ=255)
                                                                                                                                                   ctctCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGACTGACCGCGAGCCAATTGGACg   >  2:1208957/1‑150 (MQ=255)
                                                                                                                                                   ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg   <  1:3721819/150‑1 (MQ=255)
                                                                                                                                                   ctctCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAGTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACg   <  1:6719833/150‑1 (MQ=255)
                                                                                                                                                    tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc  >  2:2637432/1‑150 (MQ=255)
                                                                                                                                                    tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc  >  1:8174006/1‑150 (MQ=255)
                                                                                                                                                    tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc  >  1:4621464/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGC  >  NC_000913/164466‑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: