Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,053,909 A→C 8.1% I311S (ATT→AGT) ‡ pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,053,9090AC 8.1% 253.5 / 10.3 95I311S (ATT→AGT) ‡pepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base A (50/37);  new base C (3/5);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTCACAACCCG  >  NC_000913/3053760‑3054056
                                                                                                                                                     |                                                                                                                                                   
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)
                                         ggACGATGGGCGTTCTGAGCGATCAGTACATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTT‑CCGGACGATACAGGCGCAGGCTGGTTTCgagag                                                                                                           >  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)
                                                                                                                                     aCCAGTGACTTCCAGAATGGTAGTT‑CCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTg               <  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:1022113/5‑150 (MQ=255)
                                                                                                                                        gggaaCTTCCATTCTGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTa             >  1:728372/5‑150 (MQ=255)
                                                                                                                                         gTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTAt            <  2:1350087/150‑1 (MQ=255)
                                                                                                                                           gACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTc          >  2:982530/1‑150 (MQ=255)
                                                                                                                                             cTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcac        <  1:1306973/150‑1 (MQ=255)
                                                                                                                                             cTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcac        <  2:317685/150‑1 (MQ=255)
                                                                                                                                             cTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcac        <  1:492963/150‑1 (MQ=255)
                                                                                                                                             cTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcac        >  1:567692/1‑150 (MQ=255)
                                                                                                                                             cTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTcac        >  2:168372/1‑150 (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)
                                                                                                                                                     |                                                                                                                                                   
CATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCTGAGCGATCAGTTCATCAACATCACCTTTCAGGATGCCGAGTTTTACCAGGCCGCTAACCATGATGCGCACCACTTCACCAGTGACTTCCAGAATGGAAGTTCCCGGACGATACAGGCGCAGGCTGGTTTCGAGAGACTCCAGCACAATGTCGTAGATTTCACGCTGGGCCTGGGTGAATTTGCCGTTGACCGGGAAGGTGCGGGTAATATCGCCAGCGTAACCTTTGTATTCACAACCCG  >  NC_000913/3053760‑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: ‑