Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,053,907 G→T 6.3% L312M (CTG→ATG)  pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,053,9070GT 6.3% 262.4 / 8.0 95L312M (CTG→ATG) pepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base G (50/39);  new base T (3/3);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.03e-01

GACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTCACAACCCG  >  NC_000913/3053758‑3054056
                                                                                                                                                     |                                                                                                                                                     
gACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAg                                                                                                                                                       <  2:590109/150‑1 (MQ=255)
  cATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGaa                                                                                                                                                     >  1:1306790/1‑150 (MQ=255)
   aTCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAAt                                                                                                                                                    <  1:737043/150‑1 (MQ=255)
    tacAGTCCTACCCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCTCTAACCATGATGCGCACCACATCACCAGTGACTTCCAGAATg                                                                                                                                                   >  2:1218005/3‑150 (MQ=255)
    tCCAGTCCTAACCAGTGGCTACCGCCATGCATACAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATg                                                                                                                                                   <  2:102835/150‑1 (MQ=255)
    tCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATg                                                                                                                                                   <  2:515267/150‑1 (MQ=255)
      cAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGa                                                                                                                                                 >  2:822474/1‑150 (MQ=255)
      cAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGa                                                                                                                                                 >  2:492963/1‑150 (MQ=255)
        gTCCTAACCAGTGGCTAAGGCCATGCATATTGAAAGGACGATGGGCGTTCTGAGCGTTCAGTTCATCTACATCACCTTTCAGGATGCCGAGTTTTACCATGCCGCTAAGCATGATGCGCACCACCTCACCTGTGACTTCCAGAATGGAAg                                                                                                                                               <  2:926368/150‑1 (MQ=255)
          ccTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGtt                                                                                                                                             >  2:792974/1‑150 (MQ=255)
           cTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGa                                                                                                                                                 >  1:1100909/1‑145 (MQ=255)
           cTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGa                                                                                                                                                 <  2:1100909/145‑1 (MQ=255)
           cTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGGTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCATTCTGGAAGtca                                                                                                                                            <  1:458964/150‑3 (MQ=255)
              aCCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCATTCTGGAAGtcactg                                                                                                                                         <  2:1450131/150‑6 (MQ=255)
               ccAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCATTCTGGAAGtcactgg                                                                                                                                        <  1:55820/150‑7 (MQ=255)
                    ggCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGAt                                                                                                                                   >  1:317685/1‑150 (MQ=255)
                       tAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTGATCAGCATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGCTGCGCACCACTGCACCAGTGACTTCCAGAATTGAAGTTCACGGAGGATAcc                                                                                                                                >  2:392525/1‑149 (MQ=255)
                        aaGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAg                                                                                                                               >  2:597056/1‑150 (MQ=255)
                          ggCCATGCATTAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGAAGATACAGgc                                                                                                                             >  2:672274/1‑150 (MQ=255)
                          ggCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGgc                                                                                                                             <  1:185135/150‑1 (MQ=255)
                           gCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGgct                                                                                                                            >  1:1360260/1‑149 (MQ=255)
                           gCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGgcg                                                                                                                            <  1:855261/150‑1 (MQ=255)
                           gCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGgcg                                                                                                                            <  1:898418/150‑1 (MQ=255)
                            cgaTGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCCACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGgcgc                                                                                                                           <  1:1604221/148‑1 (MQ=255)
                            ccATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGgcgc                                                                                                                           <  1:511246/150‑1 (MQ=255)
                            ccATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCATTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGgcgc                                                                                                                           >  2:619370/1‑150 (MQ=255)
                             cATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCa                                                                                                                          <  1:792974/150‑1 (MQ=255)
                             cATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAACTTCACGGACGATACAGGCGCa                                                                                                                          >  1:1122133/1‑150 (MQ=255)
                              aTGCATAAAGAAACGACGATGGGCTTTCTGAGCGATCAGTTCAACAACATCACCTTTCAGGATTCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAg                                                                                                                         <  2:916682/150‑1 (MQ=255)
                                  aTACAGAACGGTCGATGGGCGGTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGACGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTAACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTg                                                                                                                     <  2:107314/150‑1 (MQ=255)
                                   tAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTgg                                                                                                                    <  1:50982/150‑1 (MQ=255)
                                    aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCGGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGt                                                                                                                   >  1:1142658/1‑150 (MQ=255)
                                    aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGt                                                                                                                   >  2:1588195/1‑150 (MQ=255)
                                    aaagaaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGACGCTAACCATGATGCGCACATCTTCAACAGTGACTTCCTGAATGGAAGTTCAAGGACGATACAGGCGCAGGCTGGt                                                                                                                   >  2:473140/1‑150 (MQ=255)
                                        aaagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCg                                                                                                               <  2:983481/150‑1 (MQ=255)
                                         aagGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCga                                                                                                              >  2:548950/1‑150 (MQ=255)
                                          agGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCgag                                                                                                             <  2:1654630/150‑1 (MQ=255)
                                          agGACGATGGGCGTTCTGAGCGAGCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCgag                                                                                                             <  1:1153571/150‑1 (MQ=255)
                                           ggACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAAGCATGATGCGCACCACTTCACCAGAGACTTCCAGAATGGAAGTTCCCTGACTCGACAGGCGCAGGCTGGTTTCgaga                                                                                                            >  1:68919/1‑150 (MQ=255)
                                           ggACGATGGGCGTTCTGAGCGATCAGTACATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCGGACGATACAGGCGCAGGCTGGTTTCgagag                                                                                                           >  1:1636088/1‑150 (MQ=255)
                                             aCGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGAAGGCTGGTTTCgagagc                                                                                                          >  1:1279368/1‑149 (MQ=255)
                                              cGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCAATTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCTCAGGCTGGTTTCGAGAGAc                                                                                                         <  1:920425/150‑1 (MQ=255)
                                              cGATGGGCGTTCTGACCTATCAGTTTATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGATAACCATGATGCGCACCACTTAACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGAc                                                                                                         <  1:142243/150‑1 (MQ=255)
                                                        tCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTGTCGAGAGCCTCCAGCACaa                                                                                               >  2:1304446/1‑150 (MQ=255)
                                                         cTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAAt                                                                                              >  2:967213/1‑150 (MQ=255)
                                                           gAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTCACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGCCGATACAGGCGCAGGAGGGTTTCGAGAGACTCCAGCACAATGt                                                                                            >  2:346259/1‑150 (MQ=255)
                                                           gAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGCCTCCAGCACAATGt                                                                                            >  2:1398834/1‑150 (MQ=255)
                                                             gCGATCAGTTCATCAACATCACCGTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCg                                                                                          >  2:1083114/1‑150 (MQ=255)
                                                               gATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGtt                                                                                        >  1:1252456/1‑149 (MQ=255)
                                                                aTCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAg                                                                                       <  1:539541/150‑1 (MQ=255)
                                                                   aGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGAtt                                                                                    >  1:1633243/1‑150 (MQ=255)
                                                                   aGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGAtt                                                                                    >  2:1416605/1‑150 (MQ=255)
                                                                      tcatcaACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTcc                                                                                 >  2:606114/1‑149 (MQ=255)
                                                                      tcatcaACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCa                                                                                 >  2:843314/1‑150 (MQ=255)
                                                                      tcatcaACATCACCTTTCAGGATGCAGAGTTTTACCAGGCCGCTAACCATGATCCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGTACGATACAGGAGCAGGCTGGTTTCGAGAGACTCCACCACAAAGTCGAAGATTTCa                                                                                 >  1:296843/1‑150 (MQ=255)
                                                                         tcaACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGc                                                                              >  2:1606880/1‑150 (MQ=255)
                                                                             cATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTggg                                                                          >  1:1131717/1‑150 (MQ=255)
                                                                               tCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGcc                                                                        <  1:548950/150‑1 (MQ=255)
                                                                                 aCCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTg                                                                      <  1:597056/150‑1 (MQ=255)
                                                                                   cTTTCATGATGCCGAGTTTTACCCGGCAGCTAACCATGATGCTCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGATATACTACATCACAATGTCGTAGATTTCACGCTGGGCCTggg                                                                    <  2:1078695/150‑1 (MQ=255)
                                                                                     ttCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTg                                                                  <  2:1306790/150‑1 (MQ=255)
                                                                                      tCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCCGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGa                                                                 >  2:1126325/1‑150 (MQ=255)
                                                                                       cAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCCCTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGaa                                                                <  1:1464233/150‑1 (MQ=255)
                                                                                        aGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGCAt                                                               >  1:1395847/1‑150 (MQ=255)
                                                                                         ggATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCCGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAAtt                                                              >  2:1315713/1‑150 (MQ=255)
                                                                                             gCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGACACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGcc                                                          >  2:274564/1‑150 (MQ=255)
                                                                                              ccGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCg                                                         >  1:576816/1‑150 (MQ=255)
                                                                                               cGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTCGGCCTGGGTGAATTTGCCGt                                                        >  1:1143362/1‑150 (MQ=255)
                                                                                                gAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGtt                                                       >  2:750381/1‑150 (MQ=255)
                                                                                                      tACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCgg                                                 >  1:1705793/1‑150 (MQ=255)
                                                                                                                 tAACCATGATGCGCACCACTTCACCAGTGACTTCCAGACTGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTAGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCggg                                      <  2:545821/150‑1 (MQ=255)
                                                                                                                   aCCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCGGGTTTCGACAGACTCCAGCACAATGTCGTAGATTTCAAGCTGGGCCTGGGTGAATTTGCCGTTGACCTGGAAGGTGCGGGTa                                    >  2:330104/1‑150 (MQ=255)
                                                                                                                    ccATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCTAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTaa                                   >  1:156672/1‑150 (MQ=255)
                                                                                                                              caccacTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGc                         <  2:1108774/150‑1 (MQ=255)
                                                                                                                                ccacTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGACTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGt                       >  1:1571441/1‑150 (MQ=255)
                                                                                                                                    ttCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAAcc                   <  2:985022/150‑1 (MQ=255)
                                                                                                                                    ttCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAAcc                   <  2:1142658/150‑1 (MQ=255)
                                                                                                                                      cACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCtt                 >  2:789637/1‑150 (MQ=255)
                                                                                                                                       aCCAGTGACTTCCAGAATGGTAGTTCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTg               <  2:1636088/150‑1 (MQ=255)
                                                                                                                                        ccgggaaCTTCCATTCTGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTg               >  1:1445991/7‑150 (MQ=255)
                                                                                                                                        ccAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTg               >  2:767802/1‑150 (MQ=255)
                                                                                                                                        ccAGAGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGGCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTg               <  1:657143/150‑1 (MQ=255)
                                                                                                                                          gggaaCTTCCATTCTGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTa             >  1:728372/5‑150 (MQ=255)
                                                                                                                                          gggaaCTTCCATTCTGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTa             >  1:1022113/5‑150 (MQ=255)
                                                                                                                                           gTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTAt            <  2:1350087/150‑1 (MQ=255)
                                                                                                                                             gACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTc          >  2:982530/1‑150 (MQ=255)
                                                                                                                                               cTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcac        <  2:317685/150‑1 (MQ=255)
                                                                                                                                               cTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcac        <  1:492963/150‑1 (MQ=255)
                                                                                                                                               cTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcac        >  2:168372/1‑150 (MQ=255)
                                                                                                                                               cTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcac        >  1:567692/1‑150 (MQ=255)
                                                                                                                                               cTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcac        <  1:1306973/150‑1 (MQ=255)
                                                                                                                                               cTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGAGTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcac        <  1:843314/150‑1 (MQ=255)
                                                                                                                                                ttCCAGAATGGAAGTTCTCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGGTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcaca       <  1:822474/150‑1 (MQ=255)
                                                                                                                                                  ccAGACTGGCAGCTCCCGGACGCTACAGCCCCAGGCTGGTGTCGCGAGACTCCAGCACAATGTCGTACATTTCACGCTGCGCCTGGGTGAATTTGCCGGTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTATGTATTCACAAc     <  2:1729156/150‑1 (MQ=255)
                                                                                                                                                  ccAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGCCTCCAGCAGAATGTCGTAGATTTCACGCTGGGCCTGGGTGTATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTATCCTTTGTATTCACAAc     >  1:119999/1‑150 (MQ=255)
                                                                                                                                                  ccAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCGGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTCACAAc     <  2:102159/150‑1 (MQ=255)
                                                                                                                                                     gAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTCACAACCcg  >  1:1121662/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTCACAACCCG  >  NC_000913/3053758‑3054056

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