Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,068 A→C 30.9% 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,0680AC30.9% 76.7 / 30.6 88intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base C (4/24);  total (35/53)
Fisher's exact test for biased strand distribution p-value = 9.38e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  CP009273/160920‑161214
                                                                                                                                                    |                                                                                                                                                  
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                    <  2:1740975/149‑1 (MQ=255)
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                    <  2:1740835/149‑1 (MQ=255)
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                    <  2:3809371/149‑1 (MQ=255)
 gTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                   >  1:2213951/1‑149 (MQ=255)
  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                  >  2:4453782/1‑149 (MQ=255)
  tGAGGTGAAAAAAGATATGAACGGGCGTTATCCCAAACTTGTCTGGCCGGACCACCCGGAAAATACTGCACCGACGCGACGGACTAAAATTTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGAGGGCATc                                                                                                                                                  <  1:4657917/149‑1 (MQ=255)
  tGAGGTGAAAAAAGAGATGAATGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCTACGGACTAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCGGTGGGTGAGGGCCGGGGTGAGGGCATc                                                                                                                                                  <  1:1185397/149‑1 (MQ=255)
   gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGgcagca                                                                                                                                                 >  1:2486336/1‑149 (MQ=255)
     ggTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgc                                                                                                                                               >  1:2150371/1‑149 (MQ=255)
       tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgc                                                                                                                                             >  2:3913789/1‑149 (MQ=255)
        gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                            >  1:1206467/1‑149 (MQ=255)
        gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcg                                                                                                                                            >  1:4320324/1‑149 (MQ=255)
           aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCAc                                                                                                                                         >  2:4428519/1‑149 (MQ=255)
           aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                         >  2:2636462/1‑149 (MQ=255)
           aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                         >  1:4072400/1‑149 (MQ=255)
            aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                        >  2:3837582/1‑149 (MQ=255)
            aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                        >  1:4340526/1‑149 (MQ=255)
               gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgg                                                                                                                                     >  1:4465567/1‑147 (MQ=255)
               gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgg                                                                                                                                     >  2:1557458/1‑147 (MQ=255)
                agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTCAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACggggg                                                                                                                                    >  1:462732/1‑145 (MQ=255)
                agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACggggg                                                                                                                                    >  1:462790/1‑145 (MQ=255)
                agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggg                                                                                                                                    >  1:3634637/1‑146 (MQ=255)
                 gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGtggac                                                                                                                                   >  2:3745671/1‑145 (MQ=255)
                 gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTgcac                                                                                                                                   >  2:3745719/1‑149 (MQ=255)
                 gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgggggg                                                                                                                                   >  1:2353031/1‑144 (MQ=255)
                 gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggg                                                                                                                                   >  1:1899427/1‑145 (MQ=255)
                 gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggcc                                                                                                                                   >  2:291714/1‑145 (MQ=255)
                  aTGAAAGGGCGTTATCCCAAACATGTCGGGCCGGACGACCCGGAAAATACTGCACCGACGGGACGGACGAAAAAGTTTTCTTATGGTTGTTTTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                  <  2:3309289/149‑1 (MQ=37)
                      tgggggGTTACCCAAAACTTGTCGGGCCGGACGACCGGGAAAATCCGGCACGTCCGGTTGGGATTAAAATTTTTTCTTATGGTTGTTTTTCTCCCTCTCCTTTTGGTTTTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTtcaccctca                                                                                                                              <  2:669489/147‑1 (MQ=255)
                          gCGTTATCCCAACCTTGTGTGGCCGGCCGACCCGGCAAATCCGCCCCCGACGCGCCGTCCTAAAATTTATTCGAATGGTTGGTTTTCTCCCTCCCCCTGTGGTAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTcaccctcaccct                                                                                                                          <  2:3233434/149‑1 (MQ=255)
                          gCGTAATCCCAACCATGACGGGCCGGACGACCCGGAAAATACTGCCCCGACGCGCCGGACGAAAAAGTATTCTTACGTTGGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTcaccctcaccct                                                                                                                          <  2:1728301/149‑1 (MQ=25)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATCCTGCACCGACGCGACGGCCGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGCGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTggccccccacccgca                                                                                                                        >  2:771473/1‑146 (MQ=37)
                             ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGAGGGACTAAAATGTTTTCGTATGGTTGTTTTTCCCCCGCACCCTGGGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                       <  2:3220140/149‑1 (MQ=255)
                             ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGAccccaaccccaaa                                                                                                                       >  2:1203281/1‑148 (MQ=38)
                                 cccAAACTTGGGGGGCCGGACGACCCGGCAAAAACGGCCCCGACGCGAGGGCCTAAAATTTATTCTTATGTTTGTTTTTCCCCCCCTCCCGGGGGTAGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                   <  2:2415735/149‑1 (MQ=255)
                                   cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGGCCCCCACCCCaaaacaaaa                                                                                                                 >  1:4217403/1‑142 (MQ=37)
                                      aCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGAACCACACCCCAAACCACAccca                                                                                                              >  2:2128960/1‑148 (MQ=37)
                                      aCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGGAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTGGGCCCCCACCCCAACCCACAccca                                                                                                              >  1:2832877/1‑148 (MQ=37)
                                         tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGTCCCACACCCCAAACCACAcccacca                                                                                                           >  1:5097214/1‑145 (MQ=37)
                                         tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCCAACCCTCGcccacca                                                                                                           >  2:2354370/1‑145 (MQ=37)
                                            cgggCCGGACGACCCGGCAAATACTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCGCCGTGTGGGTGCGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAggg                                                                                                        <  2:1199505/147‑1 (MQ=39)
                                                ccGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACCCCCGAACCCACTCCCACAAGGgagat                                                                                                    >  2:2087237/1‑148 (MQ=37)
                                                ccGGACGACCCGCAAAATACTGCACCGCCGCGACGGACGAAAATTTATTCTTATGGTTGGTTTACCCCCTCTCCGTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                    <  2:2355670/149‑1 (MQ=37)
                                                ccGCACGACCCGGAAAATACTGCACCGCCGCGCCGGACGAAAAAGTATTCGTATGGTGGGTTTTCACCCTCGCCCGGGGGGCGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                    <  2:4959533/149‑1 (MQ=255)
                                                   gacgccCCGGCAAATACGGCCCCGACGCGACGGACGAAAAAGTAGTGGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggg                                                                                                 <  2:3337477/149‑1 (MQ=38)
                                                   gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTACCTCAAGGGAGAgggg                                                                                                 >  2:3391605/1‑149 (MQ=255)
                                                       acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCGCCCACCCCCGCACCCACACCCGCCAGGGCGAGGGGGCCg                                                                                             >  2:1959711/1‑149 (MQ=25)
                                                       acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAACCCACACCCGAAAGGGGGAGGGGGCCg                                                                                             >  2:1196056/1‑149 (MQ=37)
                                                        cccGGCAAATACTCCACCGACGCGACGGACTAAAATGTATTCGTAGGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGa                                                                                            <  2:3357459/149‑1 (MQ=38)
                                                         ccGGCAACTCCTGCACCGCCGCGCGTGACTAAAATTTATGGTTATGGTTGTTTTGCTCCCTCTCCGTGTGGTTGTGGGCCGGGTTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                           <  2:4163758/149‑1 (MQ=255)
                                                          cgcaaaaTCCTGCCCCGCGGCGCGGGCCCAAAAATTTTTCGTTGGTTTGGTTTTCCCCCTCCCCGTGTGGGTGTGGGGCGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                          <  2:2924473/145‑1 (MQ=255)
                                                            gCAAATACGGCCCCGCCGCGCCGGACTAAAAAGTTTTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                        <  2:3238218/149‑1 (MQ=37)
                                                                ctcctccACCGCCGGGCCGGACGAAAATTTTTTGGTAGGGTTGTTTTGCTCCCCCGCCGGGGGGGTGTGGGCCGGGGGGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcac                                                                                    <  2:187942/143‑1 (MQ=255)
                                                                    ggccccgccgCGCCTGACTAAAATTTATTCGTTTGGTTGTTTTTCCCCCTCCCCCTGGGGGTGCGGGCGGGGTTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                <  1:771064/145‑1 (MQ=255)
                                                                     gCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCAACCAAAAACACACCCGCAAGGGGGAGGGGGCCCCTCGCGCACACAtt                                                                               >  2:3824759/1‑149 (MQ=21)
                                                                       acctcggcgacgGACTAAAAAGTATTGGTTTGGTTGTTTTTATCCCTCTCCGTGTGTGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                             <  2:3200815/143‑1 (MQ=255)
                                                                         cgccgcgccgGACGAAACAGTATGCGTCCGGGTGGTTTTCGCCCGCGCCCGGGGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTGTga                                                                           <  2:4639517/146‑1 (MQ=255)
                                                                               gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCAAACCCACAACCGAGTGGGGCCCGCACGAGCACCAAATTTGAGAGata                                                                     >  2:4292368/1‑149 (MQ=21)
                                                                                        aaaaGTATTCGTACGGTTGGTTTTCTCCCTAGCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                            <  2:3212199/149‑1 (MQ=255)
                                                                                            gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                        <  2:3369684/149‑1 (MQ=255)
                                                                                                 cGTATGTTTGGTTTTCGCCCTCTCCGTGTGGGTGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                   <  2:4616970/149‑1 (MQ=255)
                                                                                                  ttatggttgtttTTCGCCGTCGCCGTGGGGGTGGGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                  <  2:3580655/145‑1 (MQ=255)
                                                                                                   taggtttgtttTTCACCCTCTCCCGGGGGGAGTGGGCGGGGTTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                 <  2:1559560/144‑1 (MQ=255)
                                                                                                    aaGTTTGGTTTACACCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                <  2:4725895/149‑1 (MQ=255)
                                                                                                       gttggttTTCTCCCTCTCCGTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaat                                             <  2:3450857/149‑2 (MQ=255)
                                                                                                       gtggtttttCCCCCTCCCCCTGTGGCTGGGGGCGGGGTTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCGGTCTTGTAACagaag                                             <  2:2837833/144‑1 (MQ=255)
                                                                                                          gtttttCTCCCTCGCCGGGGGGGTGTGGGCGGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                          <  2:4198644/147‑1 (MQ=255)
                                                                                                                ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                    <  2:1289412/149‑1 (MQ=255)
                                                                                                                ctcccgctccctGGGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGaaaa                                    <  2:3419646/149‑1 (MQ=255)
                                                                                                                cgccctctccctGTGGGTGCGGGCCGGGGTGAGGGCCTCTGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                    <  2:1288787/147‑1 (MQ=255)
                                                                                                                 tccctctccgtGGGGGTGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                   <  2:1728976/149‑1 (MQ=255)
                                                                                                                   catccccctGGGGTAGCGGCCGGGGTTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                 <  1:3409889/144‑1 (MQ=255)
                                                                                                                    ctccccctGGGGGTGTGGGGGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                <  1:279967/149‑1 (MQ=255)
                                                                                                                     tctccctGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:760908/149‑1 (MQ=255)
                                                                                                                     tctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAAACCACAACCACCAGGGAGGCGGGCCCCCTCGGGCACCAATTTTTATAGATATAATCTTCTGTCTTGTTACACAAAAACAGCAAATAggg                               >  2:3831386/1‑149 (MQ=12)
                                                                                                                     tcgccctGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                               <  2:760912/146‑1 (MQ=255)
                                                                                                                      ccccctGGGGGAGAGGGCCGGGGTGCGGGCCTCACCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAACGAGCACACATTTAGAGAGATATCTTCTTCGGTCTTGTAACAGAAGAACAGAAAATCGGGc                              <  2:5091250/147‑1 (MQ=255)
                                                                                                                       tccgggtgGGTGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                             <  2:4714405/144‑1 (MQ=255)
                                                                                                                        ccctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAACTCGGGCtt                            <  2:651820/149‑1 (MQ=255)
                                                                                                                         cgggtgGGTGTGGGCGGGGTTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                           <  1:2565405/146‑1 (MQ=255)
                                                                                                                               ggtgggggCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTGCGcc                     <  1:1061565/144‑1 (MQ=255)
                                                                                                                                   ggggCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCCTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTGCGCCTGaa                 <  1:3695759/148‑1 (MQ=255)
                                                                                                                                    gggCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                <  2:733610/149‑1 (MQ=255)
                                                                                                                                    gggCCGGGGCGAGGCCCTCCGCGCGCACGCTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                <  2:3334028/149‑1 (MQ=255)
                                                                                                                                     gccgggggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAATAACAGAAAATCGGGCTTTTGCGCCTGAata               <  2:4494870/145‑1 (MQ=255)
                                                                                                                                            gTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga        <  2:1347228/149‑1 (MQ=255)
                                                                                                                                             tgtgggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag       <  2:3677045/146‑1 (MQ=255)
                                                                                                                                             tgtgggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag       <  2:3677108/146‑1 (MQ=255)
                                                                                                                                             tgcgggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag       <  2:940351/146‑1 (MQ=255)
                                                                                                                                                 ggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa   <  2:462732/149‑1 (MQ=255)
                                                                                                                                                  gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  2:4453894/149‑1 (MQ=255)
                                                                                                                                                  gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa  <  2:807426/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAA  >  CP009273/160920‑161214

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: