Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 251,998 T→C 7.3% 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 7.3% 1014.4 / 51.7 315K69R (AAA→AGA) pepDaminoacyl‑histidine dipeptidase (peptidase D)
Reads supporting (aligned to +/- strand):  ref base T (131/160);  new base C (13/10);  total (144/171)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 318 total aligned reads displayed.

AATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAAC  >  CP009273/251850‑252145
                                                                                                                                                    |                                                                                                                                                   
 aTGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTa                                                                                                                                                   <  2:4094220/150‑1 (MQ=255)
     cGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCGGTAAGACGACCGGGTTACGAt                                                                                                                                               >  1:2797292/1‑150 (MQ=255)
       ttATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGAttt                                                                                                                                             >  1:3439078/1‑150 (MQ=255)
         aTCCGCACCCCGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTc                                                                                                                                           <  1:1859702/150‑1 (MQ=255)
            cGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCAt                                                                                                                                        <  2:3257311/150‑1 (MQ=255)
             gCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCCTCTCGAGGTGGGCCTGGACGACGACCGGGCTACGATTTTCCATa                                                                                                                                       >  2:4135475/1‑150 (MQ=255)
               aCCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATAcc                                                                                                                                     <  2:2434722/150‑1 (MQ=255)
                   aGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGAAGACCGGTCTACGATTTTCCATACCTGCg                                                                                                                                 >  2:867952/1‑150 (MQ=255)
                     cGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGt                                                                                                                               >  2:3055761/1‑150 (MQ=255)
                      gtggtgCCGCGCGCTTTAACCCATTCGCCATCAATATCAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTa                                                                                                                              <  1:281732/150‑1 (MQ=255)
                      gtggtgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTa                                                                                                                              <  1:411852/150‑1 (MQ=255)
                        ggtgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGc                                                                                                                            >  1:1022850/1‑150 (MQ=255)
                        ggtgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTACATAACAGCGGTAGc                                                                                                                            >  1:4262660/1‑150 (MQ=255)
                        ggtgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTATGCGGCACCATATACAGGTGGGCCAGAAAGACCACCGCTTTAAGAATTTCCATACCTGCGGTAGc                                                                                                                            >  1:3881393/1‑150 (MQ=255)
                          tgCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAg                                                                                                                          <  1:1802801/150‑1 (MQ=255)
                              gcgcgcTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGttt                                                                                                                      <  1:1317359/150‑1 (MQ=255)
                              gcgcgcTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTGTACGATTTTCCATACCTGCGGTAGCAGGttt                                                                                                                      <  2:2018785/150‑1 (MQ=255)
                               cgcgcTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTa                                                                                                                     <  2:2955356/150‑1 (MQ=255)
                                 cgcTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTGTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACg                                                                                                                   >  2:3094030/1‑150 (MQ=255)
                                    tttAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAAt                                                                                                                <  2:2236295/150‑1 (MQ=255)
                                       aaCCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAg                                                                                                             <  2:1813665/150‑1 (MQ=255)
                                       aaCCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAg                                                                                                             <  2:3420179/150‑1 (MQ=255)
                                        aCCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAgg                                                                                                            <  2:460376/150‑1 (MQ=255)
                                        aCCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAgg                                                                                                            <  1:2945760/150‑1 (MQ=255)
                                         cccATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTATGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATAACAGCGGTAGCAGGTTTAAGAATCAGGa                                                                                                           >  1:4190083/1‑150 (MQ=255)
                                           cATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGata                                                                                                         <  1:1137678/150‑1 (MQ=255)
                                            aTTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGTTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGatat                                                                                                        <  1:1234021/150‑1 (MQ=255)
                                            aTTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGatat                                                                                                        >  2:3604641/1‑150 (MQ=255)
                                              tCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATta                                                                                                      <  1:1879762/150‑1 (MQ=255)
                                                gCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGGAGCAGGTTTACGAATCAGGATATtacc                                                                                                    >  2:2698434/1‑150 (MQ=255)
                                                 ccATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtacct                                                                                                   <  2:767198/150‑1 (MQ=255)
                                                  cATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATAGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTACATAAATTCGGTAGCAGGTTTAAGAATAAGGATATtaccta                                                                                                  >  2:427817/1‑150 (MQ=255)
                                                    tCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTATGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTAAGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATtacctacc                                                                                                >  1:2376980/1‑150 (MQ=255)
                                                      aaTATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTg                                                                                              <  2:2773779/150‑1 (MQ=255)
                                                      aaTATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTATGCGGCACCATATCGAGGTGGGCCTGTAAGAAGACCGGTTTAAGATTTTACATAACAGCGGTAGCAGGTTTAAGAATCAGGATATTACCTAACTg                                                                                              >  1:1977898/1‑150 (MQ=255)
                                                      aaTATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTATGCGGCACCAAATCGAGGTGGGCCAGTAAGAAGACCGGTTTACGATTTTCCATACCAGCGGTAGCAGGTTTACGAATCAGGATATTACCAACCTg                                                                                              >  1:2840215/1‑150 (MQ=255)
                                                       atatAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTACATACCTGCGGTAGCAGGTTTTCGAATCATGATATTACCTAACTGa                                                                                             >  2:3504461/1‑150 (MQ=255)
                                                        tataAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGAt                                                                                            <  2:685724/150‑1 (MQ=255)
                                                        tataAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGAt                                                                                            <  2:3647059/150‑1 (MQ=255)
                                                         ataAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTATGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTACATACCAGCGGTAGCAGGTTTACCAATCAGGATATTAACAACCTGATc                                                                                           >  1:623397/1‑150 (MQ=255)
                                                          taAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCAAATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTAACTACCTGATcg                                                                                          >  1:3221357/1‑150 (MQ=255)
                                                             ggCTGGATAGGTTCTTTCGTGTAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCTTATCGCGGTGGGCCTGTTAGTCGAGCGGTTTACGATTTTGCATCCCTGCGGTAGCAGGTTTACGAATCTGGATATTACCTACCGGATCGCGt                                                                                       <  2:3125087/150‑1 (MQ=255)
                                                               cTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTc                                                                                     <  2:2396581/150‑1 (MQ=255)
                                                                tGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCg                                                                                    <  1:926947/150‑1 (MQ=255)
                                                                   aTAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTATGCGGCACCAAATCGAGGTGGGCCGGTAAGACGAACCGTTTAAGATTTTACAAAACAGCGGTATCAGGTTTTAGAAAAACGATATTAACAAACAGATAGCCTTCGACa                                                                                 >  1:2975610/1‑150 (MQ=255)
                                                                    tAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACAt                                                                                <  1:3707348/150‑1 (MQ=255)
                                                                    tAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATAGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTACAAAACTGCGGTATCAGGTTTACGAATCAGGATATTAACAACCAGATCGCGTTAGACAt                                                                                >  1:210824/1‑150 (MQ=255)
                                                                      ggATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGCGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATgg                                                                              <  2:2179035/150‑1 (MQ=255)
                                                                      ggATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATgg                                                                              <  2:1871869/150‑1 (MQ=255)
                                                                        aTCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGaa                                                                            <  1:2394231/150‑1 (MQ=255)
                                                                           tttCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAAcc                                                                         >  2:1936236/1‑150 (MQ=255)
                                                                           tttCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTAAGAAACAGGATATTACCTACCAGATCGCGTTCGACATGGAAAcc                                                                         >  1:2162689/1‑150 (MQ=255)
                                                                              cGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCttt                                                                      <  1:743366/150‑1 (MQ=255)
                                                                                tGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCAAATAGAGGTGGGCCGGTAAGACGACCGGGTTAAGATTTTCCATACCAGCGGTAGCAGGTTTACGAAAAAGGATATTAACAACCAGATCCCCTTCGACAAGGAAACCATtct                                                                    >  1:1573259/1‑150 (MQ=255)
                                                                                   aGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCAGATCGCGTTCGACATGGAAACCTTTATCtt                                                                 >  2:3158557/1‑150 (MQ=255)
                                                                                    gTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCttt                                                                <  2:653390/150‑1 (MQ=255)
                                                                                     tCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTACATACCTGCGGTAGCAGGTTTAAGAATAAGGATTTTAACTACCTGATCGCGTTCGACATGGAAAACATTCTAtttt                                                               >  1:970060/1‑150 (MQ=255)
                                                                                     tCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTCTACGATTTTCCATACCTGCGGTAGCAGGTTTAAGAATCAGGATATTAACAACCTGATCGCGTTCGACATGGAAACCTTTATAtttt                                                               >  1:4356809/1‑150 (MQ=255)
                                                                                       aTGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGGTCGACAAGGAAACCTTTCTCTTTTGc                                                             >  2:1541477/1‑150 (MQ=255)
                                                                                       aTGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTAACTACCTGATCGCGTTAGACAAGGAAACCTTTCTCGTTTGc                                                             >  2:1571197/1‑150 (MQ=255)
                                                                                        tGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGcc                                                            <  2:3440496/150‑1 (MQ=255)
                                                                                         gCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGccc                                                           <  2:593077/150‑1 (MQ=255)
                                                                                         gCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGccc                                                           <  2:2653495/150‑1 (MQ=255)
                                                                                          cACGGTGTCGTTATTTTTCGGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCcca                                                          <  2:4172828/150‑1 (MQ=255)
                                                                                              gtgtCGTTATTTTTATGCGGCACCAAATAGAGGTTGGCCTGTAAGACGACCCGTCTAAGATTTTACATACCAGCGGTAGCAGGTTTACGAAAAAGGATATTAACTAACAGATAGCGTTCGACAAGGAAACCTTTCTCTTTTTCccaacca                                                      >  1:3936399/1‑150 (MQ=255)
                                                                                                  cGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACcaacaa                                                  <  1:3055761/150‑1 (MQ=255)
                                                                                                   gTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAAt                                                 >  1:2001029/1‑150 (MQ=255)
                                                                                                   gTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAAt                                                 >  1:1178367/1‑150 (MQ=255)
                                                                                                   gTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAAt                                                 >  1:2325981/1‑150 (MQ=255)
                                                                                                     tATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGt                                               <  1:3086229/150‑1 (MQ=255)
                                                                                                     tATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGt                                               >  1:4021483/1‑150 (MQ=255)
                                                                                                        ttttCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGACACCTTTCTCTTTTGCCCAACCAACAATGTAtt                                            >  1:2534345/1‑150 (MQ=255)
                                                                                                        ttttCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTAtt                                            <  1:2027593/150‑1 (MQ=255)
                                                                                                          ttCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCa                                          <  2:2627510/150‑1 (MQ=255)
                                                                                                               cGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAg                                     <  2:2340307/150‑1 (MQ=255)
                                                                                                               cGGCACCATATCGAGGTGGGCCTGCAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAg                                     >  1:3605180/1‑150 (MQ=255)
                                                                                                                ggCGCCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGt                                    >  1:1740074/1‑150 (MQ=255)
                                                                                                                ggCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGt                                    >  2:3430179/1‑150 (MQ=255)
                                                                                                                 gCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGtt                                   <  2:37892/150‑1 (MQ=255)
                                                                                                                 gCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGtt                                   >  1:4238560/1‑150 (MQ=255)
                                                                                                                  cACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTg                                  >  1:2650487/1‑150 (MQ=255)
                                                                                                                       tatCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCACCCAACAATGTATTCAGCGAGTTGCTCtt                             >  2:1421493/1‑150 (MQ=255)
                                                                                                                       tatCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCtt                             >  2:2668020/1‑150 (MQ=255)
                                                                                                                          cGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCAt                          <  2:1772926/150‑1 (MQ=255)
                                                                                                                            aGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGa                        >  1:261313/1‑150 (MQ=255)
                                                                                                                              gTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATa                      >  1:4177368/1‑150 (MQ=255)
                                                                                                                               tGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAg                     <  1:3282834/150‑1 (MQ=255)
                                                                                                                                gggCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAgg                    <  2:1169734/150‑1 (MQ=255)
                                                                                                                                  gCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGAc                  <  1:3538403/150‑1 (MQ=255)
                                                                                                                                  gCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGAc                  <  1:162909/150‑1 (MQ=255)
                                                                                                                                   ccTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACg                 <  2:370050/150‑1 (MQ=255)
                                                                                                                                   ccTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACg                 <  2:3906302/150‑1 (MQ=255)
                                                                                                                                      gTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGt              >  2:1250393/1‑150 (MQ=255)
                                                                                                                                        aaGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGa            >  1:1469740/1‑150 (MQ=255)
                                                                                                                                         aGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAg           <  1:948999/150‑1 (MQ=255)
                                                                                                                                           acgacCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGa         <  2:1660887/150‑1 (MQ=255)
                                                                                                                                           acgacCGGTTTACGACTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGa         <  2:3534551/150‑1 (MQ=255)
                                                                                                                                             gacCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAAt       >  2:252959/1‑150 (MQ=255)
                                                                                                                                               ccGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAg     <  2:2307499/150‑1 (MQ=255)
                                                                                                                                                 ggTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGaa   >  2:4213273/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATGCCGTTATCCGCACCCAGCGTGGTGCCGCGCGCTTTAACCCATTCGCCATCAATATAAGGCTGGATAGGATCTTTCGTGAAGTCATGCACGGTGTCGTTATTTTTCTGCGGCACCATATCGAGGTGGGCCTGTAAGACGACCGGTTTACGATTTTCCATACCTGCGGTAGCAGGTTTACGAATCAGGATATTACCTACCTGATCGCGTTCGACATGGAAACCTTTCTCTTTTGCCCAACCAACAATGTATTCAGCGAGTTGCTCTTCATGATAGGACGGGTGAGGAATAGAAC  >  CP009273/251850‑252145

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