Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,315,705 T→C 49.6% 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,7050TC49.6% ‑2.9 / 194.9 139K22K (AAA→AAGyciFputative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (40/30);  new base C (23/46);  total (63/76)
Fisher's exact test for biased strand distribution p-value = 6.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 144 total aligned reads displayed.

AATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCA  >  NC_000913/1315606‑1315800
                                                                                                   |                                                                                               
  tCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGctt                                                                                               <  2:2828224/100‑1 (MQ=255)
   cTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGtttt                                                                                              <  1:3048646/100‑1 (MQ=255)
    tGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGttttt                                                                                             >  1:384340/1‑100 (MQ=255)
      tCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCt                                                                                           <  1:1148998/100‑1 (MQ=255)
        cATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTgc                                                                                         >  1:941810/1‑100 (MQ=255)
        cATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTgc                                                                                         >  1:1663397/1‑100 (MQ=255)
          tGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTgcgc                                                                                       >  2:4226224/1‑100 (MQ=255)
           gAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCt                                                                                      <  2:2030592/100‑1 (MQ=255)
           gAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTGAATTTTTCATTTGATGTTGCTCTTGCGAGTTCTGCCTGTGCCCAGGTTAAGTGCTTtgctgccc                                                                                       <  1:1952946/99‑3 (MQ=255)
            aGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTg                                                                                     <  2:3395694/100‑1 (MQ=255)
             gTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTTTTTTGTTTTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTTATTGTTTTTCTTCGCTGt                                                                                    >  1:1285263/1‑100 (MQ=255)
             gTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGt                                                                                    >  1:960327/1‑100 (MQ=255)
             gTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGt                                                                                    >  2:3411470/1‑100 (MQ=255)
             gTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGt                                                                                    <  2:1111103/100‑1 (MQ=255)
             gTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGt                                                                                    >  1:3121375/1‑100 (MQ=255)
              tttCCTCGAGGTGCGCATGCAAAGCCTTACTTAACTTTTCATTTGATGTTGCGCTTTCGAGTTTTGCCAGTGCCCGCGTTAATTGTTTTTCTGCGCTGTa                                                                                   <  1:2318507/100‑1 (MQ=255)
              tttCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTa                                                                                   <  1:2992908/100‑1 (MQ=255)
               ttCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAg                                                                                  <  1:3994242/100‑1 (MQ=255)
                    cGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATc                                                                             >  1:3960927/1‑100 (MQ=255)
                    agaaaaaaaaaTGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATc                                                                             <  1:2746834/90‑1 (MQ=255)
                      aGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTTCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTg                                                                           <  1:821399/100‑1 (MQ=255)
                      aGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTg                                                                           >  2:295680/1‑100 (MQ=255)
                       ggTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAgagatcgga                                                                          >  1:1489512/1‑92 (MQ=255)
                         tGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGaaa                                                                        >  1:1479011/1‑100 (MQ=255)
                         tGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGaaa                                                                        <  2:420209/100‑1 (MQ=255)
                          gcgcATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAg                                                                       <  2:4081870/100‑1 (MQ=255)
                              aTGAAAAGCCTGACTTAGTTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAgg                                                                   <  1:499883/100‑1 (MQ=255)
                              aTGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAgg                                                                   <  2:405379/100‑1 (MQ=255)
                               tGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGt                                                                  <  2:2332755/100‑1 (MQ=255)
                                gAAAAGCCTGACTTAATTTTTCTTTTTATGTTGCTATTGCGAGTTTTGCCAATGCACGGGTTAATTGTTTTTCTACGCTGTAGGTATCTGAAAgaa                                                                     >  1:861821/1‑94 (MQ=255)
                                gAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTg                                                                 <  2:2044552/100‑1 (MQ=255)
                                gAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTg                                                                 <  2:2499729/100‑1 (MQ=255)
                                gAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTg                                                                 <  1:1041546/100‑1 (MQ=255)
                                gAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTg                                                                 <  2:465036/100‑1 (MQ=255)
                                 aaaaGACTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGTCCTGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAtgtg                                                                 <  1:658876/99‑1 (MQ=255)
                                    aGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaata                                                             <  2:3046992/100‑1 (MQ=255)
                                    aGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaata                                                             <  1:934486/100‑1 (MQ=255)
                                     gCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataa                                                            >  2:3901251/1‑100 (MQ=255)
                                      ccTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaa                                                           <  2:1462883/100‑1 (MQ=255)
                                      ccTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaa                                                           >  2:2623126/1‑100 (MQ=255)
                                      ccTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAAGTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaa                                                           <  1:3700453/100‑1 (MQ=255)
                                       cTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaat                                                          >  2:1277457/1‑100 (MQ=255)
                                        tGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGaataaata                                                         <  2:3649418/100‑1 (MQ=255)
                                         gACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATAc                                                        <  1:1755064/100‑1 (MQ=255)
                                           nttAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACAt                                                      >  1:990693/2‑100 (MQ=255)
                                           cTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACAt                                                      >  1:3714402/1‑100 (MQ=255)
                                           cTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACAt                                                      >  2:1684779/1‑100 (MQ=255)
                                            ttAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATc                                                     >  1:1843147/1‑100 (MQ=255)
                                                tttttCATTTGATGTTGCTCTTGCGAGCTTTGCCAGTGCCCGTGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCa                                                 >  1:639905/1‑100 (MQ=255)
                                                   ttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCGCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATg                                              >  2:3004434/1‑100 (MQ=255)
                                                   ttCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATg                                              <  1:4032989/100‑1 (MQ=255)
                                                    tCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATgg                                             >  1:265235/1‑100 (MQ=255)
                                                    tCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATgg                                             <  2:2300450/100‑1 (MQ=255)
                                                      aTTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTc                                           >  1:3873458/1‑100 (MQ=255)
                                                         tGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttc                                        >  2:2965574/1‑100 (MQ=255)
                                                           aTGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcat                                      <  2:3334764/100‑1 (MQ=255)
                                                            tGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcata                                     <  1:3120337/100‑1 (MQ=255)
                                                             gTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatat                                    <  2:2342/100‑1 (MQ=255)
                                                             gTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatat                                    >  2:1524378/1‑100 (MQ=255)
                                                             gTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatat                                    >  1:424873/1‑100 (MQ=255)
                                                             gTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatat                                    <  2:387482/100‑1 (MQ=255)
                                                               tGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattc                                  >  1:2381232/1‑100 (MQ=255)
                                                               tGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattc                                  >  2:1442311/1‑100 (MQ=255)
                                                               tGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattc                                  <  1:919789/100‑1 (MQ=255)
                                                                gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattca                                 >  1:137264/1‑100 (MQ=255)
                                                                gCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattca                                 <  1:2278254/100‑1 (MQ=255)
                                                                 ctctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGCAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattcat                                >  1:3005032/1‑100 (MQ=255)
                                                                 ctctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattcat                                >  2:1410791/1‑100 (MQ=255)
                                                                  tctTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGCATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCttcatattcata                               >  2:2072506/1‑100 (MQ=255)
                                                                     tGCGAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATAttt                            >  1:1152829/1‑100 (MQ=255)
                                                                      gCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATAtttt                           >  2:774799/1‑100 (MQ=255)
                                                                      gCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATAtttt                           <  2:2758346/100‑1 (MQ=255)
                                                                        gAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtc                         <  1:3406204/100‑1 (MQ=255)
                                                                        gAGTTTTGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtc                         <  2:1560375/100‑1 (MQ=255)
                                                                          gTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTtctc                       >  1:647848/1‑100 (MQ=255)
                                                                            tttGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCa                     >  2:280891/1‑100 (MQ=255)
                                                                             ttttcAGTGCCCGGGTTAATTGTTTTTCTACTATGTAGGTATCTTAAATCAGGTTTATAAATACATCTTCCATGGTCTTCATATTCATATTTTTCTCCAg                    <  1:2076975/96‑1 (MQ=255)
                                                                             ttGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAg                    <  1:3407505/100‑1 (MQ=255)
                                                                             ttGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAg                    <  2:1974285/100‑1 (MQ=255)
                                                                              tGCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGt                   <  2:3580000/100‑1 (MQ=255)
                                                                               gCCTGTGCTCGGTTTAATTGTTTTTCTGCGCTTTAAGTATCTCAAAGCAGCTGAATAAATACATCTTCAACGGTCTTCATATTCATATTTTTCTCCAGTg                  <  1:679768/100‑1 (MQ=255)
                                                                               gCCAGTGCCCGGGTTAATTGCTTTTCTGCTCTGTAGGTATCTGAAAGCATGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTg                  <  1:972440/100‑1 (MQ=255)
                                                                               gCCAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTg                  <  2:3214499/100‑1 (MQ=255)
                                                                                ccAGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGa                 <  2:4053637/100‑1 (MQ=255)
                                                                                  aGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGaaa               <  1:819018/100‑1 (MQ=255)
                                                                                  aGTGCCCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGaaa               <  1:4127306/100‑1 (MQ=255)
                                                                                    tGCTCGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTCAAAGTAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATc             <  1:472322/100‑1 (MQ=255)
                                                                                     gCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTATTCATATTCATATTTTTCTCCAGTGAAATCa            >  2:1732421/1‑100 (MQ=255)
                                                                                      cccGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCAc           >  2:3210962/1‑100 (MQ=255)
                                                                                      cccGGGTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCAc           >  1:2401117/1‑100 (MQ=255)
                                                                                        cGGGTTAATTGTTTTTCTTCGCTGTATGTATCTGACAGCAGGTGAATAAATACATCTTCAATCGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTg         <  1:2607328/100‑1 (MQ=255)
                                                                                        cGGGTTAATTGCTTTTCTGCGCTGTAGGTAACTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGAGAAATCACTg         >  1:2814944/1‑100 (MQ=255)
                                                                                          ggTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCg       <  2:229888/100‑1 (MQ=255)
                                                                                          ggTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCg       <  1:1742537/100‑1 (MQ=255)
                                                                                           gTTAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAAGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGcta      >  1:2139246/1‑98 (MQ=255)
                                                                                            ttAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAg     <  2:426859/100‑1 (MQ=255)
                                                                                            ttAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAg     <  2:941813/100‑1 (MQ=255)
                                                                                             tAATTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGc    <  2:1913586/100‑1 (MQ=255)
                                                                                              aaTTGCTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGcc   >  1:2440117/1‑100 (MQ=255)
                                                                                               aTTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCa  >  1:683884/1‑100 (MQ=255)
                                                                                                   |                                                                                               
AATCTGTCCATGAGTTTCCTCGAGGTGCGCATGAAAAGCCTGACTTAATTTTTCATTTGATGTTGCTCTTGCGAGTTTTGCCAGTGCCCGGGTTAATTGTTTTTCTGCGCTGTAGGTATCTGAAAGCAGGTGAATAAATACATCTTCAATGGTCTTCATATTCATATTTTTCTCCAGTGAAATCACTGCGAGCCA  >  NC_000913/1315606‑1315800

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: