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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,922,8990GT13.1% 691.9 / 54.0 214I66I (ATC→ATA) ptrBprotease II
Reads supporting (aligned to +/- strand):  ref base G (92/94);  new base T (18/10);  total (110/104)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 237 total aligned reads displayed.

GAGCAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTA  >  CP009273/1922752‑1923046
                                                                                                                                                   |                                                                                                                                                   
gAGCAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCg                                                                                                                                                     <  1:495635/148‑1 (MQ=255)
gAGCAATGTTTACCACTCATCCCACTCTTCACTGAATGCCGATTGTCTCTGGTAGATAGCATATTCTCAGCCTGTTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCTGAAACTTCTCGTTGCGGAATGCGGTCg                                                                                                                                                     <  2:701471/148‑1 (MQ=255)
 aGCAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCga                                                                                                                                                    <  1:902233/148‑1 (MQ=255)
 aGCAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCga                                                                                                                                                    <  2:1041912/148‑1 (MQ=255)
   cAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGACTGCGGTCgatg                                                                                                                                                  >  2:1538662/1‑148 (MQ=255)
   cAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCgatg                                                                                                                                                  <  2:1708982/148‑1 (MQ=255)
   cAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCgatg                                                                                                                                                  >  1:1874332/1‑148 (MQ=255)
   cAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCgatg                                                                                                                                                  <  2:803277/148‑1 (MQ=255)
    aaTGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCgatga                                                                                                                                                 <  2:2808823/148‑1 (MQ=255)
     aTGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCgatgat                                                                                                                                                <  1:3081618/148‑1 (MQ=255)
     aTGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCgatgat                                                                                                                                                >  1:2053118/1‑148 (MQ=255)
      tGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGAtt                                                                                                                                               >  2:3316726/1‑148 (MQ=255)
         ttCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTcc                                                                                                                                            <  1:563118/148‑1 (MQ=255)
             cACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTa                                                                                                                                        >  1:2512098/1‑148 (MQ=255)
             cACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTa                                                                                                                                        <  2:1666164/148‑1 (MQ=255)
             cACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTa                                                                                                                                        <  2:348395/148‑1 (MQ=255)
             cACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTa                                                                                                                                        <  1:2984219/148‑1 (MQ=255)
              aCTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTaa                                                                                                                                       >  2:1961628/1‑148 (MQ=255)
              aCTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTaa                                                                                                                                       >  2:2365053/1‑148 (MQ=255)
               cTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAg                                                                                                                                      <  1:129777/148‑1 (MQ=255)
               cTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAg                                                                                                                                      <  1:3527820/148‑1 (MQ=255)
                 cATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGAt                                                                                                                                    >  2:713180/1‑148 (MQ=255)
                  aTCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATg                                                                                                                                   <  1:1787261/148‑1 (MQ=255)
                       aCTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTGCTTTAAGATTCGATc                                                                                                                              >  1:370467/1‑148 (MQ=255)
                       aCTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATc                                                                                                                              <  1:3409784/148‑1 (MQ=255)
                         tctTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGAAGCGATCCt                                                                                                                            >  1:1080930/1‑148 (MQ=255)
                          ctTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTACTTTAATATGCGATCCTg                                                                                                                           >  1:3578439/1‑148 (MQ=255)
                            tCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCa                                                                                                                         <  1:1342612/148‑1 (MQ=255)
                              aCTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAg                                                                                                                       >  2:1084252/1‑148 (MQ=255)
                              aCTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAg                                                                                                                       >  2:2991990/1‑148 (MQ=255)
                               cTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAgg                                                                                                                      >  1:1452520/1‑148 (MQ=255)
                               cTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAgg                                                                                                                      <  1:1228793/148‑1 (MQ=255)
                                  aaTGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATTCGATACTGCAAGGCtt                                                                                                                   >  1:1068300/1‑148 (MQ=255)
                                  aaTGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCtt                                                                                                                   <  1:3406343/148‑1 (MQ=255)
                                   aTGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCttg                                                                                                                  >  2:360329/1‑148 (MQ=255)
                                    tGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCttgt                                                                                                                 >  2:1403144/1‑148 (MQ=255)
                                    tGCCGATTGACGCTGGTAGATAGCAAATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTACTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATACTGCAAGGCttgt                                                                                                                 >  2:1841744/1‑148 (MQ=255)
                                      ccGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCttgttg                                                                                                               <  2:221806/148‑1 (MQ=255)
                                       cGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTtgttgt                                                                                                              <  2:2023883/148‑1 (MQ=255)
                                         aTTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTACTTTAAGATGCGATCCTGCAAGGCTTGTTGTGa                                                                                                            >  1:2095742/1‑148 (MQ=255)
                                           tGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTGTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAgg                                                                                                          >  2:3332184/1‑148 (MQ=255)
                                              cGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCATTAAAAATCCATACTGCCAGGCTTGTTGTGAGGCCa                                                                                                       >  1:507520/1‑148 (MQ=255)
                                                   gTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGTTTTCCTTTTAGTTTGGATCCAGCAAGGCGTGTTGTGAGGCCATCAcc                                                                                                  >  1:1339894/1‑148 (MQ=255)
                                                    tAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGAGATCCTGCAAGGCTTGTTGTGAGGCCATCAccc                                                                                                 >  2:2314686/1‑148 (MQ=255)
                                                      gATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGa                                                                                               >  2:2743780/1‑148 (MQ=255)
                                                      gATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGa                                                                                               <  1:2628081/148‑1 (MQ=255)
                                                             tatTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGt                                                                                        >  2:103967/1‑148 (MQ=255)
                                                              atTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTa                                                                                       >  2:3709969/1‑148 (MQ=255)
                                                                  acaGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTa                                                                                   <  1:7215/148‑1 (MQ=255)
                                                                   caGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTAt                                                                                  <  2:3620842/148‑1 (MQ=255)
                                                                       cTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGTTTTCCTTTAAGATTCGATCCAGCAAGGCTTGTTGTGAGGCCCACACCCCATGACCGTAACTAtttta                                                                              >  1:1749020/1‑147 (MQ=255)
                                                                         gggtCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCtt                                                                            <  2:1507076/145‑1 (MQ=255)
                                                                         ggTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCtt                                                                            <  2:1882807/148‑1 (MQ=255)
                                                                          gTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTGCTTTTAGATGCGATCCTGCCAGGCGTTGTGTGTGGCCATCACCCCATGACCGTAACAATTttgttg                                                                           >  1:3261442/1‑148 (MQ=255)
                                                                          gTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCTATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCttg                                                                           <  1:1180611/148‑1 (MQ=255)
                                                                           ttCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCttgt                                                                          <  2:3465089/148‑1 (MQ=255)
                                                                             cATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGAATTCCTTTAAGATGCGATCCATCAAGGCGTGTTGTGAGGCCAACACACGATGACCGTAACAATTTTCttgttt                                                                        >  2:3519745/1‑147 (MQ=255)
                                                                              aTAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGc                                                                       <  1:2179446/148‑1 (MQ=255)
                                                                               tAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGGTGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCa                                                                      <  2:1190520/148‑1 (MQ=255)
                                                                               tAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGAGGGTTTCCGTTAATATGCGTTACTGCCAGGGTTGTTGTGTGGCCAACACCCCATGACCGGCACTATTTTCTTGTTGCa                                                                      >  1:2253546/1‑148 (MQ=255)
                                                                                aaaTATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGTTTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAg                                                                     >  1:2865161/1‑148 (MQ=255)
                                                                                aaaTATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTACTTTAAGATTCGATCCAGCAAGGGTTGTTTTGAGGCCATCACCCCATGACCGTAAATATTTTATTGTTGGAg                                                                     >  1:836982/1‑148 (MQ=255)
                                                                                 aaTATGCCGATAGCGGTTGCCATTTTTGATGTTGGGCGCTGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAgg                                                                    <  2:2509355/148‑1 (MQ=255)
                                                                                 aaTATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAgg                                                                    <  2:535053/148‑1 (MQ=255)
                                                                                     tGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCAGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGt                                                                >  1:2829120/1‑148 (MQ=255)
                                                                                         gATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAg                                                            <  2:2341758/148‑1 (MQ=255)
                                                                                            aGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGAc                                                         <  1:2365053/148‑1 (MQ=255)
                                                                                              cGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACtt                                                       <  1:2003591/148‑1 (MQ=255)
                                                                                              cGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACtt                                                       <  2:2055244/148‑1 (MQ=255)
                                                                                              cGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACtt                                                       <  2:782775/148‑1 (MQ=255)
                                                                                                  aGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGAGGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCttg                                                   >  1:1238996/1‑146 (MQ=255)
                                                                                                     cATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTg                                                >  1:3258945/1‑148 (MQ=255)
                                                                                                     cATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTg                                                <  2:1264576/148‑1 (MQ=255)
                                                                                                     cATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTg                                                <  2:2484367/148‑1 (MQ=255)
                                                                                                      aTTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTACTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCCATGACCGTAACAATTTTATTGTTTCAGGTAGTCCAGGACTTCTTGCTga                                               >  1:3320315/1‑148 (MQ=255)
                                                                                                        ttttGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTgaga                                             >  1:3688204/1‑148 (MQ=255)
                                                                                                        ttttGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTgaga                                             <  1:1796547/148‑1 (MQ=255)
                                                                                                        ttttGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGAGGTTTTACATTAAGATGAGATCCAGCAAGGCTTGTTGTTTGGCCATCACCCCAAGACCCTAAATAATTTTTTTTTTGAGGTAGTCCAGGACTTATTGCTgaga                                             >  1:3198266/1‑148 (MQ=255)
                                                                                                         tttGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGaa                                            <  2:932773/148‑1 (MQ=255)
                                                                                                           tGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAAcg                                          >  2:622072/1‑148 (MQ=255)
                                                                                                              tGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGt                                       <  1:847363/148‑1 (MQ=255)
                                                                                                               gTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAATATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTa                                      >  2:752203/1‑148 (MQ=255)
                                                                                                                tAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATTCGATCCTGCAAGGCTTTTTTTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTAt                                     >  1:1558369/1‑148 (MQ=255)
                                                                                                                  gggCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCg                                   >  2:1689592/1‑148 (MQ=255)
                                                                                                                  gggCGCAGAAACTTCTCGTTGCGGAATGCGGTCGAGGATTTCCGTTAAGATTCGAACCTGCAAGGCGTGTTGTGAGGCCAACACCCGCTGACCGTAACTATTTTCGTGTTGCAGGTAGTCCAGGACTTCTGGCTGATAACCCGTaaag                                   >  2:93977/1‑145 (MQ=255)
                                                                                                                   ggCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCACGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGt                                  <  2:1132150/148‑1 (MQ=255)
                                                                                                                        aGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGc                             <  1:3628769/148‑1 (MQ=255)
                                                                                                                         gAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCAGCAAGGCTTGTTGTGTGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTATGGCTGAGAACGCGTATCGTCCCGCa                            >  2:2265698/1‑148 (MQ=255)
                                                                                                                          aaaCTTCTCGTTGCTGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAg                           <  1:234953/148‑1 (MQ=255)
                                                                                                                            aCTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGcc                         <  2:2711408/148‑1 (MQ=255)
                                                                                                                            aCTTCTCGTTGCGGAATGCGGTCTATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGcc                         <  2:794706/148‑1 (MQ=255)
                                                                                                                            aCTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCGTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGTGAACGCGTATCGTCCCGCAGcc                         >  2:2631979/1‑148 (MQ=255)
                                                                                                                              ttCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATTCGATCCTGCAAGGCTTGTTGGGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCag                       >  1:2170479/1‑148 (MQ=255)
                                                                                                                                     tGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAAt                <  2:1359626/148‑1 (MQ=255)
                                                                                                                                      gCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAAtt               <  1:626747/148‑1 (MQ=255)
                                                                                                                                       cGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTa              <  1:103967/148‑1 (MQ=255)
                                                                                                                                             gCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATg        <  1:208460/148‑1 (MQ=255)
                                                                                                                                                  cGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGt   >  2:2436028/1‑148 (MQ=255)
                                                                                                                                                   gatgatTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTa  <  2:2752291/148‑1 (MQ=255)
                                                                                                                                                   gatgatTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTa  <  1:985065/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAGCAATGTTTCCCACTCATCCCACTCTTCACTGAATGCCGATTGACGCTGGTAGATAGCATATTCACAGCCTGGTTCATAAATATGCCGATAGCGGTAGCCATTTTTGATGTAGGGCGCAGAAACTTCTCGTTGCGGAATGCGGTCGATGATTTCCTTTAAGATGCGATCCTGCAAGGCTTGTTGTGAGGCCATCACCCGATGACCGTAACTATTTTCTTGTTGCAGGTAGTCCAGGACTTCTGGCTGAGAACGCGTATCGTCCCGCAGCCAGTAGTAATTATCGATGCGCGTA  >  CP009273/1922752‑1923046

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: ‑