Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,315,705 T→C K22K (AAA→AAG yciF ← putative rubrerythrin/ferritin‑like metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,315,7050TC98.9% 335.1 / ‑5.1 93K22K (AAA→AAGyciFputative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (44/48);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCAT  >  NC_000913/1315611‑1315801
                                                                                              |                                                                                                
gTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGcttttn                                                                                             <  2:1887627/100‑2 (MQ=255)
gTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTc                                                                                             >  1:2089181/1‑100 (MQ=255)
  ccATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTg                                                                                           <  2:2009622/100‑1 (MQ=255)
  ccATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTg                                                                                           >  2:286854/1‑100 (MQ=255)
  ccATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTg                                                                                           >  2:486463/1‑100 (MQ=255)
   cATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTgc                                                                                          >  1:248409/1‑100 (MQ=255)
     tGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTgcgc                                                                                        >  1:1042800/1‑100 (MQ=255)
     tGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTgcgc                                                                                        <  2:1884618/100‑1 (MQ=255)
     tGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTgcgc                                                                                        <  1:1684359/100‑1 (MQ=255)
        gTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGt                                                                                     >  2:2163067/1‑100 (MQ=255)
         tttCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTa                                                                                    <  1:722180/100‑1 (MQ=255)
          ttCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAg                                                                                   >  2:196871/1‑100 (MQ=255)
           tCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAgg                                                                                  <  1:2315287/100‑1 (MQ=255)
           tCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAgg                                                                                  <  2:2143874/100‑1 (MQ=255)
             cTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTa                                                                                >  2:201370/1‑100 (MQ=255)
             cTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTa                                                                                <  1:840712/100‑1 (MQ=255)
              tCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTAt                                                                               <  1:486465/100‑1 (MQ=255)
                 aGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTg                                                                            <  1:614046/100‑1 (MQ=255)
                  ggTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCtn                                                                            <  1:416615/99‑2 (MQ=255)
                  ggTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGa                                                                           <  1:870525/100‑1 (MQ=255)
                     gcgcCTGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAg                                                                        <  1:1120121/100‑1 (MQ=255)
                     gcgcATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGaaa                                                                         >  2:220587/1‑99 (MQ=255)
                     gcgcATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAg                                                                        <  1:1089502/100‑1 (MQ=255)
                     cgatctGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAg                                                                        <  1:1113968/95‑1 (MQ=255)
                       gcATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCa                                                                      <  1:2257119/100‑1 (MQ=255)
                        cATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAg                                                                     >  2:608342/1‑100 (MQ=255)
                         aTGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAgg                                                                    <  2:411808/100‑1 (MQ=255)
                         aTGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAgg                                                                    >  2:1266851/1‑100 (MQ=255)
                          tGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGt                                                                   >  1:1219572/1‑100 (MQ=255)
                           gAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAgagatcg                                                                  >  2:1113966/1‑94 (MQ=255)
                           gAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTg                                                                  <  1:100190/100‑1 (MQ=255)
                           gAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTg                                                                  <  1:1805635/100‑1 (MQ=255)
                           gAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTg                                                                  <  2:686028/100‑1 (MQ=255)
                           gAAAAGCATGACTTAATTTTTAATTTGATGTTGCTATTGCGAGTTTTGCCAGTGACCAGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGt                                                                   >  2:1612803/1‑99 (MQ=255)
                            aaaaGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGa                                                                 >  1:2154564/1‑100 (MQ=255)
                             aaaGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaa                                                                >  2:71461/1‑100 (MQ=255)
                             aaaGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaa                                                                <  2:1043573/100‑1 (MQ=255)
                             aaaGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaa                                                                >  2:2110314/1‑100 (MQ=255)
                             aaaGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaa                                                                <  1:705233/100‑1 (MQ=255)
                             aaaGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaa                                                                <  1:1965930/100‑1 (MQ=255)
                               aGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaatn                                                              <  1:1061452/100‑2 (MQ=255)
                                gCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataa                                                             <  1:2046413/100‑1 (MQ=255)
                                   tGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaata                                                          <  1:71463/100‑1 (MQ=255)
                                   tGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaata                                                          <  1:2110321/100‑1 (MQ=255)
                                          aTTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTc                                                   >  1:1937872/1‑100 (MQ=255)
                                           tttttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCa                                                  >  2:172869/1‑100 (MQ=255)
                                            ttttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCaa                                                 >  2:2317446/1‑100 (MQ=255)
                                                 cttccgatctGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTc                                            <  2:1507353/91‑1 (MQ=255)
                                                 aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTc                                            >  2:383147/1‑100 (MQ=255)
                                                 aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTc                                            >  2:1411540/1‑100 (MQ=255)
                                                   ttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCtt                                          <  1:481630/100‑1 (MQ=255)
                                                   ttGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCtt                                          >  2:2168460/1‑100 (MQ=255)
                                                    tGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttc                                         >  1:569539/1‑100 (MQ=255)
                                                       tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcata                                      >  2:1136673/1‑100 (MQ=255)
                                                       tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcata                                      >  2:459535/1‑100 (MQ=255)
                                                       tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGCTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcata                                      >  2:1430232/1‑100 (MQ=255)
                                                          tGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattc                                   <  2:1848514/100‑1 (MQ=255)
                                                          tGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattc                                   <  1:2168467/100‑1 (MQ=255)
                                                           gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTcagatcggaag                                  >  1:1507357/1‑90 (MQ=255)
                                                           gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattca                                  >  2:1793452/1‑100 (MQ=255)
                                                            ctctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattcat                                 >  2:2133088/1‑100 (MQ=255)
                                                              ctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTtcatattcatat                               >  1:816242/1‑100 (MQ=255)
                                                                tGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGATAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATAttt                             <  2:248407/100‑1 (MQ=255)
                                                                tGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATAttt                             >  1:153159/1‑100 (MQ=255)
                                                                 gCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATAtttt                            >  2:1610142/1‑100 (MQ=255)
                                                                 gCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGCATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATAtttt                            >  2:1991399/1‑100 (MQ=255)
                                                                  cGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTAAATGGTCTTCATATTCATAttttt                           >  1:1598547/1‑100 (MQ=255)
                                                                   tccgatctccAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtc                          <  2:1136910/92‑1 (MQ=255)
                                                                    aGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtct                         <  2:1042798/100‑1 (MQ=255)
                                                                     gTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtctc                        >  1:2037982/1‑100 (MQ=255)
                                                                      ttttGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTAGCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTcc                       <  2:1491739/100‑1 (MQ=255)
                                                                        ttGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAg                     <  1:1430236/100‑1 (MQ=255)
                                                                          gCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTg                   <  1:220589/100‑1 (MQ=255)
                                                                          gCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTg                   >  1:800106/1‑100 (MQ=255)
                                                                           ccAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtcagatcgga                  >  1:1136912/1‑92 (MQ=255)
                                                                                 cccGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCAc            <  1:459537/100‑1 (MQ=255)
                                                                                 cccGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCAc            >  1:27930/1‑100 (MQ=255)
                                                                                  ccGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCACTGGTCTTCATATTCATATTTTTCTCCAGTGAAATCAc            <  2:594332/99‑1 (MQ=255)
                                                                                    gggTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGc         <  2:2089174/100‑1 (MQ=255)
                                                                                     ggTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCg        >  2:1057962/1‑100 (MQ=255)
                                                                                     ggTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCg        <  2:153157/100‑1 (MQ=255)
                                                                                       ttAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAg      <  2:1878257/100‑1 (MQ=255)
                                                                                       ttAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAg      <  1:1136675/100‑1 (MQ=255)
                                                                                         taTTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGcc    <  1:434992/99‑1 (MQ=255)
                                                                                         taTTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGcc    <  1:1548745/99‑1 (MQ=255)
                                                                                         aaTTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGcc    <  1:2133095/100‑1 (MQ=255)
                                                                                         aaTTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGcc    <  2:27929/100‑1 (MQ=255)
                                                                                         aaTTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGcc    <  1:608343/100‑1 (MQ=255)
                                                                                          aTTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa   >  2:1548741/1‑100 (MQ=255)
                                                                                          aTTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa   >  2:434990/1‑100 (MQ=255)
                                                                                           ttGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCAt  >  1:1989887/1‑100 (MQ=255)
                                                                                           ttGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCAt  <  1:444073/100‑1 (MQ=255)
                                                                                           ttGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCAt  >  2:151585/1‑100 (MQ=255)
                                                                                           ttGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCAt  >  2:1380174/1‑100 (MQ=255)
                                                                                              |                                                                                                
GTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCAT  >  NC_000913/1315611‑1315801

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