Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,406,380 C→G 9.6% A88P (GCC→CCC) ‡ 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,3800CG 9.6% 223.6 / 12.7 73A88P (GCC→CCC) ‡lrhADNA‑binding transcriptional dual regulator LrhA
Reads supporting (aligned to +/- strand):  ref base C (29/37);  new base G (5/2);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

TTACGTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAG  >  NC_000913/2406232‑2406529
                                                                                                                                                    |                                                                                                                                                     
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                                                                                                            <  2:608345/150‑1 (MQ=255)
                                          gAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGt                                                                                                            >  1:1283736/1‑150 (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                                   <  1:622456/150‑1 (MQ=255)
                                                                                                                   ccTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCa                                   >  2:952695/1‑150 (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:1060637/4‑150 (MQ=255)
                                                                                                                                  gttttAATGATGAGGCCTGCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCt                    >  1:1566794/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)
                                                                                                                                                  ggCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTg    >  2:726024/1‑150 (MQ=255)
                                                                                                                                                  ggCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTg    <  1:1567295/150‑1 (MQ=255)
                                                                                                                                                    ccTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAg  >  1:851926/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTACGTTTAACGCGGACATCCAGCGCTAATTTCGGATAAACCGAACTCACGCGATTTAACAGGAAAGGTAAGATCGTATCGGCAGATTCATCTGAAGCACCGATAGTTAATACGCCCTGAAGATTACTGAACATTAATGATGAGCAGGCCTCATCATTAAAACGCAGGATTTTCCTGGCGTAACCAAGAAGTTGAATGCCATGTTCAGTTAACAGTTTGTTGCGACCGTGACGAGCGAACAGTTCTTTCCCAACGAGTTGTTCCAGACGCTGCATTTGCTGACTTACGGCGGACTGAG  >  NC_000913/2406232‑2406529

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: