Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 519,341 C→G 5.6% H243D (CAT→GAT)  argC → N‑acetyl‑gamma‑glutamyl‑phosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947519,3410CG 5.6% 307.5 / 7.1 99H243D (CAT→GAT) argCN‑acetyl‑gamma‑glutamyl‑phosphate reductase
Reads supporting (aligned to +/- strand):  ref base C (48/45);  new base G (3/3);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 102 total aligned reads displayed.

TGTTCTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTCTGCCGC  >  NC_002947/519192‑519490
                                                                                                                                                     |                                                                                                                                                     
tGTTCTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATcc                                                                                                                                                       >  1:122408/1‑150 (MQ=255)
tGTTCTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATcc                                                                                                                                                       >  2:2398411/1‑150 (MQ=255)
 gTTCTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCa                                                                                                                                                      <  1:2406266/150‑1 (MQ=255)
  aacagCGAGCCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCAt                                                                                                                                                     <  1:1243336/146‑1 (MQ=255)
    cTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGc                                                                                                                                                   <  2:180913/150‑1 (MQ=255)
    cTGCGAGGCCGGCGAAAGCATGAAAGCCAAAGTCTTACGGGTGCATGGCCTCTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGc                                                                                                                                                   <  1:1738174/150‑1 (MQ=255)
     tGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCAt                                                                                                                                                     >  1:2230542/1‑147 (MQ=255)
     tGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCAt                                                                                                                                                     <  2:2230542/147‑1 (MQ=255)
        gAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATGGATGCCAcg                                                                                                                                               <  1:1846492/150‑2 (MQ=255)
        gAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCAcc                                                                                                                                               >  1:1408627/1‑150 (MQ=255)
        gAGGCCGGCGAAAGCATGAAAGCCTACCCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCAcc                                                                                                                                               <  2:2110953/150‑1 (MQ=255)
          ggCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCt                                                                                                                                             <  2:2085398/150‑1 (MQ=255)
              ggCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGTTTCGTGGCATCCATGCCACCCTGTAc                                                                                                                                         <  1:1645384/150‑1 (MQ=255)
                    aGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAAc                                                                                                                                   >  2:1012252/1‑150 (MQ=255)
                    aGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAAc                                                                                                                                   <  1:1012252/150‑1 (MQ=255)
                    aGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGAGGATCGGTGGCATCCATGCCACCCTGTACGCGAAc                                                                                                                                   >  2:1330413/1‑150 (MQ=255)
                       gtGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGAGCCTGCGCCTAGACGCAGACAAGGCCATCGGCCTGACTTTCATGCCGCACCTGACGCCTATGATCCGTGGCGTCCATACCACCCTTTACGCGAACGTg                                                                                                                                >  1:1552755/2‑150 (MQ=255)
                       aTGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTg                                                                                                                                >  2:2058418/1‑150 (MQ=255)
                         gAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGTGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGt                                                                                                                              <  2:1346371/150‑1 (MQ=255)
                         gAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGt                                                                                                                              >  1:652206/1‑150 (MQ=255)
                          aaaGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTc                                                                                                                             <  1:1521054/150‑1 (MQ=255)
                                  cgcgGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGATCGATACCTc                                                                                                                     >  2:1641996/1‑150 (MQ=255)
                                     tgTAGAGCGCGATCGCCACCGTCCGGTTATCAGCCAGGGCTTGCGCCTGACCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTGGGt                                                                                                                  >  2:344200/2‑150 (MQ=255)
                                     ggTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGt                                                                                                                  >  1:131515/1‑150 (MQ=255)
                                        aaaGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGGGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCGATGCCACTCGGTACGCGAACGTGGTCGATACCTCGGTCGa                                                                                                               >  1:691752/1‑150 (MQ=255)
                                        aaaGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATGGATGCCAcgga                                                                                                                                             <  2:1041762/120‑4 (MQ=255)
                                        aaaGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATGGATGCCAcgga                                                                                                                                             >  1:1041762/1‑117 (MQ=255)
                                             gCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGc                                                                                                          <  2:652206/150‑1 (MQ=255)
                                               cATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTgg                                                                                                                               >  1:1153512/1‑127 (MQ=255)
                                               cATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTgg                                                                                                                               <  2:1153512/127‑1 (MQ=255)
                                               cATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAg                                                                                                        <  1:147117/150‑1 (MQ=255)
                                                aTCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCt                                                                                                                      >  2:28492/1‑135 (MQ=255)
                                                aTCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCt                                                                                                                      <  1:28492/135‑1 (MQ=255)
                                                    ccACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGgcgcc                                                                                                   >  1:136293/1‑149 (MQ=255)
                                                          gCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCga                                                                                             <  1:2058418/150‑1 (MQ=255)
                                                           ccGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATGGATGCCAcg                                                                                                                                               >  1:748708/1‑98 (MQ=255)
                                                           ccGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATGGATGCCAcg                                                                                                                                               <  2:748708/99‑2 (MQ=255)
                                                           ccGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCATGTACGCGAACGTGGTAGATACCGCGCGCGACCTGCAGGCGCTGTTCgag                                                                                            >  2:1468711/1‑150 (MQ=255)
                                                             ggAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGaa                                                                                          >  1:533249/1‑150 (MQ=255)
                                                              gagaTCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGaaa                                                                                         >  1:1939497/1‑150 (MQ=255)
                                                                    aGCCAGGTCCTGCGCCTGGCCGCGGTCAATTACATCGCCCTGACTTTCGTGCCCCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTAc                                                                                   <  1:722876/150‑1 (MQ=255)
                                                                      ccAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGc                                                                                 <  1:1755972/150‑1 (MQ=255)
                                                                       cAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGcc                                                                                >  2:133631/1‑150 (MQ=255)
                                                                        aGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCg                                                                               >  1:552728/1‑150 (MQ=255)
                                                                             cTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGc                                                                                                          >  2:2156413/1‑118 (MQ=255)
                                                                             cTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGc                                                                                                          <  1:2156413/118‑1 (MQ=255)
                                                                                 gcCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGt                                                                      >  1:2151587/1‑150 (MQ=255)
                                                                                   cTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTtc                                                                    >  2:1347996/1‑150 (MQ=255)
                                                                                     ggCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTtcgt                                                                  >  2:353413/1‑150 (MQ=255)
                                                                                       ccGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTa                                                                                                                                          <  1:2332790/76‑1 (MQ=255)
                                                                                       ccGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTa                                                                                                                                          >  2:2332790/1‑76 (MQ=255)
                                                                                             gCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGa                                                          <  1:353413/150‑1 (MQ=255)
                                                                                              cAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGAt                                                         <  2:1408627/150‑1 (MQ=255)
                                                                                              cAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGAACTGCAGAGGTTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGAt                                                         <  1:2224921/150‑1 (MQ=255)
                                                                                               aaGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATg                                                        >  2:486312/1‑150 (MQ=255)
                                                                                               aaGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATg                                                        >  2:1822145/1‑150 (MQ=255)
                                                                                                  gACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCca                                                     >  2:2134514/1‑150 (MQ=255)
                                                                                                   aCATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTc                                                                                                                 <  1:266579/89‑1 (MQ=255)
                                                                                                   aCATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTc                                                                                                                 >  2:266579/1‑89 (MQ=255)
                                                                                                    cATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTa                                                                                    <  1:708646/117‑1 (MQ=255)
                                                                                                    cATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTa                                                                                    >  2:708646/1‑117 (MQ=255)
                                                                                                     aTCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCAt                                                                                                                                                     <  1:616918/51‑1 (MQ=255)
                                                                                                     aTCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCAt                                                                                                                                                     >  2:616918/1‑51 (MQ=255)
                                                                                                     aTCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCt                                                                                                            >  1:1883344/1‑92 (MQ=255)
                                                                                                     aTCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCt                                                                                                            <  2:1883344/92‑1 (MQ=255)
                                                                                                     aTCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCcagca                                                  >  2:876663/1‑150 (MQ=255)
                                                                                                     aTCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCcagca                                                  <  1:486312/150‑1 (MQ=255)
                                                                                                     aTCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTAAAACAGGCCGTTGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATAAATAGCCCGTCGTCGTGATGCcagca                                                  >  1:150674/1‑150 (MQ=255)
                                                                                                         gCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAAc                                                                                        >  1:2332031/1‑108 (MQ=255)
                                                                                                         gCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAAc                                                                                        <  2:2332031/108‑1 (MQ=255)
                                                                                                         gCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTg                                                           >  1:2114756/1‑137 (MQ=255)
                                                                                                         gCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTg                                                           <  2:2114756/137‑1 (MQ=255)
                                                                                                         gCCGGACTTTCGTGCCGCCCTTGATGCCGATGATCCGGGGCATCCATGCCACCCTGTAGGCGAACGTGGTCAATATCTCGGTCGACCTGCTGGCGCTGTTCGAGAACCGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggca                                              <  1:344200/150‑1 (MQ=255)
                                                                                                          ccTGACTTTCGTGCCGCACCTGACGACGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggcag                                             <  2:382290/150‑1 (MQ=255)
                                                                                                            tGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGcc                                           >  2:1091604/1‑150 (MQ=255)
                                                                                                                tttCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATcc                                       <  2:1883376/150‑1 (MQ=255)
                                                                                                                tttCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATcc                                       <  2:2312037/150‑1 (MQ=255)
                                                                                                                 ttCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCAATCCg                                      <  2:1010033/150‑1 (MQ=255)
                                                                                                                     tGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGaaa                                  <  1:133631/150‑1 (MQ=255)
                                                                                                                          cACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGt                             >  2:680975/1‑150 (MQ=255)
                                                                                                                           aCCTGACGCCGATTATCCGTGGCATCCATGCCGCCCTGTACGCGAACGTGGTGGATACCTCGGTCGACCTGCAGGCGCTGGTCGAGAAACACTACGCCGATGAACCGTTCGTCGACGTTATGCCAGCAGGCAGCCATCCGGAAACCCGTa                            <  2:691752/150‑1 (MQ=255)
                                                                                                                           aCCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCga                                                                                             >  1:488648/1‑85 (MQ=255)
                                                                                                                           aCCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCga                                                                                             <  2:488648/85‑1 (MQ=255)
                                                                                                                            ccTGACGCCGATGATCCGTGGCATCCATGCCACACTCGATGGGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAg                           <  1:1330413/150‑1 (MQ=255)
                                                                                                                             cTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAgc                          >  2:292291/1‑150 (MQ=255)
                                                                                                                              tGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTGGCGTTAGGGTCAGCAAAGGTGTCGGTGGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAgcg                         >  1:891746/1‑150 (MQ=255)
                                                                                                                                   ccGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCcagcag                                                 <  1:2113122/121‑1 (MQ=255)
                                                                                                                                   ccGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCcagcag                                                 >  2:2113122/1‑121 (MQ=255)
                                                                                                                                    cGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGcgtgcgtg                   >  1:1634011/1‑150 (MQ=255)
                                                                                                                                     gatgatCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCgtgg                  >  2:2324788/1‑150 (MQ=255)
                                                                                                                                        gatCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGc               <  2:2151587/150‑1 (MQ=255)
                                                                                                                                             gTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACg          >  1:2229985/1‑150 (MQ=255)
                                                                                                                                               ggCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGaa                                                                          >  1:384963/1‑84 (MQ=255)
                                                                                                                                               ggCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGaa                                                                          <  2:384963/84‑1 (MQ=255)
                                                                                                                                                gCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGtgt       <  1:2134514/150‑3 (MQ=255)
                                                                                                                                                gCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTcc       >  2:2310662/1‑149 (MQ=255)
                                                                                                                                                  aTCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTCTGc     <  1:1822145/150‑1 (MQ=255)
                                                                                                                                                  aTCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTCTGc     >  2:1863171/1‑150 (MQ=255)
                                                                                                                                                  aTCCATGCCACCCTGTACGCGAACATGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggcag                                             >  1:2126104/1‑110 (MQ=255)
                                                                                                                                                  aTCCATGCCACCCTGTACGCGAACATGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAgcaggcag                                             <  2:2126104/110‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTTCTGCGAGGCCGGCGAAAGCATGAAAGCCTACGCGGTAAAGGGCCATCGCCACTTGCCGGAGATCAGCCAGGGCCTGCGCCTGGCCGCGGGCAAGGACATCGGCCTGACTTTCGTGCCGCACCTGACGCCGATGATCCGTGGCATCCATGCCACCCTGTACGCGAACGTGGTCGATACCTCGGTCGACCTGCAGGCGCTGTTCGAGAAACGCTACGCCGATGAACCGTTCGTCGACGTGATGCCAGCAGGCAGCCATCCGGAAACCCGTAGCGTGCGTGGTGCCAACGTCTGCCGC  >  NC_002947/519192‑519490

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