Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 251,998 T→C 8.1% K69R (AAA→AGA)  pepD ← aminoacyl‑histidine dipeptidase (peptidase D)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273251,9980TC 8.1% 749.5 / 41.4 233K69R (AAA→AGA) pepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base T (81/132);  new base C (10/9);  total (91/142)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 236 total aligned reads displayed.

CAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACA  >  CP009273/251849‑252146
                                                                                                                                                     |                                                                                                                                                    
cAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGtt                                                                                                                                                      >  2:65300/1‑150 (MQ=255)
     ccGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGa                                                                                                                                                 <  2:270110/150‑1 (MQ=255)
       gTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGAtt                                                                                                                                               <  1:1391250/150‑1 (MQ=255)
            ccGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCa                                                                                                                                          >  2:1731617/1‑150 (MQ=255)
             cGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCAt                                                                                                                                         <  2:1664932/150‑1 (MQ=255)
                 cccAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCATACCt                                                                                                                                     <  1:3293488/150‑1 (MQ=255)
                  ccAGCGTGGTGCCGCGCGCTTTAACCCATACGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTg                                                                                                                                    >  1:2477009/1‑150 (MQ=255)
                    aGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGc                                                                                                                                   >  1:1602934/1‑149 (MQ=255)
                       gtggtgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTa                                                                                                                               >  2:143813/1‑150 (MQ=255)
                        tggtgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCTTATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAg                                                                                                                              <  1:3467403/150‑1 (MQ=255)
                        tggtgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAg                                                                                                                              <  1:932934/150‑1 (MQ=255)
                        tggtgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAg                                                                                                                              <  2:881275/150‑1 (MQ=255)
                        tggtgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAg                                                                                                                              <  1:1194576/150‑1 (MQ=255)
                          gtgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCa                                                                                                                            <  2:1604711/150‑1 (MQ=255)
                          gtgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCa                                                                                                                            <  2:2178679/150‑1 (MQ=255)
                          gtgCCGCGCGCTTTAACCCAGTCGACATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGATTACGATTTTCCATACCTGCGGTAGCa                                                                                                                            <  2:2790160/150‑1 (MQ=255)
                           tgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAg                                                                                                                           <  1:1542056/150‑1 (MQ=255)
                            gCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAgg                                                                                                                          >  2:2983490/1‑150 (MQ=255)
                             ccGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGt                                                                                                                         <  1:3442761/150‑1 (MQ=255)
                              cgcgcgCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGtt                                                                                                                        >  2:3902505/1‑150 (MQ=255)
                               gcgcgcTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGttt                                                                                                                       <  2:2180699/150‑1 (MQ=255)
                                 gcgcTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTAc                                                                                                                     <  2:3744181/150‑1 (MQ=255)
                                   gcTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGa                                                                                                                   >  2:2077399/1‑150 (MQ=255)
                                    cTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGaa                                                                                                                  >  2:3298555/1‑150 (MQ=255)
                                      ttAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATc                                                                                                                <  1:1163147/150‑1 (MQ=255)
                                      ttAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATc                                                                                                                <  1:3862796/150‑1 (MQ=255)
                                      ttAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATc                                                                                                                <  1:3823286/150‑1 (MQ=255)
                                       tAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCa                                                                                                               >  2:751957/1‑150 (MQ=255)
                                        aaCCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAg                                                                                                              >  1:84330/1‑150 (MQ=255)
                                           ccATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGat                                                                                                           >  2:3434421/1‑150 (MQ=255)
                                            cATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGata                                                                                                          >  1:2237807/1‑150 (MQ=255)
                                             aTTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGatat                                                                                                         >  2:3363337/1‑150 (MQ=255)
                                              gtCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATAtt                                                                                                        <  2:4261659/149‑1 (MQ=255)
                                                  ccATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtacct                                                                                                    <  1:2666929/150‑1 (MQ=255)
                                                  ccATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtacct                                                                                                    <  1:1070254/150‑1 (MQ=255)
                                                   cATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtaccta                                                                                                   >  2:2342326/1‑150 (MQ=255)
                                                   cATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtaccta                                                                                                   >  2:37563/1‑150 (MQ=255)
                                                     tCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtacctacc                                                                                                 <  1:2292883/150‑1 (MQ=255)
                                                       aaTATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTg                                                                                               <  2:2643777/150‑1 (MQ=255)
                                                       aaTATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTg                                                                                               >  1:1148574/1‑150 (MQ=255)
                                                          ataAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATc                                                                                            >  2:3908594/1‑150 (MQ=255)
                                                          ataAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATc                                                                                            >  1:2610812/1‑150 (MQ=255)
                                                          ataAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTGTTCTGCGGCACCATCTCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATc                                                                                            <  2:1380397/150‑1 (MQ=255)
                                                           taAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATcg                                                                                           <  2:4064764/150‑1 (MQ=255)
                                                           taAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATcg                                                                                           <  2:1606049/150‑1 (MQ=255)
                                                               gCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGtt                                                                                       <  1:1835761/150‑1 (MQ=255)
                                                                cTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTc                                                                                      <  1:752807/150‑1 (MQ=255)
                                                                cTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTc                                                                                      <  1:909839/150‑1 (MQ=255)
                                                                    aTAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACa                                                                                  <  1:4088555/150‑1 (MQ=255)
                                                                        gATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGa                                                                              <  1:652898/150‑1 (MQ=255)
                                                                        gATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGa                                                                              >  1:4136439/1‑150 (MQ=255)
                                                                         aTCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGaa                                                                             <  2:2610812/150‑1 (MQ=255)
                                                                          tCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGaaa                                                                            >  1:4029524/1‑150 (MQ=255)
                                                                              tCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCtt                                                                        <  2:2995417/150‑1 (MQ=255)
                                                                              tCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCtt                                                                        >  2:3307289/1‑150 (MQ=255)
                                                                               cGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCttt                                                                       >  1:658727/1‑150 (MQ=255)
                                                                                gTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTtc                                                                      <  2:3391442/150‑1 (MQ=255)
                                                                                  gAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTtctc                                                                    >  2:2956582/1‑150 (MQ=255)
                                                                                     gTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCttt                                                                 <  2:2447410/150‑1 (MQ=255)
                                                                                       cATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTg                                                               >  1:4308113/1‑150 (MQ=255)
                                                                                       cATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTg                                                               >  1:3247427/1‑150 (MQ=255)
                                                                                       cATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTg                                                               >  1:3969853/1‑150 (MQ=255)
                                                                                       cATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTg                                                               >  1:2063138/1‑150 (MQ=255)
                                                                                          gCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGccc                                                            <  1:2491890/150‑1 (MQ=255)
                                                                                          gCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGccc                                                            <  1:4527676/150‑1 (MQ=255)
                                                                                           cACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCcca                                                           >  2:3637733/1‑150 (MQ=255)
                                                                                              ggTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGTCGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCccaacc                                                        <  1:3060745/150‑1 (MQ=255)
                                                                                              ggTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCccaacc                                                        <  1:3060723/150‑1 (MQ=255)
                                                                                              ggTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCccaacc                                                        <  2:2873273/150‑1 (MQ=255)
                                                                                               gtgtCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCccaacca                                                       >  2:2477144/1‑150 (MQ=255)
                                                                                                 gtCGTTATTTGTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCcaaccaac                                                     <  2:2720834/150‑1 (MQ=255)
                                                                                                   cGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACcaacaa                                                   >  1:3451092/1‑150 (MQ=255)
                                                                                                   cGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACcaacaa                                                   >  2:1108133/1‑150 (MQ=255)
                                                                                                    gTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAAt                                                  <  1:2733847/150‑1 (MQ=255)
                                                                                                    gTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAAt                                                  >  1:1219479/1‑150 (MQ=255)
                                                                                                     ttATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATg                                                 >  1:599866/1‑150 (MQ=255)
                                                                                                      tATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGt                                                >  2:2452765/1‑150 (MQ=255)
                                                                                                      tATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGt                                                <  1:3849280/150‑1 (MQ=255)
                                                                                                       aTTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTa                                               >  1:72778/1‑150 (MQ=255)
                                                                                                           ttatGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCa                                           <  1:3893490/147‑1 (MQ=255)
                                                                                                           ttCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCa                                           <  2:244783/150‑1 (MQ=255)
                                                                                                              tGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCg                                        >  2:2424823/1‑150 (MQ=255)
                                                                                                                  gCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGtt                                    <  2:1006427/150‑1 (MQ=255)
                                                                                                                           cGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCAt                           >  2:187469/1‑150 (MQ=255)
                                                                                                                            gAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATg                          >  2:4299162/1‑150 (MQ=255)
                                                                                                                             aGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGa                         <  1:751957/150‑1 (MQ=255)
                                                                                                                              ggTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGAt                        <  1:1600599/150‑1 (MQ=255)
                                                                                                                               gTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATa                       <  2:1148574/150‑1 (MQ=255)
                                                                                                                                   gCCTGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGAc                   <  2:2304198/150‑1 (MQ=255)
                                                                                                                                     cTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACgg                 <  2:4010903/150‑1 (MQ=255)
                                                                                                                                     cTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACgg                 <  2:156114/150‑1 (MQ=255)
                                                                                                                                      tGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACggg                <  2:3231131/150‑1 (MQ=255)
                                                                                                                                         aaGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGa             >  2:569009/1‑150 (MQ=255)
                                                                                                                                           gacgacCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAgg           <  1:2956582/150‑1 (MQ=255)
                                                                                                                                             cgacCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGaa         >  2:4543098/1‑150 (MQ=255)
                                                                                                                                             cgacCGGTTTACGATTGTCCATACCTGCGGTAGCAGGTTTACGCATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTCTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGaa         <  1:1734540/150‑1 (MQ=255)
                                                                                                                                                ccGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAg      <  2:4273970/150‑1 (MQ=255)
                                                                                                                                                  ggTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGaa    <  2:2863866/150‑1 (MQ=255)
                                                                                                                                                   gTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAAc   <  2:658727/150‑1 (MQ=255)
                                                                                                                                                    tttACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACa  >  1:2417614/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAACA  >  CP009273/251849‑252146

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