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→GGG) yigI4HBT 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

GTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATC  >  NC_000913/4004421‑4004611
                                                                                               |                                                                                                
gTAGTTCATCTTCACTGATGGTTGCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTg                                                                                             <  2:752517/100‑1 (MQ=255)
 tAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGc                                                                                            >  2:977456/1‑100 (MQ=255)
 tAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGc                                                                                            <  1:499116/100‑1 (MQ=255)
  aGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGc                                                                                            >  2:83850/1‑99 (MQ=255)
   gTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑ACCCCCCGTGCaa                                                                                          >  2:833153/1‑100 (MQ=255)
      cATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAAt                                                                                       >  1:259085/1‑100 (MQ=255)
        tCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGc                                                                                     >  2:1314136/1‑100 (MQ=255)
        tCTTCACTGATGGTTTCGTGACGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGc                                                                                     <  2:1192176/100‑1 (MQ=255)
         cTTCCCTGCTGGTTTCGTTGCGGGTTAAGGTGCTTCTTACTCACGCCAGAGCCGCGGTGTCCTCCCGCGCCGACGCAATGA‑CCCCCCCGTGCaaaa                                                                                        >  2:1702675/1‑96 (MQ=255)
         cTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCt                                                                                    >  2:1000779/1‑100 (MQ=255)
             aCTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCCTTgc                                                                                >  1:251121/1‑100 (MQ=255)
              cTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTgcg                                                                               >  1:2090615/1‑100 (MQ=255)
               tGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTgcgc                                                                              <  2:2273698/100‑1 (MQ=255)
                gATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGcc                                                                             >  2:1061688/1‑100 (MQ=255)
                gATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGcc                                                                             <  1:1074135/100‑1 (MQ=255)
                 aTGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGccc                                                                            >  2:2006328/1‑100 (MQ=255)
                 aTGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGccc                                                                            <  1:491393/100‑1 (MQ=255)
                  tGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCa                                                                           >  2:993907/1‑100 (MQ=255)
                   ggTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAg                                                                          >  2:1790935/1‑100 (MQ=255)
                   ggTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAg                                                                          >  2:1972962/1‑100 (MQ=255)
                    gTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAg                                                                          <  1:135711/99‑1 (MQ=255)
                    gTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGt                                                                         >  1:1620406/1‑100 (MQ=255)
                      ttCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTg                                                                       >  1:1664236/1‑100 (MQ=255)
                       tCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGc                                                                      >  1:2162883/1‑100 (MQ=255)
                        cGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGcc                                                                     <  1:967201/100‑1 (MQ=255)
                         gTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGccc                                                                    <  2:70899/100‑1 (MQ=255)
                         gTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGccc                                                                    >  1:1371815/1‑100 (MQ=255)
                          tGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCcca                                                                   >  1:2166877/1‑100 (MQ=255)
                           ggCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccac                                                                  >  1:1345799/1‑100 (MQ=255)
                           ggCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccac                                                                  >  1:513853/1‑100 (MQ=255)
                             cGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccacca                                                                <  2:2162876/100‑1 (MQ=255)
                               ggTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCccaccacc                                                              >  2:2027729/1‑98 (MQ=255)
                                 ttAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACcatcat                                                            >  1:30655/1‑100 (MQ=255)
                                 tcaaGGTACTTCGCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATTC‑CCCCCTCGTGCAAAATGCTTTGCGCCCAGTTGCccacc                                                                 <  2:1830406/93‑1 (MQ=255)
                                  tAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCAtt                                                           >  1:2178358/1‑100 (MQ=255)
                                   aaGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTg                                                          >  2:1879494/1‑100 (MQ=255)
                                     ggTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGc                                                        <  2:1345796/100‑1 (MQ=255)
                                       tACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCACAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTg                                                      >  1:333285/1‑100 (MQ=255)
                                       tACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTg                                                      >  1:1183647/1‑100 (MQ=255)
                                       tACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTg                                                      <  2:672734/100‑1 (MQ=255)
                                        aCTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGa                                                     >  2:1165073/1‑100 (MQ=255)
                                           tCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAACATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttt                                                  >  1:260671/1‑100 (MQ=255)
                                           tCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttt                                                  >  1:738346/1‑100 (MQ=255)
                                           tCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttt                                                  >  1:2099063/1‑100 (MQ=255)
                                              cACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGAttttt                                                >  1:2303540/1‑99 (MQ=255)
                                              cACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTa                                               <  2:1410358/100‑1 (MQ=255)
                                               aCGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTaa                                              <  1:1561544/100‑1 (MQ=255)
                                                ngCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTaaa                                             >  2:384128/2‑100 (MQ=255)
                                                 gCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAg                                            >  1:1665073/1‑100 (MQ=255)
                                                  cacaCCAGACCGGCGGCGACATCCAGCGCCGACGCAATGACCCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAg                                            <  1:1689108/100‑1 (MQ=255)
                                                  cacaCCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAgg                                           >  1:271668/1‑100 (MQ=255)
                                                  cacaCCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAgg                                           >  2:1210490/1‑100 (MQ=255)
                                                    cacCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGcc                                         <  2:1888090/100‑1 (MQ=255)
                                                    cacCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGcc                                         <  2:1127032/100‑1 (MQ=255)
                                                      ccAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAg                                       >  1:135000/1‑100 (MQ=255)
                                                        aGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCt                                     >  1:26814/1‑100 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  <  2:815555/100‑1 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  <  1:1061690/100‑1 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  <  2:251119/100‑1 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  >  2:1774559/1‑100 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  <  1:1790940/100‑1 (MQ=255)
                                                           ccGGCGGCGACATCCAGCGCCGACGAAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCtgtg                                  <  2:1048236/100‑1 (MQ=255)
                                                              gcggcgACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGa                               <  1:1847619/100‑1 (MQ=255)
                                                                 gcgACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAAct                            >  2:1716705/1‑100 (MQ=255)
                                                                 gcgACATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAAct                            >  2:390011/1‑100 (MQ=255)
                                                                     cATCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCttt                        >  1:2079105/1‑100 (MQ=255)
                                                                      aTCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCtttt                       <  1:2006335/100‑1 (MQ=255)
                                                                      aTCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCtttt                       >  1:686101/1‑100 (MQ=255)
                                                                      aTCCAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCtttt                       <  1:1662341/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGA‑CCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                     <  2:2178351/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                     <  2:2099056/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                     <  2:738345/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                     <  1:1210493/100‑1 (MQ=255)
                                                                        ccAGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTc                     <  1:1590937/100‑1 (MQ=255)
                                                                          aGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGt                   <  2:2303531/100‑1 (MQ=255)
                                                                          aGCGCCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGt                   <  1:1879500/100‑1 (MQ=255)
                                                                             gcCGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAAc                >  1:1176508/1‑100 (MQ=255)
                                                                              ccGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACg               <  1:993909/100‑1 (MQ=255)
                                                                               cGACGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGc              <  1:2027736/100‑1 (MQ=255)
                                                                                 aCGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTc            <  2:1415816/100‑1 (MQ=255)
                                                                                  cGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTcc           <  1:1972969/100‑1 (MQ=255)
                                                                                  cGCAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTcc           <  1:1124498/100‑1 (MQ=255)
                                                                                    cAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAg         <  2:2079098/100‑1 (MQ=255)
                                                                                    cAATGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAg         >  1:2088177/1‑100 (MQ=255)
                                                                                      aTGG‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGtt       <  2:2166870/100‑1 (MQ=255)
                                                                                      aTGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGtt       <  1:1702680/100‑1 (MQ=255)
                                                                                       tGA‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTc      <  1:1314139/100‑1 (MQ=255)
                                                                                         a‑CCCCCCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCa    <  2:1183644/100‑1 (MQ=255)
                                                                                           ccccccCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCa    <  2:978531/99‑1 (MQ=255)
                                                                                            ccccccGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATc  >  2:441635/1‑100 (MQ=255)
                                                                                               |                                                                                                
GTAGTTCATCTTCACTGATGGTTTCGTGGCGGGTTAAGGTACTTCCCACGCACACCAGACCGGCGGCGACATCCAGCGCCGACGCAATGA‑CCCCGCCGTGCAAAATGCTTTGCGCCCAGTTGCCCACCATCATTGGCTGATTTTTAAAGGCCAGCTGTGCGAACTCTTTTTCGTAACGCTCCAGTTCCATC  >  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: ‑