Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,922,899 G→T 5.6% I66I (ATC→ATA ptrB ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,922,8990GT 5.6% 960.3 / 19.5 250I66I (ATC→ATAptrBprotease II
Reads supporting (aligned to +/- strand):  ref base G (118/118);  new base T (8/6);  total (126/124)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 264 total aligned reads displayed.

GAGCAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGT  >  CP009273/1922752‑1923045
                                                                                                                                                   |                                                                                                                                                  
 aGCAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCTTTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCga                                                                                                                                                   <  1:1206180/148‑1 (MQ=255)
  gCAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCgat                                                                                                                                                  >  1:2187424/1‑148 (MQ=255)
   cAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCgatg                                                                                                                                                 >  2:3384410/1‑148 (MQ=255)
     aTGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCgatgat                                                                                                                                               >  1:1459792/1‑148 (MQ=255)
       gTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGAttt                                                                                                                                             <  1:2862186/148‑1 (MQ=255)
            ccACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCttt                                                                                                                                        >  2:2294964/1‑148 (MQ=255)
              aCTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTaa                                                                                                                                      >  1:898939/1‑148 (MQ=255)
               cTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAg                                                                                                                                     <  1:2091471/148‑1 (MQ=255)
                 cATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGAt                                                                                                                                   >  1:2490574/1‑148 (MQ=255)
                  aTCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATg                                                                                                                                  <  1:1351969/148‑1 (MQ=255)
                  aTCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGAGGATTTCCTTTAAGATg                                                                                                                                  >  1:459536/1‑148 (MQ=255)
                   tCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGc                                                                                                                                 >  2:1393602/1‑148 (MQ=255)
                      cACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGAGGATTTGCTTTAAGATGCGAt                                                                                                                              >  1:922667/1‑148 (MQ=255)
                           ttCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGAGGATTTACTTTAAGATGCGATCCTGc                                                                                                                         >  1:3264050/1‑148 (MQ=255)
                            tCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCa                                                                                                                        <  1:2765730/148‑1 (MQ=255)
                              aCTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAg                                                                                                                      <  2:3337554/148‑1 (MQ=255)
                               cTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATGCTGCAAgg                                                                                                                     >  2:2683083/1‑148 (MQ=255)
                               cTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAgg                                                                                                                     >  2:775041/1‑148 (MQ=255)
                                 gAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGGTTTACTTTAACATGCGATCCAGCAACGCt                                                                                                                   >  1:1542396/1‑148 (MQ=255)
                                   aTGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCttg                                                                                                                 >  2:2167561/1‑148 (MQ=255)
                                     gCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCttgtt                                                                                                               <  1:3153438/148‑1 (MQ=255)
                                        gATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTtgtg                                                                                                            <  1:60250/148‑1 (MQ=255)
                                        gATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTtgtg                                                                                                            <  1:2402834/148‑1 (MQ=255)
                                         aTTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGa                                                                                                           <  1:1044006/148‑1 (MQ=255)
                                                cTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCCGTAGCCATATTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAACGCTTGTTGGGAGGCCCTc                                                                                                    >  1:2573283/1‑148 (MQ=255)
                                                  ggTAGATAGCATATTCAGAGGCTGGTTCATAAATATGGCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCAc                                                                                                  <  2:3239936/148‑1 (MQ=255)
                                                  ggTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCGTTAATATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCAc                                                                                                  >  1:863780/1‑148 (MQ=255)
                                                     aGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCg                                                                                               >  2:2062739/1‑148 (MQ=255)
                                                      gATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAGGGCTTGTTGTGAGGCCATCACCCGa                                                                                              <  1:3138762/148‑1 (MQ=255)
                                                       aTAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGAt                                                                                             <  2:207904/148‑1 (MQ=255)
                                                       aTAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGAt                                                                                             <  2:1355998/148‑1 (MQ=255)
                                                         aGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGCTGa                                                                                           >  1:2209469/1‑148 (MQ=255)
                                                             tatTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGt                                                                                       <  2:2817261/148‑1 (MQ=255)
                                                             tatTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGAGGGTTTCCGTTAAGATGCGATCCTGCCAGGCTTGTTGTGAGGCCCTAACCCCAAGACCgg                                                                                       >  1:3410048/1‑147 (MQ=255)
                                                              atTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTa                                                                                      >  2:2418092/1‑148 (MQ=255)
                                                               ttCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTaa                                                                                     >  2:390612/1‑148 (MQ=255)
                                                               gtCACAGCCTGGTTCATAACTATGCCGCTAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTaa                                                                                     <  1:2713620/147‑1 (MQ=255)
                                                                  acaGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGCTGACCGTAACTa                                                                                  >  2:168334/1‑148 (MQ=255)
                                                                  acaGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTa                                                                                  >  2:2511774/1‑148 (MQ=255)
                                                                  acaGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCATTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTa                                                                                  >  1:2649496/1‑148 (MQ=255)
                                                                   caGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATNCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTAt                                                                                 <  1:1801204/148‑1 (MQ=255)
                                                                     gCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTAttt                                                                               >  2:3239240/1‑148 (MQ=255)
                                                                        tGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCt                                                                            >  2:2570470/1‑148 (MQ=255)
                                                                        tGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCt                                                                            >  2:433475/1‑148 (MQ=255)
                                                                         ggTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCtt                                                                           <  1:2556763/148‑1 (MQ=255)
                                                                          gTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCttg                                                                          <  1:102940/148‑1 (MQ=255)
                                                                               tAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATTCGATCCTGCAAGGCGTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCa                                                                     >  2:1270170/1‑148 (MQ=255)
                                                                                aaaTATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGTTTTCCTTTAAGATGCGATCCTGCCAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAg                                                                    >  2:2584006/1‑148 (MQ=255)
                                                                                aaaTATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAg                                                                    <  1:819526/148‑1 (MQ=255)
                                                                                  atatGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGt                                                                  <  2:649526/148‑1 (MQ=255)
                                                                                  atatGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGt                                                                  <  2:2244839/148‑1 (MQ=255)
                                                                                  atatGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGt                                                                  <  1:166916/148‑1 (MQ=255)
                                                                                  atatGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAAATATTTTCTTGTTGCAGGt                                                                  >  2:2852203/1‑148 (MQ=255)
                                                                                    atGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAg                                                                >  1:486417/1‑148 (MQ=255)
                                                                                     tGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGt                                                               <  1:2135731/148‑1 (MQ=255)
                                                                                      gCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTc                                                              <  2:1143233/148‑1 (MQ=255)
                                                                                      gCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTc                                                              <  1:1393602/148‑1 (MQ=255)
                                                                                       ccGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTcc                                                             <  2:3060670/148‑1 (MQ=255)
                                                                                          aTAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAgg                                                          >  2:1923375/1‑148 (MQ=255)
                                                                                           tAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGAGGATTTGCTTTAAGATTCGATCCGGCCAGGCTTGTTGTGAGGCCCAAACCCCATGACCGTAAAAATTTTCTTGTTTGAGGTAGTCCAGGa                                                         >  1:3064668/1‑148 (MQ=255)
                                                                                           tAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGAGGATTTCCTTTAAGATGCGCTCCTGCAAGGCTTGTTGTGAGGCCAGCACCCCCTGACCGTAACAATTTTATTGTTGCAGGTAGTCCAGGa                                                         >  1:73467/1‑148 (MQ=255)
                                                                                           tAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGAGGATTTACTTTAAGATGCGATCCTGCAAGGCGTGTTGTGAGGCCATCACCCCATGACCGTAACAATTTTATTGTTGGAGGTAGTCCAGGa                                                         >  1:1873091/1‑148 (MQ=255)
                                                                                            aGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACGTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGAc                                                        <  2:1852681/148‑1 (MQ=255)
                                                                                             gCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACt                                                       >  2:2874734/1‑148 (MQ=255)
                                                                                                gTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCt                                                    <  1:262208/148‑1 (MQ=255)
                                                                                                 tAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTg                                                   <  2:2899473/148‑1 (MQ=255)
                                                                                                   gCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGc                                                 >  2:2990710/1‑148 (MQ=255)
                                                                                                   gCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGc                                                 >  2:1209926/1‑148 (MQ=255)
                                                                                                     cATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTg                                               <  1:3242559/148‑1 (MQ=255)
                                                                                                     cATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTg                                               <  1:3360425/148‑1 (MQ=255)
                                                                                                      aTTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTga                                              >  2:143200/1‑148 (MQ=255)
                                                                                                        ttttGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTgaga                                            <  2:897832/148‑1 (MQ=255)
                                                                                                         tttGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGaa                                           <  2:2406194/148‑1 (MQ=255)
                                                                                                          ttGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAAc                                          <  2:2431049/148‑1 (MQ=255)
                                                                                                              tGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGt                                      <  2:2649496/148‑1 (MQ=255)
                                                                                                                tAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTAt                                    <  2:781450/148‑1 (MQ=255)
                                                                                                                  gggCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCg                                  <  1:1505423/148‑1 (MQ=255)
                                                                                                                  gggCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCg                                  >  2:561592/1‑148 (MQ=255)
                                                                                                                   ggCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGt                                 <  1:1377645/148‑1 (MQ=255)
                                                                                                                   ggCGCAGAAACTTATCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGt                                 <  1:191380/148‑1 (MQ=255)
                                                                                                                    gcgcAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTc                                <  1:1708099/148‑1 (MQ=255)
                                                                                                                    gcgcAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTc                                <  1:2413140/148‑1 (MQ=255)
                                                                                                                    gcgcAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTc                                <  1:228433/148‑1 (MQ=255)
                                                                                                                    gcgcAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTc                                <  1:2191137/148‑1 (MQ=255)
                                                                                                                        aGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGc                            <  2:1311401/148‑1 (MQ=255)
                                                                                                                        aGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCGTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCAACACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGc                            >  2:399974/1‑148 (MQ=255)
                                                                                                                          aaaCTTCTCGTTGCGGAATGCGGTCGATGATTTACTTTAAGATGCGATCCAGCAAGGCTTTTTGTGAGGCCAACACCCGATGACCGTAACAATTTTCTTGTTGCAGGTAGTCCAGGACTTCTTGCTGAGAACGCGTATCGTCCCGCAg                          >  1:3523365/1‑148 (MQ=255)
                                                                                                                              ttCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCag                      <  2:1542396/148‑1 (MQ=255)
                                                                                                                               tctcGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCagt                     <  2:1757744/148‑1 (MQ=255)
                                                                                                                               tctcGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCagt                     >  2:1855479/1‑148 (MQ=255)
                                                                                                                               tctcGTTGCGGAATGCGGTCGATGATTTCCATTAAGATGCGATCCTGCAAGGCGTGTTGTGAGGCCAACACCCCATGACCGGAACTATTTTCGTGTTGCAGGGAGTCCAGGACGTCTGGCGGAGAACGCGTATCGTCCCGCAGCCcgt                     >  2:330820/1‑148 (MQ=255)
                                                                                                                                  cGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCagtagt                  <  2:2195285/148‑1 (MQ=255)
                                                                                                                                   gTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAgtagta                 <  1:3342623/148‑1 (MQ=255)
                                                                                                                                    ttGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTaa                <  2:1460539/148‑1 (MQ=255)
                                                                                                                                      gCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAAtt              <  1:1297005/148‑1 (MQ=255)
                                                                                                                                       cGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTa             <  1:1732863/148‑1 (MQ=255)
                                                                                                                                         gAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATc           >  1:2672248/1‑148 (MQ=255)
                                                                                                                                         gAATGCGGTCGATGATTTACTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATAACCCGATGACCGTAACAAATTTATTTTTGCAGGTAGTCCAGGACTTCTGGCTGAGAAAGCGTATCGTACCGCAGCCAGGAGTAATTTTc           >  2:1048953/1‑148 (MQ=255)
                                                                                                                                                gTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATgcgc    >  2:1602876/1‑148 (MQ=255)
                                                                                                                                                  cGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGt  >  2:468192/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAGCAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGT  >  CP009273/1922752‑1923045

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