Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,467,120 A→G intergenic (‑129/+40) bfd ← / ← chiA bacterioferritin‑associated ferredoxin/periplasmic endochitinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,467,1200AG98.3% 396.4 / NA 117intergenic (‑129/+40)bfd/chiAbacterioferritin‑associated ferredoxin/periplasmic endochitinase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base G (54/60);  minor base T (1/1);  total (55/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

ATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTG  >  NC_000913/3467023‑3467213
                                                                                                 |                                                                                             
atGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAGATATGAATGAAAATAATGGGCTACCAAAAGTAGCCCATTGGCa                                                                                             >  1:1580117/1‑100 (MQ=255)
atGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTGTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCa                                                                                             >  2:1542146/1‑100 (MQ=255)
   aaaaTAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCaaaa                                                                                          >  1:16139/1‑100 (MQ=255)
      aTAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAAt                                                                                       <  2:508468/100‑1 (MQ=255)
       tAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATg                                                                                      <  1:1028373/100‑1 (MQ=255)
        aGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATgc                                                                                     <  2:1461805/100‑1 (MQ=255)
        aGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATgc                                                                                     <  2:635767/100‑1 (MQ=255)
        aGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATgc                                                                                     >  1:2192870/1‑100 (MQ=255)
         gTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATgcg                                                                                    <  1:2344320/100‑1 (MQ=255)
         gTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATgcg                                                                                    >  2:720952/1‑100 (MQ=255)
           tCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATgcggc                                                                                  <  1:1208728/100‑1 (MQ=255)
           tCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACAAAAAGTAGCCCATTGGCaaaaaa                                                                                        >  2:1701239/1‑94 (MQ=255)
               aTTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATa                                                                              >  2:753758/1‑100 (MQ=255)
                tttCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATAc                                                                             >  1:1720250/1‑100 (MQ=255)
                 ttCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACt                                                                            <  2:2300739/100‑1 (MQ=255)
                   cAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTgg                                                                          <  2:1000477/100‑1 (MQ=255)
                    aaTTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGa                                                                         <  2:1870930/100‑1 (MQ=255)
                      ttACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAg                                                                       >  2:2055301/1‑100 (MQ=255)
                      ttACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAg                                                                       <  1:2022114/100‑1 (MQ=255)
                      ttACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATAAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAg                                                                       <  2:1407826/100‑1 (MQ=255)
                       tACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAgg                                                                      <  1:1930258/100‑1 (MQ=255)
                       tACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAgg                                                                      <  1:1933646/100‑1 (MQ=255)
                       tACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAgg                                                                      <  1:1425436/100‑1 (MQ=255)
                       tACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAgg                                                                      >  2:1359632/1‑100 (MQ=255)
                       tACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAgg                                                                      >  2:2345327/1‑100 (MQ=255)
                        actgAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATAccttcca                                                                       <  2:2727871/95‑7 (MQ=255)
                        aCGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGt                                                                     >  1:1772689/1‑100 (MQ=255)
                           gAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAggaga                                                                   >  1:603735/1‑96 (MQ=255)
                           gAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATc                                                                  >  1:2292145/1‑100 (MQ=255)
                           gAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATc                                                                  >  2:244383/1‑100 (MQ=255)
                           gAAAGACTCTTTTTTAAGCAAAGAATAGAAATAGGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATc                                                                  <  2:1196747/100‑1 (MQ=255)
                            aaaGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCg                                                                 <  1:702035/100‑1 (MQ=255)
                              aGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGTCAAAAAATGCGGCGATACTGGAAGGTATCGcc                                                               >  2:112791/1‑100 (MQ=255)
                               gACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGcc                                                               >  1:1490502/1‑99 (MQ=255)
                               gACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCa                                                              <  2:843117/100‑1 (MQ=255)
                               gACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCa                                                              <  1:771878/100‑1 (MQ=255)
                               gACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCa                                                              <  1:2283946/100‑1 (MQ=255)
                               gACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCa                                                              >  1:39245/1‑100 (MQ=255)
                                aCTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCaa                                                             >  1:2847368/1‑100 (MQ=255)
                                 ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                            <  1:384757/100‑1 (MQ=255)
                                 ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                            >  2:2406509/1‑100 (MQ=255)
                                 ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                            <  2:707792/100‑1 (MQ=255)
                                 ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                            <  2:2658951/100‑1 (MQ=255)
                                 ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                            <  2:1251553/100‑1 (MQ=255)
                                 ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                            >  1:1374079/1‑100 (MQ=255)
                                 ctctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                            >  1:1251557/1‑100 (MQ=255)
                                 ctctTTTGTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAac                                                            >  1:1524003/1‑100 (MQ=255)
                                  tctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAaca                                                           >  2:2547724/1‑100 (MQ=255)
                                  tctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAaca                                                           >  2:2547919/1‑100 (MQ=255)
                                  tctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAaca                                                           >  1:1789958/1‑100 (MQ=255)
                                   ctTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAacac                                                          >  1:1486525/1‑100 (MQ=255)
                                      ttttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACAtt                                                       <  1:1708364/100‑1 (MQ=255)
                                      ttttAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACAtt                                                       <  1:2603178/100‑1 (MQ=255)
                                          aaGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTAtt                                                   >  2:2071905/1‑100 (MQ=255)
                                           aGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTg                                                  <  2:368334/100‑1 (MQ=255)
                                           aGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTg                                                  >  1:1131750/1‑100 (MQ=255)
                                             cAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCt                                                >  1:838056/1‑100 (MQ=255)
                                             cAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCt                                                <  1:667336/100‑1 (MQ=255)
                                              aaaGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCtt                                               >  2:2294940/1‑100 (MQ=255)
                                              aaaGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCtt                                               <  2:467717/100‑1 (MQ=255)
                                              aaaGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCtt                                               >  2:2041625/1‑100 (MQ=255)
                                              aaaGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCtt                                               <  2:620924/100‑1 (MQ=255)
                                                aGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAg                                             <  1:2313192/100‑1 (MQ=255)
                                                aGAATAGAAATAGGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAg                                             <  1:531729/100‑1 (MQ=255)
                                                 gAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGt                                            <  1:1642121/100‑1 (MQ=255)
                                                   aTAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTaa                                          <  2:2787038/100‑1 (MQ=255)
                                                     aGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAAc                                        <  2:1131746/100‑1 (MQ=255)
                                                     aGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAAc                                        >  2:1129302/1‑100 (MQ=255)
                                                        aagatGAATAAAAATAATGGGCTACCATAAGTAGCCCATTGTCAAAAGATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACgg                                      <  1:2449666/96‑1 (MQ=255)
                                                         atatGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGgcg                                    <  2:1131094/100‑1 (MQ=255)
                                                          tatGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGgcgc                                   >  2:2464542/1‑100 (MQ=255)
                                                          tatGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGgcgc                                   <  1:2222847/100‑1 (MQ=255)
                                                             gAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGaa                                >  2:868549/1‑100 (MQ=255)
                                                             gAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGTGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGa                                 <  2:1594574/99‑1 (MQ=255)
                                                               aTAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGa                              >  1:2660953/1‑100 (MQ=255)
                                                                tAAAAATAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGAt                             >  1:1122261/1‑100 (MQ=255)
                                                                  aaaaTAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGAttt                           >  2:820527/1‑100 (MQ=255)
                                                                  aaaaTAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGAtt                            >  2:724395/1‑99 (MQ=255)
                                                                   aaaTAATGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTc                          >  2:756516/1‑100 (MQ=255)
                                                                     ataatGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGc                        <  2:2289807/100‑1 (MQ=255)
                                                                     ataatGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGc                        <  2:1772676/100‑1 (MQ=255)
                                                                       aatGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTa                      <  2:838055/100‑1 (MQ=255)
                                                                         tGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTaaa                    <  2:2292132/100‑1 (MQ=255)
                                                                         tGGGCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTaaa                    <  2:16139/100‑1 (MQ=255)
                                                                           ggCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATg                  >  2:1953751/1‑100 (MQ=255)
                                                                           ggCTACCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATg                  <  1:2411573/100‑1 (MQ=255)
                                                                               aCCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTa              >  2:2232815/1‑100 (MQ=255)
                                                                               aCCAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTa              >  2:2399661/1‑100 (MQ=255)
                                                                                ccAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTAc             >  1:37444/1‑100 (MQ=255)
                                                                                 cAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACt            <  1:2406525/100‑1 (MQ=255)
                                                                                 cAAAAGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACt            <  1:1700717/100‑1 (MQ=255)
                                                                                     aGTAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCgg        <  1:2077695/100‑1 (MQ=255)
                                                                                       tAGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTg      <  1:2477853/100‑1 (MQ=255)
                                                                                        aGCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGc     <  1:1147365/100‑1 (MQ=255)
                                                                                         gCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGGCAACACATTTATTGCTTAGTAAACGGCGCGAATATTTCGCTAAATGCGTACTTCGGTGcc    <  2:1328102/100‑1 (MQ=255)
                                                                                         gCCCATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGcc    <  2:1122257/100‑1 (MQ=255)
                                                                                          cccATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGcc    <  1:928740/99‑1 (MQ=255)
                                                                                          cccATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGACAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCt   <  2:1580108/100‑1 (MQ=255)
                                                                                           ccATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTg  <  1:1736788/100‑1 (MQ=255)
                                                                                           ccATTGGCAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTg  >  2:1619480/1‑100 (MQ=255)
                                                                                                 |                                                                                             
ATGAAAATAGTTCTTATTTCAATTACGGAAAGACTCTTTTTTAAGCAAAGAATAGAAATATGAATAAAAATAATGGGCTACCAAAAGTAGCCCATTGACAAAAAATGCGGCGATACTGGAAGGTATCGCCAACACATTTATTGCTTAGTAAACGGCGCGAAGATTTCGCTAAATGCGTACTTCGGTGCCTG  >  NC_000913/3467023‑3467213

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: ‑