Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,068 A→C 17.4% 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,0680AC17.4% 107.6 / 13.7 87intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (33/38);  new base C (3/13);  total (36/51)
Fisher's exact test for biased strand distribution p-value = 5.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  CP009273/160921‑161217
                                                                                                                                                   |                                                                                                                                                      
gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTc                                                                                                                                                      >  1:1532806/1‑150 (MQ=255)
gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                      >  1:1532781/1‑150 (MQ=255)
gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                      >  1:2808784/1‑150 (MQ=255)
  gAGTTGAAAAAAGTGATGAAAGGGCGTTATCCAAAACAGGTCTGCCCGGACGACCCGGCAAATAGTGCACAGACGCGACGGACGAAAATTTATTCGTATGGTTGGTTTTCTCCCTCTCCTTGTGGGTGAGGGCCGGGGTTAGGGCATCAg                                                                                                                                                    <  2:882629/150‑1 (MQ=255)
  gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAACAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAg                                                                                                                                                    >  1:3984749/1‑150 (MQ=255)
  gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                                    >  1:3984708/1‑150 (MQ=255)
  gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                                    <  2:3923247/150‑1 (MQ=255)
    tgTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                                                  >  2:3730097/2‑150 (MQ=255)
     gTGCAATAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGTCGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                                 >  1:2180885/1‑150 (MQ=255)
         aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                             >  1:704361/1‑150 (MQ=255)
          aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                            >  1:3048403/1‑150 (MQ=255)
               agaTGACAGGGGTTTAGCCAAACCATGTGGGGCCGGCCGCCCCGTCAATTACGGCCCCGACGCGCGGGACTAAAAAGTTTTCTTATGGTTGGTTTTCTCCCTCTCCGTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcac                                                                                                                                       <  2:3540135/150‑1 (MQ=255)
               agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggac                                                                                                                                       >  1:2211188/1‑146 (MQ=255)
                gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTGAGGTGGGTTTTCTCACTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCGGCGCGCACggggtcg                                                                                                                                      >  2:3227884/1‑144 (MQ=255)
                gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGttggcc                                                                                                                                      >  2:3227829/1‑146 (MQ=255)
                 tttaacggggGTTATCCAAAACTTGTGGGGCCGGACGACCGGGAAAATACTGCCCCGCCGCGCGGGACTAAAAATTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGGGGTTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTtcaccc                                                                                                                                     <  2:3660952/144‑1 (MQ=255)
                 aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                                     >  2:736806/1‑150 (MQ=255)
                 aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaccc                                                                                                                                     >  2:2609658/1‑150 (MQ=255)
                 aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaacc                                                                                                                                     >  2:664538/1‑150 (MQ=255)
                  tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgcaccca                                                                                                                                    >  2:1137097/1‑149 (MQ=255)
                   gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggacccgc                                                                                                                                   >  2:3314716/1‑148 (MQ=255)
                   gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGtcacccac                                                                                                                                   >  2:2263109/1‑148 (MQ=255)
                     atgggCTTTATCCCAAACAGGTCTGGCCGGACCACCCGGAAAATACGGCACCGACGCGACGGACTAAAAAGTATTCGTAAGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcaccctcac                                                                                                                                 <  2:1049168/148‑1 (MQ=37)
                        ggCGTTAGCCCACACTTGTCTGGCCGGACGACCGGGAAAATACTGCACCGACGCTCCGTACTAAAATGTTTTCGTACGTTTGGTTTTCTCCCTCTCCTTGGGGGTGTGGGCCGGGGTTTGGGCATCAGCGCGCACGTTcaccctcaccct                                                                                                                              <  2:4361021/150‑1 (MQ=255)
                          cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGAGCCACACCCaaa                                                                                                                            >  2:925400/1‑150 (MQ=38)
                          cGTTACCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGTACTAAAAATTATTCGTAAGGTTGTTTTTCTCCCTCTCCCGGGGGTTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                            <  2:92818/150‑1 (MQ=255)
                            ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACCCCCCAAca                                                                                                                          >  2:3206047/1‑149 (MQ=255)
                               tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                       <  2:2928498/150‑1 (MQ=255)
                               tCCAAAACTTGTGTGGCCGTACGACCCGGAAATTAGTCCACCGACGCGACGTACTAAAATTTATTGTTATGTTTGTTTTTCTCCCTCGCCCTGGGGTTTAGGGCCGGTGTGCGGGAATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                       <  2:303423/150‑1 (MQ=255)
                               ccccAAACTTTTCTGGCCGGACTACCCGGCAAATCCTGCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                       <  2:303428/149‑1 (MQ=38)
                                   aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCGAACCCAATccc                                                                                                                   >  2:320079/1‑150 (MQ=38)
                                   aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCCCCCGAACCCACTccc                                                                                                                   >  2:2518406/1‑150 (MQ=255)
                                    aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCGCACCCGAAACCACAcccg                                                                                                                  >  2:948922/1‑149 (MQ=38)
                                    aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACAACCAAACCCCCAcccc                                                                                                                  >  2:949082/1‑149 (MQ=37)
                                      cttGTCTGGCCGCACGCCCCGGCAAATACTGCACCTACGCGACGGACGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                                <  2:4389777/148‑1 (MQ=38)
                                      cATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTGACCCTCACCCTAAACCACAcccacc                                                                                                                >  2:2884921/1‑147 (MQ=38)
                                      cATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCCCACCCCAAACCACAcccacc                                                                                                                >  2:2466792/1‑147 (MQ=37)
                                        tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCGCACCCAAACCCGCGcccgccac                                                                                                              >  2:3376635/1‑145 (MQ=38)
                                             ggcggcacgacCCGCAAAATACTGCCCCGCCGCGACGCACTAAAATTTATTCTTATGTTGGTTTTTCTCCATCTCCCTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                         <  2:836916/144‑1 (MQ=255)
                                             ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                         <  2:4126449/150‑1 (MQ=255)
                                             ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                         <  2:4126294/150‑1 (MQ=255)
                                             ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                                         <  2:4126317/150‑1 (MQ=255)
                                               ccGGACGACCCGGCAAATACTGCCCCGACGCGAGGGACTAAAAATTATTCTTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGTGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                       <  2:1946230/150‑1 (MQ=38)
                                                 ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCACACCCTAACCCAAAGCCGCAAGGGAGAgggg                                                                                                     >  2:2846760/1‑150 (MQ=38)
                                                   acgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGACCCCCCCCCAAAACCACAACCCCAAGGGAgggggggc                                                                                                   >  1:1484618/1‑148 (MQ=25)
                                                    cgacCCGGAAAATACTGCCCCGACGCGCCGGACTAAAATGTATTCGTATGGTGGTTTTTCCCCCTCCCCCGGTGGGATCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                  <  1:3853670/150‑1 (MQ=255)
                                                        ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCGAACCCCCTCCCCCAAGGGGGCggggggcgctg                                                                                              >  2:4278955/1‑143 (MQ=37)
                                                         cgcaaaaTACTGCCCCGACGCGCCGGACTAAAAAGTATACTTCAGGTTGTTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                             <  2:3137300/146‑1 (MQ=37)
                                                             aaaTACTCCCCCGCCGCGACGGACGAAAAATTATTCGTACGGTTGGTTTCCCCCCACTCCCTGTGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                        <  2:1375861/150‑1 (MQ=255)
                                                                  cTGCACCGACGCTCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCATCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                   <  2:331003/150‑1 (MQ=255)
                                                                   ggCACCGAGGCGACGGACGAAAAAGATTACGTCTGTTTGGTTTAACCACACGATGTGTGGTAGTGGGCTGGGGATTGGGACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                  <  2:1200551/149‑1 (MQ=255)
                                                                     cACCGCCGCGAGGCACTAAAATTTTTTCGTATGGTTGTTTTTCTCCCTCTCGCGGTGGGAGGGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  1:1765755/150‑1 (MQ=255)
                                                                      aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                               <  2:4252329/150‑1 (MQ=255)
                                                                          aggcgccgGACGAAAAATTATTGGTAAGGTTGTTTTTCTCCCTCTCCTTGTGGTAGTGGGCCGGGTTTAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                           <  2:2203225/148‑1 (MQ=255)
                                                                           cgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGAGCCCCACCCCAAACCACAACCACCAGGGGGAGGGGGGCGGACGAGCCCAAAAGTTgaggga                                                                          >  1:4549361/1‑150 (MQ=18)
                                                                             cgacgGACGAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGATTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGata                                                                        <  2:2570418/150‑1 (MQ=39)
                                                                                 ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTTGACCCACACCCCAAACCCCAACCACAAGGGGGAGGGGACCGCTCGAGCACAAAAATTTATAGATAtatt                                                                    >  2:1767203/1‑150 (MQ=21)
                                                                                  gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGCCCCACACCCAAACCCACAACCACAATGGCGAGGGGACCGAACGAGCACAAAAATTGAGTTATAtatta                                                                   >  1:478596/1‑149 (MQ=21)
                                                                                  gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACACGCACCCAAAACCACACCCACCAAGGCGAGGGGGCCCAGCGAGCACAAAATTTGAGAGATttcttc                                                                   >  2:4431998/1‑150 (MQ=21)
                                                                                  gacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGGCCCAAACCCAAAAACACACCCACAAATGAGAGGGGCCCCCACGAGCACACAAATTGATAGatatataa                                                                   >  1:478595/1‑146 (MQ=12)
                                                                                           tttttCTTAGGGTTGTTTTTCCCCCCCTCCCTGTGGTTGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:1326392/147‑1 (MQ=255)
                                                                                            tATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTCACCCGAACCCAATACCACAATGGGGGGGGGACCCATAGAGGCCCAAATTTTTGAGATATATTATTCTGTCtgg                                                         >  2:3025940/1‑148 (MQ=16)
                                                                                             aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                        <  2:3828904/150‑1 (MQ=255)
                                                                                                 ttagggttggttTTCCCCCTCGCCGTGTGGGTGGGGGGGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                    <  1:4335182/146‑1 (MQ=255)
                                                                                                 gTATGGTTGTTTTTCCCCCTCCCCCGGTGGGTGTGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                    <  1:4335180/150‑1 (MQ=255)
                                                                                                      gttggttTCCTCCCGGCCCGTGTGGGTGCGGGCGGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                               <  1:1801445/150‑1 (MQ=255)
                                                                                                        tggttTTCCCCCCCTCCCGGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                             <  2:323212/150‑1 (MQ=255)
                                                                                                                tccctctccctGGGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                     <  2:1987483/150‑1 (MQ=255)
                                                                                                                 ccctcgcccgGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  2:51845/150‑1 (MQ=255)
                                                                                                                  cagctccgtGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  1:658367/147‑1 (MQ=255)
                                                                                                                     ctccctGTGGTCGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                <  2:3606276/150‑1 (MQ=255)
                                                                                                                      tccctGTGGGCGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                               <  1:1921504/150‑1 (MQ=255)
                                                                                                                        cgtgtgGGAGTGGCCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                             <  1:3954712/148‑1 (MQ=255)
                                                                                                                        cctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                             <  1:2603525/150‑1 (MQ=255)
                                                                                                                        cctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTACTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                             <  2:3470301/150‑1 (MQ=255)
                                                                                                                         ctGTGGTAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                            <  2:621394/150‑1 (MQ=255)
                                                                                                                              ggAGAGGGCGGGGTTCTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  2:1977252/150‑1 (MQ=255)
                                                                                                                                agtgggCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                     <  2:4261250/147‑1 (MQ=255)
                                                                                                                                agcgggCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                     <  1:2518406/147‑1 (MQ=255)
                                                                                                                                agagGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                     <  2:1452387/150‑1 (MQ=255)
                                                                                                                                     cccGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                <  2:1468206/149‑1 (MQ=255)
                                                                                                                                           gtgtGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag          <  1:925400/150‑1 (MQ=255)
                                                                                                                                           gtgtGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag          <  2:4021162/150‑1 (MQ=255)
                                                                                                                                               gggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  1:3074461/150‑1 (MQ=255)
                                                                                                                                               gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  1:3074465/150‑1 (MQ=255)
                                                                                                                                               gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  1:4011894/150‑1 (MQ=255)
                                                                                                                                                ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa     <  2:3999821/150‑1 (MQ=255)
                                                                                                                                                   aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  1:2787864/150‑1 (MQ=255)
                                                                                                                                                   aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  2:1484618/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  CP009273/160921‑161217

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: