Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,597,846 (T)5→4 100% intergenic (‑286/‑138) livK ← / → panM leucine transporter subunit/PanD autocleavage accelerator, panothenate synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,597,8420T.96.2% 485.2 / 4.3 113intergenic (‑282/‑142)livK/panMleucine transporter subunit/PanD autocleavage accelerator, panothenate synthesis
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (62/46);  total (65/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 120 total aligned reads displayed.

TATTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTAAAT  >  NC_000913/3597745‑3597937
                                                                                                 |                                                                                               
tATTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttt                                                                                               >  1:1935739/1‑100 (MQ=255)
 aTTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAtttt                                                                                              <  2:2076995/100‑1 (MQ=255)
  ttGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttttc                                                                                             >  2:2231702/1‑99 (MQ=255)
  ttGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttttc                                                                                             >  2:795091/1‑99 (MQ=255)
    gatagataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttttcat                                                                                           >  2:134552/1‑97 (MQ=255)
    gatagataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttttcat                                                                                           >  2:2480813/1‑97 (MQ=255)
    gatagataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACAttttcat                                                                                           <  1:2448241/100‑4 (MQ=255)
      tagataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtc                                                                                        <  1:826949/100‑1 (MQ=255)
        gataAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATtctc                                                                                      >  2:534571/1‑100 (MQ=255)
          taAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCATAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTa                                                                                    >  2:2309795/1‑100 (MQ=255)
          taAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTa                                                                                    <  2:1221862/100‑1 (MQ=255)
            aaCAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAAt                                                                                  >  2:535293/1‑100 (MQ=255)
            aaCAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAAt                                                                                  >  2:440706/1‑100 (MQ=255)
             aCAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATg                                                                                 <  2:2213248/100‑1 (MQ=255)
              cAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGt                                                                                >  1:1651193/1‑100 (MQ=255)
              cAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGt                                                                                >  1:1524001/1‑100 (MQ=255)
               aGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACACAGGGCGGAATAATTCACTATCATTCAGGGGCTTATGCTGGACATTTTCATTctcttctg                                                                                 >  2:2344629/1‑94 (MQ=255)
                gAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGttt                                                                              >  1:292487/1‑100 (MQ=255)
                gAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGttt                                                                              >  1:1143059/1‑100 (MQ=255)
                gAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGttt                                                                              >  1:1460212/1‑100 (MQ=255)
                gAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACGTTTTCATTCTCTAATGttt                                                                              <  1:646701/100‑1 (MQ=255)
                 aaaaaaGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAGTAATTCACAAACAATCAGGGGATTATGCTGGAAATTTTCAATCTCTAATGtttt                                                                             >  2:1825906/1‑100 (MQ=255)
                  aaaaaGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTa                                                                            >  2:268548/1‑100 (MQ=255)
                    aaaGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAAt                                                                          >  2:2261260/1‑100 (MQ=255)
                    aaaGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATACACTATCATACAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTGAt                                                                          >  2:2492398/1‑100 (MQ=255)
                     aaGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTaa                                                                           >  2:109589/1‑98 (MQ=255)
                      aGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAAttt                                                                        <  1:130957/100‑1 (MQ=255)
                         gCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGt                                                                     <  1:1802306/100‑1 (MQ=255)
                          ccTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTa                                                                    >  1:1732803/1‑100 (MQ=255)
                          ccTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTa                                                                    <  1:714812/100‑1 (MQ=255)
                          ccTTTGTCAACATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTa                                                                    >  1:952897/1‑100 (MQ=255)
                           cTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTaa                                                                   >  1:2495338/1‑100 (MQ=255)
                           cTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTaa                                                                   <  2:730162/100‑1 (MQ=255)
                             ttGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAatt                                                                 >  2:397356/1‑100 (MQ=255)
                              tctCAGCATAAACTAGCAGCATAATGGTCCGGATAATTCCCTAACATCAAGGCTATTATGTTAGACATTTTCATTCCCTAATTTTTTAATTTTGTAatta                                                                <  2:245510/98‑1 (MQ=255)
                              tGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAatta                                                                <  2:53042/100‑1 (MQ=255)
                              tGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAat                                                                  >  1:2462956/1‑98 (MQ=255)
                               gTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAattat                                                               >  2:1919333/1‑100 (MQ=255)
                               gTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAattat                                                               >  1:989176/1‑100 (MQ=255)
                                tCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAattatt                                                              <  1:2225335/100‑1 (MQ=255)
                                tCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAattatt                                                              <  1:1752222/100‑1 (MQ=255)
                                tCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAattatt                                                              <  2:240504/100‑1 (MQ=255)
                                 cAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAattatt                                                              >  2:544503/1‑99 (MQ=255)
                                 cAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTg                                                             >  1:577209/1‑100 (MQ=255)
                                 cAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTg                                                             >  1:512905/1‑100 (MQ=255)
                                   gCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTg                                                             >  1:612325/1‑98 (MQ=255)
                                      tAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGtt                                                        >  1:2282971/1‑100 (MQ=255)
                                      tAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGtt                                                        >  1:1621336/1‑100 (MQ=255)
                                      tAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGtt                                                        >  2:1441327/1‑100 (MQ=255)
                                      tAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGtt                                                        >  2:1431813/1‑100 (MQ=255)
                                       aaaaTAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTa                                                       >  1:2291716/1‑100 (MQ=255)
                                         aaTAACGGCACAAAGGGCGGAATAATTCACTGTCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTaaa                                                     >  1:1491938/1‑100 (MQ=255)
                                         aaTAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTaaa                                                     >  2:1295592/1‑100 (MQ=255)
                                         aaTAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTaa                                                      >  2:1464196/1‑99 (MQ=255)
                                            aaCGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTaaaaaa                                                  <  1:1683885/100‑1 (MQ=255)
                                              cGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAtt                                                <  2:1570399/100‑1 (MQ=255)
                                              cGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAAtt                                                <  1:1578528/100‑1 (MQ=255)
                                               ggCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTa                                               >  2:1936977/1‑100 (MQ=255)
                                                gCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTaa                                              >  1:1667673/1‑100 (MQ=255)
                                                 cacaAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAAt                                             >  1:1928217/1‑100 (MQ=255)
                                                  acaAAGGGCGGAATAATTTACTATCATTCAGTTTATTATGTTGGATATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATc                                            <  2:1752539/100‑1 (MQ=255)
                                                  acaAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATc                                            <  2:544401/100‑1 (MQ=255)
                                                     aaGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTaaagatc                                         >  2:1610588/1‑95 (MQ=255)
                                                          cGGAATAATTCACTGTCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCa                                    <  2:1491927/100‑1 (MQ=255)
                                                          cGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCa                                    >  1:1223019/1‑100 (MQ=255)
                                                            gAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTTATTATTGCTGTTAAAAAATTAATCACCTGCCaaa                                  >  2:1286961/1‑100 (MQ=255)
                                                              ataatTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAg                                <  1:1246093/100‑1 (MQ=255)
                                                               taatTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGa                               >  1:986287/1‑100 (MQ=255)
                                                               taatTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGa                               <  1:1550286/100‑1 (MQ=255)
                                                                 atTCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGaaa                             >  2:429807/1‑100 (MQ=255)
                                                                  ttCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAAt                            >  2:960420/1‑100 (MQ=255)
                                                                  ttCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAAt                            >  2:1752132/1‑100 (MQ=255)
                                                                   tCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTACAAAATTAATCACCTGCCAAAAGaaa                             >  2:1828109/1‑98 (MQ=255)
                                                                   tCACTATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATa                           <  1:709921/100‑1 (MQ=255)
                                                                       tATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaaaaa                       <  1:1712217/100‑1 (MQ=255)
                                                                       tATCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATaaaaa                       <  1:661168/100‑1 (MQ=255)
                                                                         tCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAag                     >  2:717058/1‑100 (MQ=255)
                                                                         tCATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAag                     >  1:1921666/1‑100 (MQ=255)
                                                                          cATTCAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAaga                    >  2:1432525/1‑100 (MQ=255)
                                                                           attggggggATTGTGCTGGACGTTTTAATTCGGTGATGTTTTGAGTTTGGAATTATTGCAGTTACAAAGTTAATCACCTGCCACAAGAATTAAAAAag                     >  2:156977/6‑98 (MQ=255)
                                                                           aTTCAGGGGATTATGCTGGACATTTTCATTATCTAATGTTTTAATTTTGTAATTATTGCTGTTATAAAATTAATCACCTGCCAAAAGAAGTAAAAAaggg                   >  1:2020344/1‑98 (MQ=255)
                                                                            ttCAGGGGATTATGCTTGATATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAagag                   <  2:1868763/99‑1 (MQ=255)
                                                                              cAGGGGATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAgaga                  >  2:1579041/1‑98 (MQ=255)
                                                                               agtaaaTTATGCTGGATATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTTATGTTAAAAAATTAATCACCTGTCAAAATAAATAAAAAAGAGAAAg               <  2:1598216/95‑1 (MQ=255)
                                                                               aGGGGATTATGCTGGACATTGTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAg               >  2:2021452/1‑100 (MQ=255)
                                                                                ggggATTATGCTGGGTATTTTTCATTCTATAATGTTTTAACTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAg               <  2:301307/100‑1 (MQ=255)
                                                                                ggggATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGc              >  2:1761831/1‑100 (MQ=255)
                                                                                 gggATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGcc             >  2:1273492/1‑100 (MQ=255)
                                                                                  ggATTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCt            >  1:2450133/1‑100 (MQ=255)
                                                                                    aTTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTcc          <  2:2147687/100‑1 (MQ=255)
                                                                                    aTTATGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTcc          <  1:121704/100‑1 (MQ=255)
                                                                                    aTTATGCTGGACATGTTCATTCTCTAATTTTTTAATTTTGTAAGAAAAAAAAAAAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTcc          <  2:818818/100‑1 (MQ=255)
                                                                                       aTGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGAt       >  1:1616963/1‑100 (MQ=255)
                                                                                        tGCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGAtt      >  1:1743509/1‑100 (MQ=255)
                                                                                         gCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCATCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGAtt      >  2:1385937/1‑99 (MQ=255)
                                                                                         gCTGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTa     <  1:2374447/100‑1 (MQ=255)
                                                                                           tGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTaaa   <  1:2231716/100‑1 (MQ=255)
                                                                                           tGGACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTaaa   <  1:2033635/100‑1 (MQ=255)
                                                                                            ggACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTAAat  <  1:489995/100‑1 (MQ=255)
                                                                                            ggACATTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTAAat  <  1:440710/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TATTGATAGATAAACAGAAAAAAGTGCCTTTGTCAGCATAAAATAACGGCACAAAGGGCGGAATAATTCACTATCATTCAGGGGATTATGCTGGACATTTTTCATTCTCTAATGTTTTAATTTTGTAATTATTGCTGTTAAAAAATTAATCACCTGCCAAAAGAAATAAAAAAGAGAAAGCCTCCGATTAAAT  >  NC_000913/3597745‑3597937

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: