Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,466,858 G→C A45G (GCC→GGC)  bfd ← bacterioferritin‑associated ferredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,466,8580GC97.1% 357.5 / 0.4 102A45G (GCC→GGC) bfdbacterioferritin‑associated ferredoxin
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base C (53/46);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 102 total aligned reads displayed.

ACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTT  >  NC_000913/3466763‑3466955
                                                                                               |                                                                                                 
aCGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCg                                                                                               >  1:1965757/1‑100 (MQ=255)
aCGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCg                                                                                               >  1:1362569/1‑100 (MQ=255)
 cGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCgg                                                                                              <  2:1007484/100‑1 (MQ=255)
   gCTTACGAGTAAGTCAAAAAACACCCTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGac                                                                                            <  1:1631578/100‑1 (MQ=255)
    cTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGaca                                                                                           <  1:186915/100‑1 (MQ=255)
     ttACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACACCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGacac                                                                                          <  1:406274/100‑1 (MQ=255)
      tACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGAcaca                                                                                         <  2:55785/100‑1 (MQ=255)
      tACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGAcaca                                                                                         <  2:688032/100‑1 (MQ=255)
        cGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACAtt                                                                                       >  1:1573535/1‑100 (MQ=255)
          aGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTa                                                                                     >  1:702253/1‑100 (MQ=255)
             aaGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCg                                                                                  >  2:663731/1‑100 (MQ=255)
              aGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGc                                                                                 >  2:6423/1‑100 (MQ=255)
              aGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGc                                                                                 >  1:1597802/1‑100 (MQ=255)
              aGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGc                                                                                 <  2:1163156/100‑1 (MQ=255)
              aGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGc                                                                                 <  2:1597790/100‑1 (MQ=255)
               gTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCa                                                                                >  1:1175452/1‑100 (MQ=255)
                 cAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACt                                                                              <  2:760351/100‑1 (MQ=255)
                 cAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACt                                                                              >  2:106400/1‑100 (MQ=255)
                 cAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACt                                                                              >  1:1660397/1‑100 (MQ=255)
                 cAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACt                                                                              >  1:1892824/1‑100 (MQ=255)
                  aaaaaaCACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTg                                                                             <  2:842351/100‑1 (MQ=255)
                   aaaaaCACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTg                                                                             >  1:152155/1‑99 (MQ=255)
                    aaaaCACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGAt                                                                           >  2:23880/1‑100 (MQ=255)
                       acacacTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGc                                                                        <  1:160049/100‑1 (MQ=255)
                       acacacTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGc                                                                        <  1:550097/100‑1 (MQ=255)
                          cacTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGa                                                                     >  2:1687452/1‑100 (MQ=255)
                            cTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGAcc                                                                   >  2:395858/1‑100 (MQ=255)
                              ccGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCgg                                                                 >  1:1161551/1‑100 (MQ=255)
                                gTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGaa                                                               <  1:140775/100‑1 (MQ=255)
                                   tATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATaa                                                            >  1:1190509/1‑100 (MQ=255)
                                   tATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATaa                                                            >  1:1533015/1‑100 (MQ=255)
                                    aTGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATaaa                                                           >  1:943421/1‑100 (MQ=255)
                                    aTGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATaaa                                                           >  1:943259/1‑100 (MQ=255)
                                     tGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAt                                                          >  1:1926992/1‑100 (MQ=255)
                                      gCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtt                                                         <  2:321206/100‑1 (MQ=255)
                                      gCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtt                                                         >  2:337341/1‑100 (MQ=255)
                                      gCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtt                                                         >  1:1105687/1‑100 (MQ=255)
                                       cGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAttt                                                        <  2:1161545/100‑1 (MQ=255)
                                        ggACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCCCTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAttt                                                        <  2:1104864/99‑1 (MQ=255)
                                        ggACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAtttt                                                       <  2:1573525/100‑1 (MQ=255)
                                         gACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAAttttt                                                      <  1:106401/100‑1 (MQ=255)
                                            tCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTa                                                   <  2:1892806/100‑1 (MQ=255)
                                             ccTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTaa                                                  <  2:1783900/100‑1 (MQ=255)
                                              cTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTaa                                                  <  2:770221/99‑1 (MQ=255)
                                              cTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAt                                                 >  1:1823152/1‑100 (MQ=255)
                                              cTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAt                                                 >  1:945625/1‑100 (MQ=255)
                                              cTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAt                                                 >  2:697158/1‑100 (MQ=255)
                                               ttAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAAtt                                                >  2:2019447/1‑100 (MQ=255)
                                                tAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGTACACATTTACTGCACTGATTGCCGACCGGAATAAATTTTTTTAATTg                                               <  2:718061/100‑1 (MQ=255)
                                                tAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTg                                               >  1:747352/1‑100 (MQ=255)
                                                 aaaCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGc                                              <  1:1816201/100‑1 (MQ=255)
                                                  aaCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCt                                             >  2:944576/1‑100 (MQ=255)
                                                   tcTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTg                                            <  1:28890/99‑1 (MQ=255)
                                                   aCTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCt                                             <  2:712923/99‑1 (MQ=255)
                                                    cTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTgg                                           <  1:399147/100‑1 (MQ=255)
                                                    cTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTgg                                           <  2:1926974/100‑1 (MQ=255)
                                                    cTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTga                                           >  2:28890/1‑99 (MQ=255)
                                                      ccGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGaa                                         >  1:622373/1‑100 (MQ=255)
                                                      ccGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCCGCATGAACATATTTACCGCACTGATTGCCGACCGGAATACATTTTTTTAAtt                                                <  2:271807/93‑1 (MQ=255)
                                                       cGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGTCGACCGGAATAAATTTTTTTAATTGCTGGAAc                                        <  1:1317826/100‑1 (MQ=255)
                                                          cAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAg                                     >  2:530268/1‑100 (MQ=255)
                                                           aGCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGAAACGAGt                                    >  2:1101974/1‑100 (MQ=255)
                                                            gCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTg                                   <  2:1208658/100‑1 (MQ=255)
                                                            gCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCCCTGATTGCCGCCCGGAATAAATTTTTTTAATTGCTGGAACGAGTg                                   >  1:650261/1‑100 (MQ=255)
                                                            gCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTg                                   <  2:1190503/100‑1 (MQ=255)
                                                            gCTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCCACCGGAATAAATTTTTTTAATTGCTGGAACGAGTg                                   >  2:1770050/1‑100 (MQ=255)
                                                             cTGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGa                                  <  2:1965739/100‑1 (MQ=255)
                                                              tGCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAg                                 <  2:171188/100‑1 (MQ=255)
                                                               gCATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAgg                                >  1:1007760/1‑100 (MQ=255)
                                                                cATTAATTCATCCTCCATCACTTCCCGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAggg                               <  1:1665027/100‑1 (MQ=255)
                                                                cATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAggg                               <  2:2041915/100‑1 (MQ=255)
                                                                cATTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAggg                               >  2:1520533/1‑100 (MQ=255)
                                                                 aTTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAgggg                              <  2:1660384/100‑1 (MQ=255)
                                                                 aTTAATTCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAgggg                              >  1:951762/1‑100 (MQ=255)
                                                                 aTTAATTCATCCTCCATCACTGCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAgggg                              <  1:550331/100‑1 (MQ=255)
                                                                      ttCATCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAAc                         <  1:6423/100‑1 (MQ=255)
                                                                         aTCCTCCATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTgg                      >  2:937123/1‑100 (MQ=255)
                                                                          tcctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgc                     <  1:23880/100‑1 (MQ=255)
                                                                          tcctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgc                     <  2:271789/100‑1 (MQ=255)
                                                                          tcctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgc                     >  2:521520/1‑100 (MQ=255)
                                                                           cctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcg                    >  1:1610484/1‑100 (MQ=255)
                                                                            cttccgatcTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcgc                   <  2:390714/92‑1 (MQ=255)
                                                                            ctccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGgcgc                   >  1:273938/1‑100 (MQ=255)
                                                                              ccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGcac                 >  1:151475/1‑100 (MQ=255)
                                                                              ccATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGcac                 <  2:1175446/100‑1 (MQ=255)
                                                                               cATCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGcaca                >  1:1087986/1‑100 (MQ=255)
                                                                                 tCACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGc              >  1:921986/1‑100 (MQ=255)
                                                                                  cACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGcc             <  2:1007755/100‑1 (MQ=255)
                                                                                  cACTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGcc             >  2:749190/1‑100 (MQ=255)
                                                                                   aCTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCt            >  1:240372/1‑100 (MQ=255)
                                                                                   aCTTCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCt            <  2:704135/100‑1 (MQ=255)
                                                                                      tCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGcagatcggaag         >  1:390715/1‑91 (MQ=255)
                                                                                      tCACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGa          <  2:650259/99‑1 (MQ=255)
                                                                                       naCGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTg           >  2:1849056/2‑97 (MQ=255)
                                                                                       cACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACg        <  2:1362560/100‑1 (MQ=255)
                                                                                       cACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACg        <  1:697160/100‑1 (MQ=255)
                                                                                       cACGCGCGCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACg        >  2:1928799/1‑100 (MQ=255)
                                                                                         cgcgcgCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGaa      >  2:335525/1‑100 (MQ=255)
                                                                                         cgcgcgCCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGaa      >  2:224223/1‑100 (MQ=255)
                                                                                             cgGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAAtttt  <  2:1040892/100‑1 (MQ=255)
                                                                                               |                                                                                                 
ACGGCTTACGAGTAAGTCAAAAAACACACTCCGTTTATGCGGACTCCTTAAACTCCGGCAGCTGCATTAATTCATCCTCCATCACTTCACGCGCGGCCCGGACACATTTACCGCACTGATTGCCGACCGGAATAAATTTTTTTAATTGCTGGAACGAGTGAGGGGAAAACTGGCGCACAGCCTGACGAATTTT  >  NC_000913/3466763‑3466955

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