Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,068 A→C 28.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
*CP009273161,0680AC28.5% 107.6 / 28.3 111intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (29/49);  new base C (6/27);  total (35/76)
Fisher's exact test for biased strand distribution p-value = 7.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

GTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAA  >  CP009273/160921‑161212
                                                                                                                                                   |                                                                                                                                                 
gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                 <  2:2926989/150‑1 (MQ=255)
gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                 <  2:1872150/150‑1 (MQ=255)
 tGAGGTTCAAAAAGAGATTAATGGGCTTTATCCCAAAATTGTCTGGCCGTACGACCCGGCAAATACTGCACCGACGCGACGGACTAAAATTTATTCGTTTGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCa                                                                                                                                                <  1:302692/150‑1 (MQ=37)
 tGAGGTGCAAAAAGAGATGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                <  1:4273065/150‑1 (MQ=255)
 tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCa                                                                                                                                                >  1:28374/1‑150 (MQ=255)
  gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTccg                                                                                                                                               >  2:4227836/1‑148 (MQ=255)
  gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTccg                                                                                                                                               >  2:1650116/1‑148 (MQ=255)
   aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgc                                                                                                                                              >  2:67546/1‑150 (MQ=255)
      tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCcgcgcg                                                                                                                                           >  2:1354127/1‑150 (MQ=255)
        cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                         >  2:3548399/1‑150 (MQ=255)
        cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                         >  2:5116542/1‑150 (MQ=255)
         aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                        >  2:4171896/1‑150 (MQ=255)
         aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                        >  2:4736459/1‑150 (MQ=255)
         aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                        >  1:5088696/1‑150 (MQ=255)
            aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGtt                                                                                                                                     >  1:1369515/1‑150 (MQ=255)
             agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgc                                                                                                                                    >  1:3426917/1‑148 (MQ=255)
             agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtc                                                                                                                                    >  1:3627687/1‑150 (MQ=255)
                  ttatcgggCTTTATCCAATATATGTGGGGCCGGACTACCCGAAAAATACTGCACCTACGCTAGGTCCTAATAATTATGTGTATGGTGTTTTTTCTCCCTCGCAGTGTGGTCTTGGGCGGGGTTTGGGGACTCAGCGCGCACGTtcaccct                                                                                                                               <  2:2256333/145‑1 (MQ=255)
                    aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccctcc                                                                                                                             >  2:3525278/1‑149 (MQ=255)
                    aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccaca                                                                                                                             >  1:2620936/1‑150 (MQ=255)
                     aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccaccc                                                                                                                            >  2:4850895/1‑146 (MQ=255)
                     aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcacccaccc                                                                                                                            >  2:4850490/1‑146 (MQ=255)
                        ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTcaccctcaccct                                                                                                                         <  2:3760345/150‑1 (MQ=255)
                        ggCGTTATCCCAAACATGCCTGGCCGGACGCCCCGGCAAATACTGCACCGACGCGAGGGACTAAAAAGTATTCGTAAGGTTGGTTTTCGCCCTCTCCGTGGGGTAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTcaccctcaccct                                                                                                                         <  2:3168533/150‑1 (MQ=255)
                           gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGACAAAGTATTCGTCAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGGCGCCCGCACCCCAAc                                                                                                                      >  2:4955818/1‑150 (MQ=37)
                           gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtgacccccccccccac                                                                                                                      >  2:4955846/1‑146 (MQ=255)
                             tATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCCCCCCAAccc                                                                                                                    >  1:4402033/1‑150 (MQ=38)
                                cccAACCCTGTCGGGCCGGACGACCCGGCAAATACTGCACCGCCGCGACGGACGAAAAAGTATTCGAAAGGTTGGTTTTCGCCCTCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                 <  2:4716410/150‑1 (MQ=38)
                                cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGAACCACACCCAAAaccaca                                                                                                                 >  1:2456368/1‑147 (MQ=37)
                                 ccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGCACCACACCCCAAcccccac                                                                                                                >  2:1778930/1‑146 (MQ=37)
                                      cATGTCTGGCCGGACGACCCGGCAAATACGGCCCCGACGCGCCGGACGAAAAATTATTCGTACGGTTGGTTTTCTCCCGCTCCCGGTGGGTGAGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                           <  2:673732/150‑1 (MQ=255)
                                        tGTCTGGCCGGACGCCCCGGAAAATACTGCACCGACGCGCCGGACGAAAAATTTTTCTTATGGTTGTTTTTCTCCCTCCCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                         <  2:2282094/150‑1 (MQ=37)
                                                cgccccccccGGAAAATCCTGCACCGACGCGAGGGACTAAAATGTATTCTTAGGGTTGTTTTACCCCCTCACCCTGTGGTAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAggg                                                                                                 <  2:1784971/143‑1 (MQ=255)
                                                 tgacgacCCGCCAAATACTGCACCCCCGCGACGGACTAAAAAGTATTCGTATGGTTGTTTTACTCCCACTCCGTGTGGTAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggg                                                                                                <  2:1078/149‑1 (MQ=255)
                                                     gccccGCCAAATCCTGCACCCACGCCACGCCCGAAAACGTTTGCGTACGTGTGGTTTTCACCCTCTCCCTGTGGGCGAGGCCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                            <  2:1010337/148‑1 (MQ=255)
                                                     gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCACTCCCAACAGGGCGGGGGGGCCg                                                                                            >  2:2743526/1‑150 (MQ=25)
                                                     gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGACCCACCCCCAAACCCCCACCCACAAGGGAGAGGGGGccc                                                                                            >  2:2602775/1‑149 (MQ=37)
                                                     gacCCGGCAAATACGGCACCGACGCGACGGACTAAAATTTTTTCTTAAGTTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                            <  1:5146411/150‑1 (MQ=37)
                                                      ccccGGCAAATACCGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTGTTCGCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGa                                                                                           <  2:2923162/149‑1 (MQ=255)
                                                      acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCTGGGTGAGGGCATCAGCGCGCACGTGCAAACACACCCTAACCCACTCCCGCTAGGGAGAGGGGCCCgc                                                                                           >  2:3237423/1‑149 (MQ=37)
                                                        cgggaaaaCACTGCCCCGACGCGCCGGAGTAAAACGTATTCGTCGGGTTGGTTTTCGCCCCCCCCCGGTGGGCGCGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                         <  2:1995798/145‑1 (MQ=255)
                                                        ccGGCCACTACGGCACCGACGCGCCGGACGAAAACGTATTCTTAAGGTTGGTTTTCTCCCACTCCTTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                         <  2:1996034/150‑1 (MQ=38)
                                                        ccGCCAACTACTGCCCCGACGCGCGGGACGAAAAACTATTCGTCAGGTTGGTTGTCGCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCCCGTACACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                         <  2:3194541/150‑1 (MQ=37)
                                                         cGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCGCCCTCTCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                        <  2:1942812/150‑1 (MQ=255)
                                                            cAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCCCCCCAACCCTCACCCGCCAGGGGGGGGGGGCCgcacgcgc                                                                                     >  2:3667724/1‑143 (MQ=25)
                                                               atacgccACCGCCGCGCGGGACCAAAATTTATTCGTACGGTTGGTTTTCGCCCTCCCCCGGGGGTTGGGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcaca                                                                                  <  1:1135749/144‑1 (MQ=255)
                                                                  cTGCACCGACGCGACGCACGAAAATTTTTTCTTATGTTTGTTTTTCGCCCTCACCCTGTGGGCTTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                               <  1:2475672/150‑1 (MQ=255)
                                                                    gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCTCACCCAAACCCACACCCACAAGCGTGAGGGGACCGAGCGAGCCCaaaaat                                                                             >  2:4216038/1‑147 (MQ=25)
                                                                         gaagcgacgCACGATAATGTATTCTTATGTTTGTTTTTCTCCCTCTCCCTGGGGGAGTGGGGGGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACACTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                        <  2:453514/147‑1 (MQ=37)
                                                                          ccgcgccgGACTAAAAAGTTTCCTTCTGGTTGTTTTCCCCCCGCCCCGTGTGGTGGTGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                       <  1:1563818/149‑1 (MQ=255)
                                                                             cgccggacgAAAATTTATTCTTATGTTTGTTTTTCCCCCTCCCCCTGTGGGTGTGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGata                                                                    <  1:350851/147‑1 (MQ=255)
                                                                             cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGata                                                                    <  2:5114760/150‑1 (MQ=255)
                                                                                 ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                <  2:4988914/150‑1 (MQ=255)
                                                                                  gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                               <  2:3854143/150‑1 (MQ=255)
                                                                                   actaaaaaGTTTTCGTATGGTTGGTTTTCTCCCTCGCCCTGTGGGTGTGGGCGGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttct                                                              <  2:4892229/147‑1 (MQ=255)
                                                                                     taaaaaTTATTCTTACGTTGGGTTTTCCCCCTCTCCCGTTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                            <  2:1475967/149‑1 (MQ=255)
                                                                                     taaaaaTTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                            <  2:1475962/149‑1 (MQ=255)
                                                                                     gAAAAAGTATTCGTAAGGTTGTTTTTCTCCCGCTCCCTTTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                            <  2:2059144/150‑1 (MQ=255)
                                                                                        aaaTTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                         <  2:194504/150‑1 (MQ=255)
                                                                                          aGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                       <  2:760271/150‑1 (MQ=255)
                                                                                           gTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                      <  2:2115928/150‑1 (MQ=255)
                                                                                           gTATTCGTACGGTGGGTTTTCACCCTCTCCCTGTGGGAGCGGGCCGGGGTGCGGCCATCAGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTGCTGTCtt                                                      <  2:290105/150‑1 (MQ=255)
                                                                                           gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                      <  2:1268064/150‑1 (MQ=255)
                                                                                            tATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                     <  2:3081630/150‑1 (MQ=255)
                                                                                             aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACTCCCACCAGGGAGGGGGGGCCCCTCGTGGACCAAATTTGTGAGATATATTATTATGtttttt                                                    >  2:384078/1‑148 (MQ=16)
                                                                                              ttCGTATGGTTGTTTTCCCCCCTCGCCCTGGGGGAGAGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                   <  1:1138450/150‑1 (MQ=255)
                                                                                               tCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                  <  2:3867470/150‑1 (MQ=255)
                                                                                                cttatggttggttTTCCCCCTCGCCGTGGGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                 <  1:509520/145‑1 (MQ=255)
                                                                                                cGTAAGGTGGGTTTTCACCCGCGCCCTGGGGGAGCGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                 <  2:506539/150‑1 (MQ=255)
                                                                                                        tgttttaCTCCCACTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                         <  2:1420119/147‑1 (MQ=255)
                                                                                                         gtttttCCCCCCCGCCGGGTGGGTGTGGGGGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                        <  1:3072477/148‑1 (MQ=255)
                                                                                                          gttgacccccgctcccgGTGGGAGAGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                       <  2:710938/143‑1 (MQ=255)
                                                                                                          gttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCACCCGCACCCGAACCCACACCCCCAAGGGTGTGGGGGCCGGTCGAGCACCAAATTTTTGAGTTATATTATTATGTGTTGTTTCAGAAGGACac                                       >  2:2131081/1‑149 (MQ=12)
                                                                                                                 ccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACTGAAAATCg                               <  2:5039025/150‑1 (MQ=255)
                                                                                                                   ctctccgtGTGGGTGCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                             <  2:481167/150‑1 (MQ=255)
                                                                                                                   ctctccgtGGGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                             <  2:452020/150‑1 (MQ=255)
                                                                                                                   ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                             <  2:392034/150‑1 (MQ=255)
                                                                                                                   ctcgccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                             <  2:3815986/150‑1 (MQ=255)
                                                                                                                     ctccctGTGGGAGAGGGCGGGGGGGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                           <  2:4375950/150‑1 (MQ=255)
                                                                                                                     cccccgGTGGGAGTGGGCCGGGGTTAGGCACTCATCGCGCCCGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                           <  2:4375887/148‑1 (MQ=255)
                                                                                                                         ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGGACCCCAACCGAACACAAACCCACCAGGGTGGGGGGGCCCATACAGCACAAAATTTGTGAGATATAATCTTGTTTTTTTTAACATAAGAACCGAAACTACGGGTtttt                        >  2:364128/1‑150 (MQ=255)
                                                                                                                         ctGGGGGAGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                       <  2:1442544/150‑1 (MQ=255)
                                                                                                                          tgtgggtgtgggCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                      <  1:4789850/150‑1 (MQ=255)
                                                                                                                          ggtgGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                      <  2:1863878/149‑1 (MQ=255)
                                                                                                                          ggtgGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                      <  2:1863904/149‑1 (MQ=255)
                                                                                                                           gtgGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                     <  1:2962550/150‑1 (MQ=255)
                                                                                                                             gggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACAACCAAACCCACTCCCACAAGGGGGGGGGGGCCGATCGAGGACCAAATTTGAGAGATAAATTCTTCTGTGGTGTTACCGTAGAAAAGCTAATCGGGGGTTTGGGcc                   >  2:4272237/1‑150 (MQ=255)
                                                                                                                             gggAGAGGGCCGGGGCGAGGGCATCCGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc                   <  2:3090427/150‑1 (MQ=255)
                                                                                                                              ggAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                  <  1:2668174/150‑1 (MQ=255)
                                                                                                                                 gtgggCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa               <  1:4323043/148‑1 (MQ=255)
                                                                                                                                  agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat              <  2:2680899/150‑1 (MQ=255)
                                                                                                                                   gggCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata             <  2:1555550/150‑1 (MQ=255)
                                                                                                                                    ggCGGGGTTGAGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat            <  1:4850490/150‑1 (MQ=255)
                                                                                                                                    ggCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat            <  1:4850895/150‑1 (MQ=255)
                                                                                                                                    ggCGGGGGGGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat            <  1:1870870/150‑1 (MQ=255)
                                                                                                                                         ggttgtgggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg       <  1:3929016/144‑1 (MQ=255)
                                                                                                                                         gggTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg       <  1:3929174/150‑1 (MQ=255)
                                                                                                                                            tGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACACATTTCGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga    <  1:2632185/150‑1 (MQ=255)
                                                                                                                                             gAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa   <  2:853012/150‑1 (MQ=255)
                                                                                                                                              ggggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa  <  1:67546/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAA  >  CP009273/160921‑161212

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: