Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,068 A→C 11.3% 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,0680AC11.3% 258.4 / 8.2 126intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (36/73);  new base C (3/14);  total (39/87)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  CP009273/160920‑161217
                                                                                                                                                     |                                                                                                                                                     
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      <  2:183020/150‑1 (MQ=255)
cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                      <  2:183025/150‑1 (MQ=255)
  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAATAGTATTCGTAAGTTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                    >  1:318333/1‑150 (MQ=255)
  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                    <  1:1458907/150‑1 (MQ=255)
  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                    <  1:1458899/150‑1 (MQ=255)
  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                    <  1:2823482/150‑1 (MQ=255)
  tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCAAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                    <  2:1958662/150‑1 (MQ=255)
         cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                             >  1:875219/1‑150 (MQ=255)
         cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                             >  2:1615904/1‑150 (MQ=255)
         cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCa                                                                                                                                             >  2:1617287/1‑150 (MQ=255)
           aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                           >  1:2584608/1‑150 (MQ=255)
             aatataTGAACGGGGGTTATCCCAAACTTGTCGGGCGGGACGACCGGGAAAATACTGCACCGCCGCGACGTACTAAAATTTATTCGTATGGTTGTTTTTCTCCCTCTCCTTGTGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGtt                                                                                                                                         <  2:809570/145‑1 (MQ=255)
              agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgg                                                                                                                                        >  1:2417273/1‑148 (MQ=255)
                   tGAATGGGCGTTATCCAAACCATGTCTGGCCGGACGACCCGGCAAATCCGGCACCGCCGCGACGGACGAAAATGTATTCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTtcaccct                                                                                                                                   <  2:2729983/150‑1 (MQ=37)
                   tGAACGGGCGTTATCCCAAACATGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTATGGGTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTtcaccct                                                                                                                                   <  2:2729994/150‑1 (MQ=38)
                    gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctc                                                                                                                                  >  1:4047527/1‑150 (MQ=255)
                    gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctc                                                                                                                                  >  1:4047505/1‑150 (MQ=255)
                    gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTATCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaaccac                                                                                                                                  >  2:1757881/1‑148 (MQ=255)
                            tttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                          <  2:2226211/149‑1 (MQ=255)
                            gtaaTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                          <  2:297051/147‑1 (MQ=255)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                          <  2:296926/150‑1 (MQ=255)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTGCCCCCACAACCCAAc                                                                                                                          >  1:1715873/1‑150 (MQ=38)
                            gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCCCACCCaaaa                                                                                                                          >  1:958371/1‑149 (MQ=38)
                             ttATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                         <  2:2332315/150‑1 (MQ=255)
                               aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCCCCCAAAacca                                                                                                                       >  1:2077064/1‑149 (MQ=38)
                               aTCCCAAACATGTCGGGCCGTACGACCCGGAAAATACTGCACCGACGCGACGTACGAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCGTGTGGGTGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCct                                                                                                                       <  2:997357/150‑1 (MQ=38)
                                ccccAAACTTGTCGGGCGGGACGACCCGGAAAAACCTCCACCGACGCCACGTACTAAAAATTATTCTTATGGTGGTTTTTCGCCCTCTCCCTGTGGGAGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                      <  1:307401/149‑1 (MQ=255)
                                 cccAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCAAAcccaca                                                                                                                     >  2:856692/1‑147 (MQ=38)
                                     aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTGGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCt                                                                                                                 <  2:2630880/150‑1 (MQ=255)
                                            cTGGCCGTACGACCCGGAAATTAGTGCTCCGACGCTTGGGACTAAAATTTTGTGGTATGTTTGTTTTTCTCCCTCTCCGTGTGGGAGAGGTCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTTACCCTCTCCCTCAAGGga                                                                                                          <  2:48075/150‑1 (MQ=255)
                                            cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                          <  2:1751830/150‑1 (MQ=255)
                                            cTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                          <  2:4014339/150‑1 (MQ=255)
                                               gCCGGACGACCCGGCAAATACTCCACCGACGCGACGGACGAAAAAGTATGCGTAAGGTTGGTTTTGTCCCTGTCCCTGTGGGTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                       <  1:1810995/150‑1 (MQ=39)
                                                  ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCGAAACCACTACCTCAAGGGAGAgggg                                                                                                    >  2:1088188/1‑150 (MQ=37)
                                                   gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                                   <  2:616875/150‑1 (MQ=255)
                                                     cgacCCGGCAAATACTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                 <  2:62461/150‑1 (MQ=255)
                                                     cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCACACCCAAACCCGCACCCGCCAGGGAGAGGGGTcc                                                                                                 >  2:4321546/1‑150 (MQ=37)
                                                     cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTAAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                 <  1:1267283/150‑1 (MQ=255)
                                                     cgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTAAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                 <  1:2771870/150‑1 (MQ=255)
                                                       ccccGGCAAATACTGCCCCTACGCGACGGACTAAAAATTTTTCGTAAGGTTGGTTTTCCCCCTCACCCTTTGGGAGTGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGa                                                                                               <  1:1965274/149‑1 (MQ=255)
                                                       acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGGGGGCCCCCACCCCAAACCACACCCACAACGGGGGGGGGTCCgg                                                                                               >  1:2318402/1‑149 (MQ=25)
                                                       acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCCAAACCCCACCCCCAAGGGAGTGGGGTcccg                                                                                               >  1:2320251/1‑148 (MQ=25)
                                                         ccGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGTTTTACTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                             <  1:3221345/150‑1 (MQ=255)
                                                                  ccggCACCGACGCGACGGACGAAAACGTAGGCGTACGGTTGGGTTTCGCCCTCTCCCTGTGGCAGAGGGCCGGGGTGCGGGCCTCCGCGCGCACGGTCACCCTCCCCCTAACCCTCTCCCTCACGGGAGCGGGGACCGAGCGAGCACaaa                                                                                    <  2:4393356/147‑1 (MQ=255)
                                                                      cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTCACCCGCCCCCTAACCCACACCCACAAGGGCGAGGGGGCCGGTCGAGCCCAAAttat                                                                                >  2:3924175/1‑148 (MQ=37)
                                                                         cgacgcgacgGACGAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                             <  1:4153515/150‑1 (MQ=255)
                                                                         cgacgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                             <  1:4153574/150‑1 (MQ=255)
                                                                           acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAGCCCACATTTTgggag                                                                           >  2:2395761/1‑150 (MQ=38)
                                                                           acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCACACCCAAAACCACAACCTCAAGGGAGATGGGCACGATCGAGCCCAAATTTTgagag                                                                           >  2:2395754/1‑150 (MQ=25)
                                                                           acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGAGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                           <  2:3762727/150‑1 (MQ=255)
                                                                              cgacgGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGata                                                                        <  1:3702724/150‑1 (MQ=255)
                                                                                acggacgAAATAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                      <  2:3516277/150‑1 (MQ=255)
                                                                                acggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                      <  2:3516402/150‑1 (MQ=255)
                                                                                acggacgAAAAAGTATGCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                      <  2:3516280/150‑1 (MQ=255)
                                                                                acggacgAAAAAGCCATCGTAAGGTTGGGTTTCGCCCTCTCCCGGGGGGAGAGGGCCGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCCCCCTCACCCTCTCCCTCAAGGGAGAGGCGCCCGATCGAGCACAAATTTTGAGAGATAtc                                                                      <  2:3516281/150‑1 (MQ=255)
                                                                                 cgtacgAAAAATTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                     <  2:1242487/147‑1 (MQ=38)
                                                                                 cggactAAAATGTATTCTTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                     <  2:1241692/150‑1 (MQ=255)
                                                                                 cggactAAAAATTTTTCTTATGGTTGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                     <  2:1242435/150‑1 (MQ=38)
                                                                                 cggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCTAACCCGCTCCCACCACGGAGAGGGGCCCGATCGAGCACAAAATTTGAGAGatatat                                                                     >  2:3949544/1‑148 (MQ=37)
                                                                                 cggacgAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                     <  2:3900943/150‑1 (MQ=255)
                                                                                  ggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                    <  2:580069/150‑1 (MQ=255)
                                                                                      gAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCACACCCAAACCCACACCCACAAGGGAGAGGGGACCCATAGAGCACAAAATTTTAGAGATAtcatatta                                                                >  1:2612097/1‑144 (MQ=25)
                                                                                       aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                               <  1:1709669/150‑1 (MQ=255)
                                                                                        aaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCTCACCCGAACCCTCTCCCTCAACGGAGAGGGGACCGATCGAGCCCCAATTTTGAGAGATATATTCTTCTg                                                              >  2:2528384/1‑150 (MQ=38)
                                                                                           aGTATTGGTATGGTTTGTTTTCTCCCTCTCCCGGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAAGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGATGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                           <  1:3906651/150‑1 (MQ=38)
                                                                                               ttCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                       <  2:132866/150‑1 (MQ=255)
                                                                                                 cGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                     <  2:969456/150‑1 (MQ=255)
                                                                                                  ttatgtttggttTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                    <  1:655149/144‑1 (MQ=255)
                                                                                                   tacggttggttTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACag                                                   <  1:3275004/147‑1 (MQ=255)
                                                                                                     tggttggttTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTATCTTGTAACagaa                                                 <  1:2848725/149‑1 (MQ=255)
                                                                                                      tggtggttTTCCTCCCTCTCCGTGTGGGAGAGGGCGGGGGTTAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                 <  1:795763/147‑1 (MQ=255)
                                                                                                       gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTTGACCCTCACCCGAACCCACTCCCGCCAGGGAGAGGGGCCCGCGAGAGCCCCAAGTTTGGGGGTGATCTTCTTCTGTCTTGTAACaggaga                                               >  2:3000434/1‑150 (MQ=21)
                                                                                                       gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCGCACCCAAACCCAAACCCACAAGGGCGAGGGGTCCGATCGAGCCCACAATTTGAGAGATATATTATTATTTATTTTAAAataaga                                               >  1:3839190/1‑150 (MQ=16)
                                                                                                          ggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCGCACCCAAACCCACAACCTCAAGGGAGAGGGGACCGATCGAGCACAAAATTTTAGTGATATCATCTTCTGTCGTGTAACAgaagaaaa                                            >  1:3706570/1‑148 (MQ=25)
                                                                                                                actccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                                       <  1:1700846/149‑1 (MQ=255)
                                                                                                                 ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTACCCCAAAACCACTACCACAAGGGTGAGGGGACCGATAGAGGACAAATTTTGAGAGATATCTTCTTCTGTATTGTAAAAGAATAACAGTAAAt                                      >  2:4232030/1‑150 (MQ=25)
                                                                                                                 ctccctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                      <  2:2284711/150‑1 (MQ=255)
                                                                                                                    cctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                   <  1:2320013/150‑1 (MQ=255)
                                                                                                                       ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                <  2:4047527/150‑1 (MQ=255)
                                                                                                                       ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                <  2:4047505/150‑1 (MQ=255)
                                                                                                                       cgccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                <  2:4047791/148‑1 (MQ=255)
                                                                                                                        tccgtGTGGTTGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                               <  1:3473715/150‑1 (MQ=255)
                                                                                                                        tccgtGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                               <  1:3473688/150‑1 (MQ=255)
                                                                                                                         cccgGTGGGAGCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                              <  1:3146750/150‑1 (MQ=255)
                                                                                                                         cccgGTGGGAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                              <  2:1110507/150‑1 (MQ=255)
                                                                                                                           ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                            <  2:3985831/150‑1 (MQ=255)
                                                                                                                           ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                            <  2:3985944/150‑1 (MQ=255)
                                                                                                                            tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                           <  2:716565/150‑1 (MQ=255)
                                                                                                                               gggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTGACCCTCACCCAAACCCAAACCCGCAAAGGGGAGGGGCCCCATAGAGCACAAAATTTGAGAGATATATTATTGTTTCTTTTAACAAAATCAAAGAAAATAGGGCTTTTTCGcc                        >  2:4349160/1‑150 (MQ=255)
                                                                                                                                ggAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:1243020/150‑1 (MQ=255)
                                                                                                                                   gtgggCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                    <  1:3640766/148‑1 (MQ=255)
                                                                                                                                   gtgggCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                    <  1:3640765/148‑1 (MQ=255)
                                                                                                                                       gcgggggTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                <  2:958371/147‑1 (MQ=255)
                                                                                                                                             gTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag          <  1:856692/150‑1 (MQ=255)
                                                                                                                                                cgggCATCAGCGCGCACGTTCACCCTAACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa       <  1:1983200/149‑1 (MQ=255)
                                                                                                                                                aGGGCATCAGCGCGCACGTTCACCCGCACCCAAAACCACACCCAAAAGGGAGAGGGGACCCATCGAGCACAAAATTTTAGATATATATTAATATGTATTTTAAAATAAAAAAAGAAAATAGGGCTTTTTAGCCAGAAAAATGAGGAGaaa       >  2:309824/1‑150 (MQ=255)
                                                                                                                                                   gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg    <  2:629341/150‑1 (MQ=255)
                                                                                                                                                    cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc   <  1:2038038/150‑1 (MQ=255)
                                                                                                                                                    cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc   <  2:1689866/150‑1 (MQ=255)
                                                                                                                                                     aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  2:2020724/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  CP009273/160920‑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: