Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,004,515 G→C G61G (GGC→GGG yigI ← 4HBT thioesterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,5150GC97.7% 319.1 / ‑0.6 88G61G (GGC→GGGyigI4HBT thioesterase family protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (44/42);  total (45/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

GTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATC  >  NC_000913/4004421‑4004611
                                                                                               |                                                                                                
gTAGTTCATCTTCACTGATGGTTGCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTg                                                                                             <  2:752517/100‑1 (MQ=255)
 tAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGc                                                                                            >  2:977456/1‑100 (MQ=255)
 tAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGc                                                                                            <  1:499116/100‑1 (MQ=255)
  aGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGc                                                                                            >  2:83850/1‑99 (MQ=255)
   gTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGAACCCCCCGTGCaa                                                                                          >  2:833153/1‑100 (MQ=255)
      cATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAAt                                                                                       >  1:259085/1‑100 (MQ=255)
        tCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGc                                                                                     >  2:1314136/1‑100 (MQ=255)
        tCTTCACTGATGGTTTCGTGACGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGc                                                                                     <  2:1192176/100‑1 (MQ=255)
         cTTCCCTGCTGGTTTCGTTGCGGGTTAAGGTGCTTCTTACTCACGCCAGAGCCGCGGTGTCCTCCCGCGCCGACGCAATGACCCCCCCGTGCaaaa                                                                                        >  2:1702675/1‑96 (MQ=255)
         cTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCt                                                                                    >  2:1000779/1‑100 (MQ=255)
             aCTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCCTTgc                                                                                >  1:251121/1‑100 (MQ=255)
              cTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTgcg                                                                               >  1:2090615/1‑100 (MQ=255)
               tGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTgcgc                                                                              <  2:2273698/100‑1 (MQ=255)
                gATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGcc                                                                             >  2:1061688/1‑100 (MQ=255)
                gATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGcc                                                                             <  1:1074135/100‑1 (MQ=255)
                 aTGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGccc                                                                            >  2:2006328/1‑100 (MQ=255)
                 aTGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGccc                                                                            <  1:491393/100‑1 (MQ=255)
                  tGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCa                                                                           >  2:993907/1‑100 (MQ=255)
                   ggTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAg                                                                          >  2:1790935/1‑100 (MQ=255)
                   ggTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAg                                                                          >  2:1972962/1‑100 (MQ=255)
                    gTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAg                                                                          <  1:135711/99‑1 (MQ=255)
                    gTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGt                                                                         >  1:1620406/1‑100 (MQ=255)
                      ttCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTg                                                                       >  1:1664236/1‑100 (MQ=255)
                       tCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGc                                                                      >  1:2162883/1‑100 (MQ=255)
                        cGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGcc                                                                     <  1:967201/100‑1 (MQ=255)
                         gTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGccc                                                                    <  2:70899/100‑1 (MQ=255)
                         gTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGccc                                                                    >  1:1371815/1‑100 (MQ=255)
                          tGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCcca                                                                   >  1:2166877/1‑100 (MQ=255)
                           ggCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccac                                                                  >  1:1345799/1‑100 (MQ=255)
                           ggCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccac                                                                  >  1:513853/1‑100 (MQ=255)
                             cGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccacca                                                                <  2:2162876/100‑1 (MQ=255)
                               ggTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccaccacc                                                              >  2:2027729/1‑98 (MQ=255)
                                 ttAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACcatcat                                                            >  1:30655/1‑100 (MQ=255)
                                 tcaaGGTACTTCGCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATTCCCCCCTCGTGCAAAATGCTTTGCGCCCAGTTGCccacc                                                                 <  2:1830406/93‑1 (MQ=255)
                                  tAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCAtt                                                           >  1:2178358/1‑100 (MQ=255)
                                   aaGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTg                                                          >  2:1879494/1‑100 (MQ=255)
                                     ggTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGc                                                        <  2:1345796/100‑1 (MQ=255)
                                       tACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCACAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTg                                                      >  1:333285/1‑100 (MQ=255)
                                       tACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTg                                                      >  1:1183647/1‑100 (MQ=255)
                                       tACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTg                                                      <  2:672734/100‑1 (MQ=255)
                                        aCTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGa                                                     >  2:1165073/1‑100 (MQ=255)
                                           tCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAACATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttt                                                  >  1:260671/1‑100 (MQ=255)
                                           tCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttt                                                  >  1:738346/1‑100 (MQ=255)
                                           tCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttt                                                  >  1:2099063/1‑100 (MQ=255)
                                              cACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttttt                                                >  1:2303540/1‑99 (MQ=255)
                                              cACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTa                                               <  2:1410358/100‑1 (MQ=255)
                                               aCGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTaa                                              <  1:1561544/100‑1 (MQ=255)
                                                ngCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTaaa                                             >  2:384128/2‑100 (MQ=255)
                                                 gCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAg                                            >  1:1665073/1‑100 (MQ=255)
                                                  cacaCCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAg                                            <  1:1689108/100‑1 (MQ=255)
                                                  cacaCCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAgg                                           >  1:271668/1‑100 (MQ=255)
                                                  cacaCCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAgg                                           >  2:1210490/1‑100 (MQ=255)
                                                    cacCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGcc                                         <  2:1888090/100‑1 (MQ=255)
                                                    cacCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGcc                                         <  2:1127032/100‑1 (MQ=255)
                                                      ccAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAg                                       >  1:135000/1‑100 (MQ=255)
                                                        aGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCt                                     >  1:26814/1‑100 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  <  2:815555/100‑1 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  <  1:1061690/100‑1 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  <  2:251119/100‑1 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  >  2:1774559/1‑100 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  <  1:1790940/100‑1 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGAAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  <  2:1048236/100‑1 (MQ=255)
                                                              gcggcgACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGa                               <  1:1847619/100‑1 (MQ=255)
                                                                 gcgACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAAct                            >  2:1716705/1‑100 (MQ=255)
                                                                 gcgACATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAAct                            >  2:390011/1‑100 (MQ=255)
                                                                     cATCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCttt                        >  1:2079105/1‑100 (MQ=255)
                                                                      aTCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCtttt                       <  1:2006335/100‑1 (MQ=255)
                                                                      aTCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCtttt                       >  1:686101/1‑100 (MQ=255)
                                                                      aTCCAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCtttt                       <  1:1662341/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                     <  2:2178351/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                     <  2:2099056/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                     <  2:738345/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                     <  1:1210493/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                     <  1:1590937/100‑1 (MQ=255)
                                                                          aGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGt                   <  2:2303531/100‑1 (MQ=255)
                                                                          aGCGCCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGt                   <  1:1879500/100‑1 (MQ=255)
                                                                             gcCGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAAc                >  1:1176508/1‑100 (MQ=255)
                                                                              ccGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACg               <  1:993909/100‑1 (MQ=255)
                                                                               cGACGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGc              <  1:2027736/100‑1 (MQ=255)
                                                                                 aCGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTc            <  2:1415816/100‑1 (MQ=255)
                                                                                  cGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTcc           <  1:1972969/100‑1 (MQ=255)
                                                                                  cGCAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTcc           <  1:1124498/100‑1 (MQ=255)
                                                                                    cAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAg         <  2:2079098/100‑1 (MQ=255)
                                                                                    cAATGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAg         >  1:2088177/1‑100 (MQ=255)
                                                                                      aTGGCCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGtt       <  2:2166870/100‑1 (MQ=255)
                                                                                      aTGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGtt       <  1:1702680/100‑1 (MQ=255)
                                                                                       tGACCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTc      <  1:1314139/100‑1 (MQ=255)
                                                                                         aCCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCa    <  2:1183644/100‑1 (MQ=255)
                                                                                           ccccccCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCa    <  2:978531/99‑1 (MQ=255)
                                                                                            ccccccGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATc  >  2:441635/1‑100 (MQ=255)
                                                                                               |                                                                                                
GTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATC  >  NC_000913/4004421‑4004611

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