Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,052,257 C→T 16.8% intergenic (+17/+32) gnsA → / ← yccM putative phosphatidylethanolamine synthesis regulator/putative 4Fe‑4S membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,052,2570CT16.8% 293.8 / 45.7 125intergenic (+17/+32)gnsA/yccMputative phosphatidylethanolamine synthesis regulator/putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base C (57/46);  new base T (9/12);  total (67/58)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

Only 100 of 135 total aligned reads displayed.

AGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCC  >  NC_000913/1052162‑1052356
                                                                                               |                                                                                                   
agTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTg                                                                                                 >  2:3773949/1‑100 (MQ=255)
   ctctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGa                                                                                              >  1:2491775/1‑100 (MQ=255)
     ctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTCTGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGa                                                                                              >  1:43914/1‑98 (MQ=255)
     ctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAgg                                                                                            >  2:1259627/1‑100 (MQ=255)
     ctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAgg                                                                                            >  2:3666892/1‑100 (MQ=255)
     ctGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAgg                                                                                            >  2:2761410/1‑100 (MQ=255)
       gAAATTCGCTTCACCGCACAAAAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGt                                                                                           <  1:2254992/99‑1 (MQ=255)
        aaaTTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGa                                                                                         <  2:3740824/100‑1 (MQ=255)
         aaTTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGAt                                                                                        <  2:538327/100‑1 (MQ=255)
         aaTTCGCTTCACCGCAAGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGAt                                                                                        <  1:3321556/100‑1 (MQ=255)
          aTTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATg                                                                                       <  2:1621159/100‑1 (MQ=255)
             cGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTc                                                                                    <  1:3541885/100‑1 (MQ=38)
             cGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTc                                                                                    <  1:115071/100‑1 (MQ=38)
             cGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAAAGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTc                                                                                    >  1:1866286/1‑100 (MQ=255)
             cGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATAAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTc                                                                                    >  2:906255/1‑100 (MQ=255)
              gCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTca                                                                                   <  2:80732/100‑2 (MQ=38)
              gCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGACGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTcc                                                                                   >  2:1793252/1‑100 (MQ=255)
              gCTTCACCGCACGAGAAAAAATGACCGGGCTTAAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTc                                                                                    >  2:1755992/1‑99 (MQ=255)
              gCTTCACCGCACGAGAAAAAATGAACGTGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTca                                                                                   <  1:261076/100‑2 (MQ=37)
               cTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTcac                                                                                  >  2:3413006/1‑98 (MQ=38)
               cTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTcac                                                                                  <  2:3587956/100‑3 (MQ=38)
               cTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTcac                                                                                  <  2:1253531/100‑3 (MQ=38)
               cTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCg                                                                                  >  2:32129/1‑100 (MQ=255)
               cTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCg                                                                                  <  2:1462685/100‑1 (MQ=255)
               cTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCg                                                                                  >  2:134057/1‑100 (MQ=255)
                ttCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTcaca                                                                                 <  2:1495673/100‑4 (MQ=38)
                ttCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCTCGCACGGGGTGATGTTcaca                                                                                 <  1:3412986/100‑4 (MQ=37)
                ttCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGc                                                                                 <  1:2952159/100‑1 (MQ=255)
                 tCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCt                                                                                <  2:2491787/100‑1 (MQ=255)
                   aCCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGTACATCACCCCGTGCGAGGTGATGTTCCGCttg                                                                              <  1:208636/100‑1 (MQ=255)
                     cGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCttgtt                                                                            <  2:1277086/100‑1 (MQ=255)
                       cACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGc                                                                          <  1:1731025/100‑1 (MQ=255)
                         cGAGAAAAACTGACCGGGCTTGAAAGATATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCt                                                                         >  1:1828381/1‑99 (MQ=255)
                         cGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTa                                                                        >  2:1281376/1‑100 (MQ=255)
                         cGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTa                                                                        <  1:1072119/100‑1 (MQ=255)
                          gagaAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTaa                                                                       >  2:2745516/1‑100 (MQ=255)
                            ggaaaaaTGACCGGGCTTTAAAGTTATAATGTCAAATTCAAAATAATGTGATTTTGTGAACATCTCCTCGTACGGGGTGATGTTcaca                                                                                 <  1:3276781/86‑4 (MQ=25)
                            gaAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAAtt                                                                     >  1:990856/1‑100 (MQ=255)
                            gaAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAAtt                                                                     >  2:3190298/1‑100 (MQ=255)
                             aaaaaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAAttt                                                                    <  1:2828955/100‑1 (MQ=255)
                              aaaaaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTa                                                                   <  2:3501072/100‑1 (MQ=255)
                                 aaTGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTg                                                                >  1:990905/1‑100 (MQ=255)
                                   tGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGAc                                                              <  2:278002/100‑1 (MQ=255)
                                   tGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGAc                                                              >  1:1607043/1‑100 (MQ=255)
                                     aCCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCa                                                            <  2:2074018/100‑1 (MQ=255)
                                      ccGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGAcc                                                             >  1:2583279/1‑98 (MQ=255)
                                      ccGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCaa                                                           >  2:2703041/1‑100 (MQ=255)
                                           cTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTTTGCACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGCCCAATCAtt                                                      <  1:2276506/100‑1 (MQ=255)
                                           cTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCAtt                                                      >  2:1301134/1‑100 (MQ=255)
                                           cTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCAtt                                                      >  1:3300445/1‑100 (MQ=255)
                                              gAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGgc                                                   <  2:1353873/100‑1 (MQ=255)
                                              gAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGgc                                                   <  2:709118/100‑1 (MQ=255)
                                              gAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGGGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGgc                                                   <  1:2062554/100‑1 (MQ=255)
                                               aaaGTTATGATTTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGgc                                                   <  2:3807907/99‑1 (MQ=255)
                                               aaaGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCACGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGgcg                                                  >  1:1816044/1‑100 (MQ=255)
                                                aaGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGgcgc                                                 >  2:2233623/1‑100 (MQ=255)
                                                    tATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTtgt                                             >  2:224810/1‑100 (MQ=255)
                                                       gttcacaaaatcacattATTTTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGa                                          >  2:2930390/6‑100 (MQ=38)
                                                        ttcacaaaatcacattATTTTGATTTTGTGAACATCACCCCGTGCGAGGCGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGaa                                         >  2:1477336/5‑100 (MQ=37)
                                                         tgTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGaa                                         >  2:1925395/1‑99 (MQ=255)
                                                           tcaaaatcaaaaTAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTa                                      <  1:2658049/100‑1 (MQ=255)
                                                           tcaaaatcaaaaTAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTa                                      >  2:3002054/1‑100 (MQ=255)
                                                            cgaaatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTaa                                     >  1:3689493/3‑100 (MQ=255)
                                                             aaaatcaaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAg                                    <  2:3105537/100‑1 (MQ=255)
                                                                  caaaatAATGTGATTTTGTGAACATTACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCg                               >  1:3153323/1‑100 (MQ=255)
                                                                   aaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCgg                              <  2:1030339/100‑1 (MQ=255)
                                                                   aaaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCgg                              <  1:3044741/100‑1 (MQ=255)
                                                                    aaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGt                             >  1:1968740/1‑100 (MQ=255)
                                                                    aaatAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGt                             <  1:1016286/100‑1 (MQ=255)
                                                                     aatagtGGGATTTTGTGAACATGACGCCGTGCGGGGTGATGGTCCGCTGGTGGTTAAGTTGGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCTGTa                            >  1:2221593/1‑100 (MQ=255)
                                                                         atGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAAt                        >  2:1423612/1‑100 (MQ=255)
                                                                         atGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAAt                        >  2:2626280/1‑100 (MQ=255)
                                                                         atGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAAt                        <  2:1750044/100‑1 (MQ=255)
                                                                           gtgATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtc                      >  1:209143/1‑100 (MQ=255)
                                                                            tgATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcc                     >  2:3054138/1‑100 (MQ=255)
                                                                            tgATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcc                     <  2:874029/100‑1 (MQ=255)
                                                                             gATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcct                    <  2:3605391/100‑1 (MQ=255)
                                                                             gATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCCAATCGATCGCCCAATCATTGGCGCTTGTGGAATTCAGCGTCGGTACAATtcc                     >  1:3304198/1‑99 (MQ=255)
                                                                              aTTTTGTGAACATCACCCCGTGCGAGGTGATGTTCTGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcctc                   >  2:2376603/1‑100 (MQ=255)
                                                                               ttttGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATtcctcc                  <  2:2310778/100‑1 (MQ=255)
                                                                                tttGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCg                 <  1:948907/100‑1 (MQ=255)
                                                                                tttGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCg                 >  2:532691/1‑100 (MQ=255)
                                                                                 ttGTGAACATCACCCCGTGCGAGGTGATGTTCTGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGTAATTAAGCGTCGGTACAATTCCTCCg                 <  1:3663431/99‑1 (MQ=255)
                                                                                 ttGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCgg                <  2:3333047/100‑1 (MQ=255)
                                                                                  tgtgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGc               >  1:1974168/1‑100 (MQ=255)
                                                                                  gcggAACATCACCTCGCACGGGGGGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGc               >  2:1257568/4‑100 (MQ=37)
                                                                                   gtgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCa              >  1:1438488/1‑100 (MQ=255)
                                                                                    tgAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCAc             >  2:2660109/1‑100 (MQ=255)
                                                                                    ggAACATCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCAc             >  2:2890221/2‑100 (MQ=38)
                                                                                     gAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAACTTAGTGACCAATCATTTGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCAcc            >  1:2554633/1‑100 (MQ=255)
                                                                                       aCATCACCTCGCACGGGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCgg          >  2:2324340/1‑100 (MQ=38)
                                                                                        cATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCggg         >  1:3192563/1‑100 (MQ=255)
                                                                                            aCCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCttt     >  2:3214271/1‑100 (MQ=255)
                                                                                             ccccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCTCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTg    >  2:764067/1‑100 (MQ=255)
                                                                                             ccccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTg    >  1:1518312/1‑100 (MQ=255)
                                                                                              cccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGc   >  1:3460568/1‑100 (MQ=255)
                                                                                              cccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGc   <  2:3605635/100‑1 (MQ=255)
                                                                                              cccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGc   <  1:2626267/100‑1 (MQ=255)
                                                                                              cccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGc   >  2:2508954/1‑100 (MQ=255)
                                                                                               ccGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGcc  <  2:2315453/100‑1 (MQ=255)
                                                                                               |                                                                                                   
AGTCTCTGAAATTCGCTTCACCGCACGAGAAAAAATGACCGGGCTTGAAAGTTATGATGTCAAAATCAAAATAATGTGATTTTGTGAACATCACCCCGTGCGAGGTGATGTTCCGCTTGTTGCTAATTTAGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCC  >  NC_000913/1052162‑1052356

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: