Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,784,169 G→A 21.7% G380D (GGT→GAT)  fadK → short chain acyl‑CoA synthetase, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,784,1690GA21.7% 167.4 / 38.5 87G380D (GGT→GAT) fadKshort chain acyl‑CoA synthetase, anaerobic
Reads supporting (aligned to +/- strand):  ref base G (35/33);  new base A (9/10);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

TCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATGAAGAAGGCTGGTA  >  NC_000913/1784070‑1784268
                                                                                                   |                                                                                                   
tcGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAgg                                                                                                     <  1:1976790/100‑1 (MQ=255)
    tgatCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGc                                                                                                 >  2:1671005/1‑100 (MQ=255)
    tgatCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGc                                                                                                 <  2:2494173/100‑1 (MQ=255)
    tgatCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGc                                                                                                 <  1:1716061/100‑1 (MQ=255)
        ccTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAg                                                                                             <  2:1872205/100‑1 (MQ=255)
         cTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAgg                                                                                            <  1:74732/100‑1 (MQ=255)
          tgtGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGt                                                                                           <  1:1574407/98‑1 (MQ=255)
           ttGTCGCGCTTTATGCACACCGATGGTTACTCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTg                                                                                          <  1:577678/100‑1 (MQ=255)
           ttGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTg                                                                                          <  1:8795/100‑1 (MQ=255)
           ttGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTg                                                                                          >  1:167057/1‑100 (MQ=255)
               cgcgCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGCTTAAAGTGGTCGATGACGCACGCACGACCTTACCGCCAGGTTGCGAAGGTgaa                                                                                        >  1:1956909/1‑98 (MQ=255)
                  gcTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGc                                                                                   >  1:545422/1‑100 (MQ=255)
                   cTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGcc                                                                                  <  2:1096257/100‑1 (MQ=255)
                     ttATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTc                                                                                >  1:1577806/1‑100 (MQ=255)
                        tGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTcgcg                                                                             >  1:2277713/1‑100 (MQ=255)
                          cacaCCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTg                                                                           >  2:2402171/1‑100 (MQ=255)
                           acacCGATGGTTACGCTGCCGCAGGTGTAGAGACGAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTg                                                                           <  2:1173557/99‑1 (MQ=255)
                             acCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGcc                                                                        <  2:1430208/100‑1 (MQ=255)
                              ccGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGccc                                                                       <  2:2413265/100‑1 (MQ=255)
                               cGATGGTTCCGCTGCCGCAGGTGTAGGGATTTAAGTGGTCGATGACGCACGTAAGAGCTTGCCGCCAGGTTGCGGAGGTGATGAAGCCTCGCGGGGcc                                                                        <  1:70134/98‑1 (MQ=255)
                                gATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTTGTCGATGATGTACGCAAGAGCTTACCGCCAGGGTGCGAAGGTGAAGAAGCCTCGCGTGGcct                                                                       <  1:1326113/98‑2 (MQ=255)
                                  tGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAt                                                                   >  1:567273/1‑100 (MQ=255)
                                  tGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAt                                                                   <  1:557015/100‑1 (MQ=255)
                                  tGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAt                                                                   >  2:569134/1‑100 (MQ=255)
                                   gggtACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAtg                                                                  <  1:252237/97‑1 (MQ=255)
                                    gTTACGCTGCCGCAGGTGTAGAGCTTAAAGTGGGCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAACATCAAGTGGCCCCAAtgt                                                                 >  1:1722180/1‑100 (MQ=255)
                                     ttACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAtgtg                                                                <  1:2046290/100‑1 (MQ=255)
                                     ttACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAAtgtg                                                                >  2:209773/1‑100 (MQ=255)
                                        cGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGttt                                                             >  1:1979392/1‑100 (MQ=255)
                                        cGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGttt                                                             >  2:1480748/1‑100 (MQ=255)
                                         gCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTa                                                            <  1:1366714/100‑1 (MQ=255)
                                            gCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCATGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATgg                                                         <  2:1900582/100‑1 (MQ=255)
                                            gCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATgg                                                         <  2:2093288/100‑1 (MQ=255)
                                            gCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATgg                                                         >  1:398478/1‑100 (MQ=255)
                                            gCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATg                                                          <  2:2404167/99‑1 (MQ=255)
                                             ccGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATggg                                                        <  2:413415/100‑1 (MQ=255)
                                             ccGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCGCGCGTGGCCCCAATGTGTATATggg                                                        >  1:1948170/1‑100 (MQ=255)
                                              ctttccagtAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATgggg                                                       <  1:1140473/93‑1 (MQ=255)
                                              cGCAGGTGTAGAAAACAGAGCTGTCCATATCACACGCAAAACACTACCGCCAGGTTGCGAAGGTCAAAAAACCACGCGTGGCCCCAACGTGTATATgtg                                                        >  1:1057493/1‑97 (MQ=255)
                                               gCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATgggg                                                       <  1:66160/99‑1 (MQ=255)
                                               gCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGAGTTTATgggg                                                       >  1:289510/1‑99 (MQ=255)
                                                  ggTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTAtt                                                   >  1:1893885/1‑100 (MQ=255)
                                                   gtgtAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTAttt                                                  <  2:2223291/100‑1 (MQ=255)
                                                    tgtAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTAtttt                                                 <  1:1174467/100‑1 (MQ=255)
                                                    tgtAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTAtttt                                                 >  1:835742/1‑100 (MQ=255)
                                                      tAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTtga                                               >  2:2344746/1‑100 (MQ=255)
                                                      tAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTtga                                               <  2:786667/100‑1 (MQ=255)
                                                       agagATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTtgat                                              >  2:1695245/1‑100 (MQ=255)
                                                       agagATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTtgat                                              <  2:1689262/100‑1 (MQ=255)
                                                       agagATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTtgat                                              >  2:1371909/1‑100 (MQ=255)
                                                       agagATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTtgat                                              >  2:1606420/1‑100 (MQ=255)
                                                        gagaTTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTtgatg                                             <  2:1069917/100‑1 (MQ=255)
                                                          gaTTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGaa                                           >  2:2186632/1‑100 (MQ=255)
                                                            ttAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAAcc                                         >  2:2397127/1‑100 (MQ=255)
                                                              aaaGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGGGGTTATGGGGTATTTTGATGAACCTg                                       >  2:647838/1‑100 (MQ=255)
                                                               aaGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTATATGGGGTATTTTGATGAAACTGa                                      >  1:2096230/1‑100 (MQ=255)
                                                                aGTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGaa                                     <  1:1518742/100‑1 (MQ=255)
                                                                 gTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAAt                                    <  1:1352886/100‑1 (MQ=255)
                                                                 gTGGTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAAt                                    >  2:1960201/1‑100 (MQ=255)
                                                                   ggTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTa                                  <  2:1174908/100‑1 (MQ=255)
                                                                   ggTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTa                                  <  2:127483/100‑1 (MQ=255)
                                                                   ggTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTa                                  <  1:379619/100‑1 (MQ=255)
                                                                    gTCGATGACGCACGCAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTaa                                 >  2:373562/1‑100 (MQ=255)
                                                                     tCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAAc                                >  2:1380015/1‑100 (MQ=255)
                                                                     tCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAAc                                <  1:1538042/100‑1 (MQ=255)
                                                                      cGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAAcc                               >  2:2230717/1‑100 (MQ=255)
                                                                       gATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCg                              >  1:1762584/1‑100 (MQ=255)
                                                                       gATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCg                              >  1:399186/1‑100 (MQ=255)
                                                                        aTGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCTTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAAGTaa                                 >  1:551907/1‑96 (MQ=255)
                                                                           acgcacgcAAGACCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCg                          >  1:535401/1‑100 (MQ=255)
                                                                            cgcacgcaAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGt                         <  1:555993/100‑1 (MQ=255)
                                                                             gcacgcaAGACCTTACCGCCAGGTTGCGTAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGt                         <  1:542983/99‑1 (MQ=255)
                                                                                cgcaAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGccc                     >  2:28416/1‑100 (MQ=255)
                                                                                  caAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTg                   <  2:1622433/100‑1 (MQ=255)
                                                                                  caAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTg                   >  2:158363/1‑100 (MQ=255)
                                                                                   aaGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTgg                  <  2:1762593/100‑1 (MQ=255)
                                                                                     gACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGAt                <  1:2402156/100‑1 (MQ=255)
                                                                                      aCCTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATg               <  1:1865832/100‑1 (MQ=255)
                                                                                        cTTACCGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATgaa             <  1:1878939/100‑1 (MQ=255)
                                                                                         ttACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATgaag            <  2:1930674/100‑1 (MQ=255)
                                                                                            cccccAGATTGCGAAGGTGAAGATGCCGCGCGTGGCCCCCATGCGTATATGGGGTCTTTTGATGAACCGGAATTAACCGCCCGTGCCCTGGATgaaaacg         >  1:808660/4‑96 (MQ=255)
                                                                                            ccGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATGaagaag         >  2:1075044/1‑100 (MQ=255)
                                                                                            ccGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCATGGATGaagaag         >  1:2266023/1‑100 (MQ=255)
                                                                                            ccGCCAGATTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATGaagaag         >  1:2209197/1‑100 (MQ=255)
                                                                                              gCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATGAAGAAGGc       >  1:923727/1‑100 (MQ=255)
                                                                                              gCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATGAAGAAGGc       >  1:1077362/1‑100 (MQ=255)
                                                                                                cAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATGAAGAAGGCTg     >  2:2265472/1‑100 (MQ=255)
                                                                                                 aGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATGAAGAAGGCTgg    <  2:2531753/100‑1 (MQ=255)
                                                                                                   gTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATGAAGAAGGCTGGTa  <  1:1380007/100‑1 (MQ=255)
                                                                                                   gTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATGAAGAAGGCTGGTa  <  1:1348673/100‑1 (MQ=255)
                                                                                                   gTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATGAAGAAGGCTGGTa  <  2:632137/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TCGATGATCCTTTGTCGCGCTTTATGCACACCGATGGTTACGCTGCCGCAGGTGTAGAGATTAAAGTGGTCGATGACGCACGCAAGACCTTACCGCCAGGTTGCGAAGGTGAAGAAGCCTCGCGTGGCCCCAATGTGTTTATGGGGTATTTTGATGAACCTGAATTAACCGCCCGTGCCCTGGATGAAGAAGGCTGGTA  >  NC_000913/1784070‑1784268

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