Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 161,068 A→C 23.5% intergenic (+47/‑149) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273161,0680AC23.5% 124.1 / 20.1 117intergenic (+47/‑149)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (19/69);  new base C (7/22);  total (26/91)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

GTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  CP009273/160921‑161217
                                                                                                                                                    |                                                                                                                                                     
gTGCGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                     <  2:850459/150‑1 (MQ=255)
gTGAGGTGCAAAAAGAGATGAACGGGCGTTATCCCACACAGGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATc                                                                                                                                                     <  2:4631767/150‑1 (MQ=255)
  gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                                   >  2:1158257/1‑150 (MQ=255)
  gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                                   >  2:2520510/1‑150 (MQ=255)
    ggTTCAAAAAGTTATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                                                 <  2:1444744/150‑1 (MQ=38)
       gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                                              >  2:3552684/1‑150 (MQ=255)
       gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                                              >  2:3552490/1‑150 (MQ=255)
          aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACg                                                                                                                                           >  1:2504644/1‑150 (MQ=255)
           aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGt                                                                                                                                          >  1:1370075/1‑150 (MQ=255)
             atataTGAACGGGCGTTATCCCAAACTTGTCTGGCCGGACGACCCGGAAAATACGGCACCGACGCGACGGACGAAAAAGTATTCTTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTtc                                                                                                                                        <  2:3973678/146‑1 (MQ=37)
             agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGtgc                                                                                                                                        >  2:3282262/1‑148 (MQ=255)
                taTGAATGGGCGTTATCCCAAAATTGTGTGGCGGGCCGACCCGGAAAATACTGCACCGACGCGACGGACTAAAATTTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                     <  2:2587458/149‑1 (MQ=25)
                    aaaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccccca                                                                                                                                 >  2:2830759/1‑150 (MQ=255)
                        ggCTTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAATTATTCTTAAGGTTTGTTTTCTCCCTCGCCCTGGGGTTGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTcaccctcaccct                                                                                                                             <  2:471079/150‑1 (MQ=38)
                         gCGTTATCCCAAACATGTCTGGCCGGACCACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTa                                                                                                                            <  2:1079391/150‑1 (MQ=255)
                           gTTATCCCAAACATGTTTGGCGGTACGACCCGTAAAATCCTGCACCGCCGCTCGGGACTAAAATTTTTGCTTATGGTTGTTTTTCTCCCTCTCCCTGTGGTTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                          <  2:473193/150‑1 (MQ=19)
                           gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTGGGCCCCCACCCAAAc                                                                                                                          >  1:2561945/1‑150 (MQ=38)
                           gTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGTACTAAAAATTATTCGTATGTTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                          <  2:4404315/150‑1 (MQ=38)
                           gTTAACCCAAACTTTACTGGCGGCACGCCCCGCAAATTACTGCACGGAGGGGAGGCAGTAAAATTATTTCTTACGGTTGTTTTTCACCCTCTCCCTGTGGGATTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                          <  2:473061/150‑1 (MQ=255)
                            ttATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCCCCCAAAcc                                                                                                                         >  1:3920737/1‑150 (MQ=38)
                                   aaaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTGGCCCCGCCCCCAAACCCACAccc                                                                                                                  >  1:3289158/1‑150 (MQ=37)
                                    aaCATGTCTGGCCGTACGACCCGGCAAATACTTCACCCACGCGACGGACTAAAAAGTATTCTTAAGTTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCt                                                                                                                 <  2:1488630/150‑1 (MQ=38)
                                    aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGAACCGCACCCAAACCCACAcccc                                                                                                                 >  2:2424539/1‑149 (MQ=37)
                                         gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCACACCCAAACCCACACCCACaacg                                                                                                            >  2:612980/1‑148 (MQ=38)
                                         gTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCCCCCACACCCAAAACCACACCCACAAgg                                                                                                            >  2:3751019/1‑150 (MQ=37)
                                          tCTGGCCGTACGACCCGGCAATTACTGCACCGACGCGACGGACGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCACTCCCTGTGGTAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAggg                                                                                                           <  2:1519331/150‑1 (MQ=255)
                                          tCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATCCGTAAGGTTGGTTTTCCCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAggg                                                                                                           <  2:1519289/150‑1 (MQ=255)
                                          tCTGGCCGGACGACCCGCCAATTCCTGCACCCACGGTACGGACGAAAAAGTTTTCGTATGGTTGGTTTTCTCCCTCACCGTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAggg                                                                                                           <  1:1148072/150‑1 (MQ=255)
                                           cgggCCGGACGACCCGGAAAATCCTGCACCGACGCGACGTACTAAAAAGTTTTCGTAAGGTTGGTTTTCCCCCTCTCCCTGTGGGTTCGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                          <  1:2188124/148‑1 (MQ=255)
                                            tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCACGTTGACCCGCCCCCAAAACCCCAACCAAAAGGgag                                                                                                         >  1:4214284/1‑150 (MQ=37)
                                            tGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCGAACCCACAACCGCAAgggcg                                                                                                         >  1:4214302/1‑148 (MQ=37)
                                             ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCTGGGTGAGGGCATCAGCGCGCACGTTCGCCCACAACCAAACCCACACCCACaacgtcag                                                                                                        >  2:3638127/1‑144 (MQ=37)
                                              gccggccgCCCCGGCAAATACTGCCCCGACGCGCGGGACTAAAATTTATTCTTATGGTTGTTTTTCACCCCCTCCCGGGGGGTGTGGGGCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                       <  1:705830/150‑1 (MQ=255)
                                               ccGGCCGCCCCGGCAAATCCTGCCCCGACGCGACGGACGAAAAATTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGATAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgg                                                                                                      <  2:2061323/150‑1 (MQ=38)
                                                    cgccccGGAAAATACTGCCCCGCCGCCACGGCCCAAAAATTATTCTAAAGGTTGGTTTACGCCCCCTCCCTGTGGGAGAGGGCCGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                 <  2:72752/147‑1 (MQ=255)
                                                    cgacCCGGAAAATACGGCACCGACGCGACGGACGAAAAATTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGCGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                 <  2:220236/150‑1 (MQ=38)
                                                     gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCTACTCCCTCAAGGGAGAGGGGACCg                                                                                                >  1:143209/1‑150 (MQ=255)
                                                     gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGGGCCCACCCCCCCACCCCCACCCACCAGGGGGCGGGGGCCg                                                                                                >  1:2907892/1‑150 (MQ=25)
                                                     gacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCCCCCCAACACACAACCACAAGGGCGTGGGGGCCg                                                                                                >  1:1307791/1‑150 (MQ=25)
                                                      ccccGGAAAATCCGGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGTTTTTCCCCCTCTCCCTGTGGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGa                                                                                               <  2:1830063/149‑1 (MQ=255)
                                                      acCCGGCAAACACGGCACCGCCGCGCCGGACGAAAAATTATTCGAACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGa                                                                                               <  2:3923802/150‑1 (MQ=37)
                                                        ccGGCAAATACTGCCCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                             <  2:1036636/150‑1 (MQ=255)
                                                          ggCAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGATAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                           <  2:3018745/150‑1 (MQ=255)
                                                            cAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCCCCCCAACCCACACCCGCCACGGAGAGGGGGCCGCTCGAGc                                                                                         >  2:508232/1‑150 (MQ=25)
                                                              aaTCCTGCACCGACGCGCCGGACTAAAACTTATTCGTACGGTTGTTTTACCCCCTCGCCCTGGGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcac                                                                                       <  1:3775166/150‑1 (MQ=255)
                                                               acacggCACCGACGCGAGGTACGAAAAAGTATTCTTACGGTTGGTTTTCTCCCTCTCCCGGGGGGTGAGGGCCGGGGGGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcaca                                                                                      <  2:667049/145‑1 (MQ=255)
                                                               aTACTTCCCTGACGCGACGGACTAAAAAGTTTTCGTAAGGTTGTTTTTGTCCCTCTCCCTGTGGTAGGGGGCCGGGTTTTGGGACTCCGCGCGCACGTTCTCCCTCACCCTAACCCTCTCCCTCAAGGGATAGGGGTACGATCGAGcaca                                                                                      <  2:667020/150‑1 (MQ=255)
                                                               aTACTGCCCCCACGCCACGGACCAAAAAGTATTCGTAAGGTTGTTTTTCTCCCTAGCCCAGTGGGATAGGCCCGGGGTTAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcaca                                                                                      <  2:2549459/150‑1 (MQ=255)
                                                               aTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcaca                                                                                      <  2:2601129/150‑1 (MQ=255)
                                                                   tGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCCCCCGCACCCAAACCCAAACCCACAAGGGGGAGGGGGCCGGTCGAGgacacaag                                                                                  >  2:884299/1‑146 (MQ=25)
                                                                    gccccGAGGCTAGGGACTAAAATTTTTTCTTATGTTTGGTTTTCGCCCCCCCCTTGGGGTATGGGGCGGGGTTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                 <  1:4131339/147‑1 (MQ=255)
                                                                    gccccGACGCGCCGGACTAAAATTTTTTGGTATGGTTGGTTTTCTCCCCCCCCGGGGGGGTGCGGGCGGGGGTGCGGCCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                 <  1:4131457/147‑1 (MQ=255)
                                                                     ccccgccgCGACGGCCGAAAAATTATGGTTATGGTTGGTTTTCGCCCACTCCCTGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                <  1:956237/148‑1 (MQ=255)
                                                                         gacgcgcggGACTAAAATGTTTTCGTATGGTTGTTTTTCCCCCTCTCCCTGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                            <  1:2579948/150‑1 (MQ=255)
                                                                            gcgacgGACTAAAATGTTTTCGTATGGTTGGTTTTCGCCCCCGCCGTGTGGGTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                         <  1:2799960/150‑1 (MQ=255)
                                                                              gacggactAAAAAGATTTCGTAGGTTTGGTTTTCTCCATCTCCCTGTGGGATCGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                                       <  1:4088766/150‑1 (MQ=255)
                                                                               ccggccgAAAACTTATAGGTACGTTTGTTTTCCTCCCCCTCCGTGTGGTTGTGGGCGGGGTTTCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                      <  2:774717/149‑1 (MQ=255)
                                                                                    cgAAATGGTATTCGTTCGGGTTGTGTTTTCCGTGGCCCTGTGGGATTGGGCGGGGGTGAGGGCCTCAGCGCGCATGTTCACCCTCCCCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGTTCGAGCACAAATTTTGAGAGATAtcttctt                                                                 <  1:576112/150‑1 (MQ=255)
                                                                                     gAAAAAGTATTCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttc                                                                <  2:1476918/150‑1 (MQ=255)
                                                                                      aaaaaGTATTCGTAAGGTTGGTTTTCCCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                               <  2:1677868/150‑1 (MQ=255)
                                                                                       aaatttATTCGTATGGTTGGTTTTCTCCCTCTCCCGGTGGGAGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTg                                                              <  1:2269479/145‑1 (MQ=255)
                                                                                           gTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGGGGGAGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:3063444/150‑1 (MQ=255)
                                                                                           gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:668001/150‑1 (MQ=255)
                                                                                           gTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:1020305/150‑1 (MQ=255)
                                                                                           gTATTCGTAAGGTTGGTTTACTCCCACTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:668020/150‑1 (MQ=255)
                                                                                            gttttCTTATGGTTGTTTTTCTCCCTCTCCCTGGGGGATCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  1:3063446/147‑1 (MQ=38)
                                                                                                    acgtttgtttTTCCCCCACCCCCTGGGGTTTTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                  <  2:24540/146‑1 (MQ=255)
                                                                                                       tttgtttTTCCCCCTCACCTTGTGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                               <  1:259336/149‑1 (MQ=255)
                                                                                                           gttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCCCACCCAAAACCACTCCCACCACGGGGAGGGGGCCGATCGGGCACAAATTTTGTGCTTTATATTGTTAAGTGTTGTAAAAGAATAACac                                           >  2:2964286/1‑149 (MQ=12)
                                                                                                              ttCCCCCTCTCCGTGTGGGTGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATTTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                        <  1:4286739/150‑1 (MQ=255)
                                                                                                                  ccctccccctGTGGTTTTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  2:1169401/150‑1 (MQ=255)
                                                                                                                    ctctccctGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                  <  2:3406521/150‑1 (MQ=255)
                                                                                                                     tctccctGTGGTAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGc                                 <  1:287887/150‑1 (MQ=255)
                                                                                                                     tccccgtgtgGGAGTGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGc                                 <  1:288003/144‑1 (MQ=255)
                                                                                                                      ccccgtgtgGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                <  1:4072410/145‑1 (MQ=255)
                                                                                                                       tccctGGGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                               <  1:1064196/150‑1 (MQ=255)
                                                                                                                        ccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGAGATCTTCTTCTGTCTTGTAACAGAAGAACCGAAAATCGGGCttt                              <  2:2918812/150‑1 (MQ=255)
                                                                                                                        ccctGGGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACACGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                              <  2:2918802/150‑1 (MQ=255)
                                                                                                                         ccggtgGGAGCGGGCGGGGTTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                             <  2:1932658/147‑1 (MQ=255)
                                                                                                                          ctGTGGGTGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                            <  1:3402903/150‑1 (MQ=255)
                                                                                                                          cggggggAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                            <  2:3735575/146‑1 (MQ=255)
                                                                                                                            gtgGGAGTGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                          <  2:794348/150‑1 (MQ=255)
                                                                                                                             tgGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                         <  2:4341993/150‑1 (MQ=255)
                                                                                                                             tgGGAGAGGGCCGGGGGGCGGGCCTCAGCGCGCACCTTCACCCTCACCCTACCCCTCCCCCTCAAGGGAGAGGGCACCGATCGAGCACAAATTTTGAGAGATATCTTCGTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcgc                         <  2:4342524/150‑1 (MQ=255)
                                                                                                                                gagaGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCACCCTCACCCGAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAACCTCGGGCTTTTGCGCCTg                      <  2:1113948/150‑1 (MQ=255)
                                                                                                                                    gggCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                  <  2:2131095/150‑1 (MQ=255)
                                                                                                                                    gggCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata                  <  2:3659366/150‑1 (MQ=255)
                                                                                                                                      gCCGGTGTGAGTGCATCAGCGCGCACGTTCACCCTCTCACTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                <  1:4160974/150‑1 (MQ=255)
                                                                                                                                      gCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                <  2:4088711/150‑1 (MQ=255)
                                                                                                                                        gggggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc              <  1:2132113/149‑1 (MQ=255)
                                                                                                                                          gggTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg            <  1:3564089/150‑1 (MQ=255)
                                                                                                                                          gggTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACATAAAATCGGGCTTTTGCGCCTGAATATTGCgg            <  2:1064192/150‑1 (MQ=255)
                                                                                                                                           ggTGTGGGCATCAGCGCGCACGTTCACCATCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGTTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGTGga           <  2:3648713/150‑1 (MQ=255)
                                                                                                                                            gTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgag          <  2:4503499/150‑1 (MQ=255)
                                                                                                                                             tgtgggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGgaga         <  2:1147464/147‑1 (MQ=255)
                                                                                                                                              gcgggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa        <  2:2807302/148‑1 (MQ=255)
                                                                                                                                               aGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaa       <  1:4101321/150‑1 (MQ=255)
                                                                                                                                                gggCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  1:1236660/150‑1 (MQ=255)
                                                                                                                                                gggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa      <  1:3112617/150‑1 (MQ=255)
                                                                                                                                                    aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAATTTGAGAGATATCTTCTTCGGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  2:2510604/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: