Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 456,588 G→T 82.6% D749Y (GAC→TAC)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273456,5880GT82.6% 790.6 / 182.4 337D749Y (GAC→TAC) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (25/34);  new base T (137/141);  total (162/175)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 340 total aligned reads displayed.

CCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGC  >  CP009273/456439‑456737
                                                                                                                                                     |                                                                                                                                                     
ccGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTg                                                                                                                                                       >  1:4181883/1‑150 (MQ=255)
   aTGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACa                                                                                                                                                    <  2:1324255/150‑1 (MQ=255)
       ggCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGt                                                                                                                                                >  2:636052/1‑150 (MQ=255)
           aTGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAAtt                                                                                                                                            <  2:2411423/150‑1 (MQ=255)
            tGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTg                                                                                                                                           >  1:1575966/1‑150 (MQ=255)
              aCCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAAGTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGcc                                                                                                                                         <  2:1435281/150‑1 (MQ=255)
               ccGGTGAGATCACTCTTCGTGGTCATGTACTTCCGATCGGTTGTTTGAAAGAAAATCTCCTGGCATCGCATAGCTGCGTGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCg                                                                                                                                        <  2:185471/150‑1 (MQ=255)
               ccGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCg                                                                                                                                        <  2:175874/150‑1 (MQ=255)
               ccGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCg                                                                                                                                        <  2:1150424/150‑1 (MQ=255)
                cGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGa                                                                                                                                       <  1:1704482/150‑1 (MQ=255)
                 ggTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGAt                                                                                                                                      <  2:3593457/150‑1 (MQ=255)
                 ggTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGAt                                                                                                                                      <  2:5094433/150‑1 (MQ=255)
                 ggTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGAt                                                                                                                                      <  2:4545690/150‑1 (MQ=255)
                    gagaTCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTg                                                                                                                                   <  1:5381403/150‑1 (MQ=255)
                     agaTCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTgg                                                                                                                                  >  2:3704834/1‑150 (MQ=255)
                       aTCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGAc                                                                                                                                <  2:3859666/150‑1 (MQ=255)
                           ctctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGTGATTCCTTACATCGTAATTGCCGATATGGACATTc                                                                                                                            >  2:3507093/1‑150 (MQ=255)
                            tctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCa                                                                                                                           <  1:1918571/150‑1 (MQ=255)
                             ctGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCAt                                                                                                                          >  1:601712/1‑150 (MQ=255)
                                     tCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGaa                                                                                                                  >  1:77210/1‑150 (MQ=255)
                                      cAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAg                                                                                                                 >  1:409065/1‑150 (MQ=255)
                                      cAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAg                                                                                                                 >  1:1024898/1‑150 (MQ=255)
                                       atgTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTTGCATCGCATCGCGGCGGGATTAAAACAGTGCTTATTCCGTTAGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAgc                                                                                                                <  2:4704925/148‑1 (MQ=255)
                                        ggTACTTCCGATCGGTGGTTTGAAAGAAAAACTCCTGTCAGCGCATCGCGGCGGGATTAATTCAGTGCTAATTCCGTTCGAAAATTAACGCGATCTGGAAGAGATTCCTTACATCTTTATTGCCGATCTGGACATTCATCCTGTGAAgcg                                                                                                               <  1:135670/150‑1 (MQ=255)
                                        ggTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAgcg                                                                                                               <  2:711805/150‑1 (MQ=255)
                                           aCTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCAt                                                                                                            <  2:2827521/150‑1 (MQ=255)
                                           aCTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCAt                                                                                                            <  1:146408/150‑1 (MQ=255)
                                             tGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTg                                                                                                          <  2:1945771/150‑1 (MQ=255)
                                                 gATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGa                                                                                                      <  1:3161740/150‑1 (MQ=255)
                                                   tCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAg                                                                                                    >  2:1502554/1‑150 (MQ=255)
                                                     ggtggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGtt                                                                                                  <  2:297685/150‑1 (MQ=255)
                                                      gtggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTc                                                                                                 >  2:1008449/1‑150 (MQ=255)
                                                       tggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCt                                                                                                <  1:3288413/150‑1 (MQ=255)
                                                        ggtTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTg                                                                                               >  2:5210664/1‑150 (MQ=255)
                                                          tttGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGAc                                                                                             <  2:5519598/150‑1 (MQ=255)
                                                           ttGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGAct                                                                                            >  1:2035640/1‑150 (MQ=255)
                                                           ttGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGAct                                                                                            >  2:3741398/1‑150 (MQ=255)
                                                             gtaagataATCTCCTGGCAGCGCATCGCGGCGGTATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGTATGTGATTCCTTACAACGTAATTTCCGTTCTTGACATTCATCCTGTTAAGAGCTTTGAGGAAGTTCTGActct                                                                                          <  1:3982762/148‑1 (MQ=255)
                                                             gaaagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGActct                                                                                          <  2:1221553/150‑1 (MQ=255)
                                                               aagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTgg                                                                                        >  1:1169552/1‑150 (MQ=255)
                                                               aagaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTgg                                                                                        <  1:4771188/150‑1 (MQ=255)
                                                                agaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGgc                                                                                       <  1:3388417/150‑1 (MQ=255)
                                                                 gaaaAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGgcg                                                                                      >  2:1431825/1‑150 (MQ=255)
                                                                    aaaCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTg                                                                                   >  2:5461021/1‑150 (MQ=255)
                                                                    aaaCTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTg                                                                                   <  1:2508032/150‑1 (MQ=255)
                                                                           tGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATg                                                                            <  1:5302214/150‑1 (MQ=255)
                                                                           tGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATg                                                                            >  1:3324685/1‑150 (MQ=255)
                                                                            ggCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGACGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGa                                                                           >  2:3897331/1‑150 (MQ=255)
                                                                            ggCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGa                                                                           <  1:3476788/150‑1 (MQ=255)
                                                                            ggCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTAGAAAATAATCGCGTTCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCATAATGa                                                                           <  2:3460152/150‑1 (MQ=255)
                                                                             gCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGaa                                                                          <  1:3792967/150‑1 (MQ=255)
                                                                              cAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAAc                                                                         >  2:1113707/1‑150 (MQ=255)
                                                                              cAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAAc                                                                         <  1:394784/150‑1 (MQ=255)
                                                                                 cgcATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGt                                                                      >  2:3321823/1‑150 (MQ=255)
                                                                                   cATCGCGGCGGGATTAAAACCGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCt                                                                    <  1:5244156/150‑1 (MQ=255)
                                                                                   cATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGACGTTCTGACTCTGGCGCTGCAAAATGAACCGTCt                                                                    >  2:3951885/1‑150 (MQ=255)
                                                                                   cATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCt                                                                    <  2:5366306/150‑1 (MQ=255)
                                                                                     tCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTgg                                                                  >  2:2751878/1‑150 (MQ=255)
                                                                                      cgcgGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGt                                                                 <  1:4921349/150‑1 (MQ=255)
                                                                                        cggcggGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTAt                                                               <  1:3533998/150‑1 (MQ=255)
                                                                                         ggcggGATTAAAACATTGCTAATTCCCTTCGAAAATAAACGAGTTCTGGAAGAGATTCCTTACAACGTAATTGCCGATATGGTCATTAATACTTTGAAGCGCATTGAGTAAGTTCTTTCTCTTGCGCTGCAAAATGAACCGTCTGGTATg                                                              <  2:5542495/150‑1 (MQ=255)
                                                                                            gggATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAg                                                           <  2:854849/150‑1 (MQ=255)
                                                                                              taTTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAATCGCATTGAGGAAGTTATGACTCTGGCGATGCAAAATGAACCGTCTGTTATGAAGGt                                                         >  1:2587648/2‑150 (MQ=255)
                                                                                                 tAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTtgt                                                      >  2:1743515/1‑150 (MQ=255)
                                                                                                  aaaaCAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTtgtg                                                     >  1:974811/1‑150 (MQ=255)
                                                                                                    aaCAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGAc                                                   >  1:2722427/1‑150 (MQ=255)
                                                                                                       aGTGCTAAGTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGc                                                <  1:3726478/150‑1 (MQ=255)
                                                                                                        gTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCa                                               >  2:573732/1‑150 (MQ=255)
                                                                                                          gCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCaaa                                             >  2:4296566/1‑150 (MQ=255)
                                                                                                            tAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAAt                                           <  1:2207253/150‑1 (MQ=255)
                                                                                                            tAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAAt                                           >  2:4938/1‑150 (MQ=255)
                                                                                                             aaTTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATa                                          <  2:1579683/150‑1 (MQ=255)
                                                                                                                tCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTg                                       <  1:908955/150‑1 (MQ=255)
                                                                                                                  cGTTCGAAAATAACCGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGAc                                     <  1:5481330/150‑1 (MQ=255)
                                                                                                                   gTTCGAAAATAAACGCGTTCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGAGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGAcc                                    <  2:4985845/150‑1 (MQ=255)
                                                                                                                    ttCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCt                                   <  1:2410447/150‑1 (MQ=255)
                                                                                                                     tCGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTc                                  <  2:3676038/150‑1 (MQ=255)
                                                                                                                      cGAAAATAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTcg                                 >  2:191666/1‑150 (MQ=255)
                                                                                                                       gAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTcgc                                <  2:255706/150‑1 (MQ=255)
                                                                                                                            tAAACGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCaaa                           >  1:2916665/1‑150 (MQ=255)
                                                                                                                             aaaCGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCaaaa                          <  1:4630682/150‑1 (MQ=255)
                                                                                                                              aaCGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAAt                         <  2:745549/150‑1 (MQ=255)
                                                                                                                              aaCGCGATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAAt                         >  1:1835534/1‑150 (MQ=255)
                                                                                                                               aCGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATg                        <  2:503505/150‑1 (MQ=255)
                                                                                                                                  cgATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCAc                     >  1:3740601/1‑150 (MQ=255)
                                                                                                                                  cgATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCAc                     >  1:4891158/1‑150 (MQ=255)
                                                                                                                                  cgATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCAc                     <  2:2350355/150‑1 (MQ=255)
                                                                                                                                  cgATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCAc                     >  1:2182189/1‑150 (MQ=255)
                                                                                                                                   gATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACt                    <  1:4700123/150‑1 (MQ=255)
                                                                                                                                   gATCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACt                    <  1:2823126/150‑1 (MQ=255)
                                                                                                                                    aTCTGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACta                   >  2:4319470/1‑150 (MQ=255)
                                                                                                                                     tCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCATAATGCACtaa                  <  1:4439749/150‑1 (MQ=255)
                                                                                                                                       tGGAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACtaata                <  1:1431825/150‑1 (MQ=255)
                                                                                                                                       tGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTACGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACtaata                <  2:4181902/150‑1 (MQ=255)
                                                                                                                                         gAAGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATaaa              >  2:199869/1‑150 (MQ=255)
                                                                                                                                          aaGAGATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATaaaa             >  2:2775900/1‑150 (MQ=255)
                                                                                                                                           agagATTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATaaaaa            >  2:3123635/1‑150 (MQ=255)
                                                                                                                                             agagTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACa          <  2:513845/150‑1 (MQ=255)
                                                                                                                                              gaTTCCTTACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAg         >  2:2140579/1‑150 (MQ=255)
                                                                                                                                                  ccTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCt     >  1:3063732/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGATGTGGCAATGACCGGTGAGATCACTCTGCGTGGTCAGGTACTGCCGATCGGTGGTTTGAAAGAAAAACTCCTGGCAGCGCATCGCGGCGGGATTAAAACAGTGCTAATTCCGTTCGAAAATAAACGCGATCTGGAAGAGATTCCTGACAACGTAATTGCCGATCTGGACATTCATCCTGTGAAGCGCATTGAGGAAGTTCTGACTCTGGCGCTGCAAAATGAACCGTCTGGTATGCAGGTTGTGACTGCAAAATAGTGACCTCGCGCAAAATGCACTAATAAAAACAGGGCTGGC  >  CP009273/456439‑456737

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: