Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,406,377 A→C 9.7% C89G (TGC→GGC) ‡ lrhA ← DNA‑binding transcriptional dual regulator LrhA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,406,3770AC 9.7% 202.2 / 14.2 72C89G (TGC→GGC) ‡lrhADNA‑binding transcriptional dual regulator LrhA
Reads supporting (aligned to +/- strand):  ref base A (28/37);  new base C (5/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

TTACGTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGA  >  NC_000913/2406232‑2406524
                                                                                                                                                 |                                                                                                                                                   
ttACGTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGcc                                                                                                                                                 <  1:269143/150‑1 (MQ=255)
 tACGTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCt                                                                                                                                                >  2:974808/1‑150 (MQ=255)
   cGTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCtca                                                                                                                                              >  2:1275408/1‑150 (MQ=255)
   cGTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCtca                                                                                                                                              <  1:906657/150‑1 (MQ=255)
        aaCGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCAtt                                                                                                                                         >  2:622456/1‑150 (MQ=255)
         aCGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTa                                                                                                                                        <  1:752195/150‑1 (MQ=255)
           gcgGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGGCCTGCTCATCATTaat                                                                                                                                      <  2:1081435/150‑2 (MQ=255)
           gcgGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGGCCTGCTCATCATTaat                                                                                                                                      <  1:168224/150‑2 (MQ=255)
           gcgGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTaaa                                                                                                                                      <  1:844232/150‑1 (MQ=255)
                cATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCa                                                                                                                                 >  1:569941/1‑150 (MQ=255)
                 aTCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAg                                                                                                                                <  2:1468114/150‑1 (MQ=255)
                  tCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAg                                                                                                                                <  2:991566/149‑1 (MQ=255)
                  tCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAg                                                                                                                                >  1:991566/1‑149 (MQ=255)
                    cAGCGCTAATTTCGGATAAACCTAACTCACGCGATTTAACAGGAAACGTAAGATCGTATCTGCAGAATAATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGAt                                                                                                                             <  1:447619/150‑1 (MQ=255)
                           aaTTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTg                                                                                                                      <  2:1049791/150‑1 (MQ=255)
                             tttCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGc                                                                                                                    <  1:1031442/150‑1 (MQ=255)
                              ttCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCg                                                                                                                   <  2:1599983/150‑1 (MQ=255)
                               tCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGt                                                                                                                  >  1:740883/1‑150 (MQ=255)
                                 ggATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTaa                                                                                                                >  2:1048914/1‑150 (MQ=255)
                                          gAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGt                                                                                                       >  1:1283736/1‑150 (MQ=255)
                                          gAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGt                                                                                                       <  2:608345/150‑1 (MQ=255)
                                               cACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAAt                                                                                                  >  2:641803/1‑150 (MQ=255)
                                               cACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAAt                                                                                                  <  1:236355/150‑1 (MQ=255)
                                                aCGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATg                                                                                                 <  1:1048914/150‑1 (MQ=255)
                                                   cgATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCa                                                                                              <  2:260857/150‑1 (MQ=255)
                                                        tAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTc                                                                                         <  1:641803/150‑1 (MQ=255)
                                                           cAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGt                                                                                      <  1:539578/150‑1 (MQ=255)
                                                               aaaGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAAc                                                                                  <  2:569941/150‑1 (MQ=255)
                                                                aaGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACa                                                                                 >  2:838330/1‑150 (MQ=255)
                                                                  ggTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGt                                                                               >  2:429360/1‑150 (MQ=255)
                                                                    tAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGttt                                                                             >  1:953030/1‑150 (MQ=255)
                                                                      aGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTttgt                                                                           <  2:740883/150‑1 (MQ=255)
                                                                         tCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGc                                                                        <  2:8384/150‑1 (MQ=255)
                                                                          cGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCg                                                                       >  2:820112/1‑150 (MQ=255)
                                                                              tCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCg                                                                   >  2:444107/1‑150 (MQ=255)
                                                                               cGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGt                                                                  <  1:1405942/150‑1 (MQ=255)
                                                                                ggCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTg                                                                 >  2:704780/1‑150 (MQ=255)
                                                                                   aGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACg                                                              <  2:612844/150‑1 (MQ=255)
                                                                                   aGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACg                                                              >  1:364559/1‑150 (MQ=255)
                                                                                      ttCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGc                                                           <  1:409517/150‑1 (MQ=255)
                                                                                      ttCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGc                                                           <  2:1508157/150‑1 (MQ=255)
                                                                                       tCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCg                                                          <  2:469583/150‑1 (MQ=255)
                                                                                        cATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGATGAATGCCATGTTCAGTTAACAGATTGTTGCGACCGTGACGAGCGa                                                         >  2:1572209/1‑150 (MQ=255)
                                                                                         aTCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGaa                                                        >  2:16558/1‑150 (MQ=255)
                                                                                          tCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAAc                                                       <  1:444107/150‑1 (MQ=255)
                                                                                               aGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTc                                                  >  2:143118/1‑150 (MQ=255)
                                                                                                  aCCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCttt                                               <  1:1040088/150‑1 (MQ=255)
                                                                                                  aCCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCttt                                               <  2:1470905/150‑1 (MQ=255)
                                                                                                    cGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTcc                                             <  1:256657/150‑1 (MQ=255)
                                                                                                       tAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCaa                                          >  2:1567295/1‑150 (MQ=255)
                                                                                                           tAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAg                                      >  2:12479/1‑150 (MQ=255)
                                                                                                           tAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAg                                      >  1:1481589/1‑150 (MQ=255)
                                                                                                              tACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAgttg                                   >  2:1113600/1‑150 (MQ=255)
                                                                                                                  cccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTcc                               <  1:1355727/150‑1 (MQ=255)
                                                                                                                   ccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCa                              >  2:952695/1‑150 (MQ=255)
                                                                                                                   ccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCa                              <  1:622456/150‑1 (MQ=255)
                                                                                                                       aaGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTGAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACg                          <  2:706254/150‑1 (MQ=255)
                                                                                                                            tACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCa                     <  2:415781/150‑1 (MQ=255)
                                                                                                                            tACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCa                     <  1:820112/150‑1 (MQ=255)
                                                                                                                            tACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCa                     <  1:794618/150‑1 (MQ=255)
                                                                                                                                 cgttttAATGATGAGGCCTGCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGCGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGCCGCTGCATTTGc                >  2:1430153/5‑150 (MQ=255)
                                                                                                                                 cgttttAATGATGAGGCCTGCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGc                >  2:1435874/5‑150 (MQ=255)
                                                                                                                                  gttttAATGATGAGGCCTGCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCt               >  1:1566794/4‑150 (MQ=255)
                                                                                                                                  gttttAATGATGAGGCCTGCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCt               >  1:1060637/4‑150 (MQ=255)
                                                                                                                                     ttAATGATGAGGCCTGCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGAc            >  2:38716/1‑150 (MQ=255)
                                                                                                                                       aaTGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACtt          <  1:704780/150‑1 (MQ=255)
                                                                                                                                        atgatgAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTa         >  2:718385/1‑150 (MQ=255)
                                                                                                                                         tgatgaGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTAc        >  2:1132515/1‑150 (MQ=255)
                                                                                                                                            tgaGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTAcggc     <  1:838330/150‑1 (MQ=255)
                                                                                                                                            tgaGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTAcggc     <  2:953030/150‑1 (MQ=255)
                                                                                                                                              aGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTAcggcgg   >  1:698320/1‑150 (MQ=255)
                                                                                                                                               gCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGa  >  1:307293/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTACGTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGA  >  NC_000913/2406232‑2406524

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