Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,089 A→G 5.8% intergenic (‑47/+12) alsB ← / ← alsR D‑allose ABC transporter periplasmic binding protein/d‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,0890AG 5.8% 337.2 / 7.7 104intergenic (‑47/+12)alsB/alsRD‑allose ABC transporter periplasmic binding protein/d‑allose‑inducible als operon transcriptional repressor; autorepressor; repressor of rpiR
Reads supporting (aligned to +/- strand):  ref base A (49/49);  new base G (4/2);  total (53/51)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

Only 100 of 107 total aligned reads displayed.

GAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGC  >  NC_000913/4311941‑4312235
                                                                                                                                                    |                                                                                                                                                  
gagGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGaa                                                                                                                                                    >  2:146007/1‑149 (MQ=255)
gagGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGaa                                                                                                                                                    >  4:371017/1‑149 (MQ=255)
gagGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGaa                                                                                                                                                    >  3:53050/1‑149 (MQ=255)
    gTTTTCAATACGGCAGCATATTCGGAGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                                                                                >  3:65208/1‑149 (MQ=255)
    gTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                                                                                <  3:508324/149‑1 (MQ=255)
    gTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                                                                                >  4:538592/1‑149 (MQ=255)
     ttttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                                                                               >  3:494053/1‑149 (MQ=255)
      ttttAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCt                                                                                                                                              >  3:68728/1‑149 (MQ=255)
      tttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCt                                                                                                                                              <  2:826612/149‑1 (MQ=255)
          aaTACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAc                                                                                                                                                >  3:272389/1‑143 (MQ=255)
          aaTACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTAc                                                                                                                                                <  4:272389/143‑1 (MQ=255)
          aaTACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                                                                                          >  1:741779/1‑149 (MQ=255)
          aaTACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                                                                                          >  1:143307/1‑149 (MQ=255)
            tACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCt                                                                                                                                        >  2:594273/1‑149 (MQ=255)
              cGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGCGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                                                                                      <  2:537082/149‑1 (MQ=255)
                   gCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACGTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCa                                                                                                                                 >  2:800175/1‑149 (MQ=255)
                     atatTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgc                                                                                                                               <  3:201623/149‑1 (MQ=255)
                      tatTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcg                                                                                                                              >  3:26724/1‑149 (MQ=255)
                        ttCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcg                                                                                                                            >  1:484643/1‑149 (MQ=255)
                          cggcggCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGcc                                                                                                                          <  3:699633/149‑1 (MQ=255)
                               gcagcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                                                                                     >  2:875580/1‑149 (MQ=255)
                                cagcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTga                                                                                                                    <  3:92307/149‑1 (MQ=255)
                                cagcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTga                                                                                                                    <  1:579526/149‑1 (MQ=255)
                                 agcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgag                                                                                                                   <  1:594273/149‑1 (MQ=255)
                                  gcaAAACCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTATCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACTTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                                                                                  >  2:862791/1‑149 (MQ=255)
                                   caAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaa                                                                                                                 <  1:71344/149‑1 (MQ=255)
                                    aaaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                                                                                >  1:669559/1‑149 (MQ=255)
                                    aaaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                                                                                >  2:110401/1‑149 (MQ=255)
                                    aaaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                                                                                >  3:394266/1‑149 (MQ=255)
                                      aaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTGCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAgaagaa                                                                                                              >  4:708251/1‑149 (MQ=255)
                                       aGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGaaa                                                                                                             <  4:171618/149‑1 (MQ=255)
                                        gcgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAt                                                                                                            <  4:189055/149‑1 (MQ=255)
                                        gcgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAt                                                                                                            >  3:702248/1‑149 (MQ=255)
                                           cTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                                                                                         <  2:128258/149‑1 (MQ=255)
                                           cTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                                                                                         <  2:202116/149‑1 (MQ=255)
                                           cTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                                                                                         <  4:615041/149‑1 (MQ=255)
                                            tGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaac                                                                                                        <  1:700299/149‑1 (MQ=255)
                                                tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGTCTGGTGAGAAGAAATCAACAAtt                                                                                                    <  4:172271/149‑1 (MQ=255)
                                                   caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcg                                                                                                 >  1:142986/1‑149 (MQ=255)
                                                    aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                                                                                <  1:494509/149‑1 (MQ=255)
                                                      caTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCg                                                                                              >  3:404717/1‑149 (MQ=255)
                                                       aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                                                                                             >  1:609183/1‑149 (MQ=255)
                                                           ggCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                                                                                >  4:750067/2‑126 (MQ=255)
                                                           aGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                                                                                <  3:750067/126‑1 (MQ=255)
                                                             cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                                                                                       <  2:321951/149‑1 (MQ=255)
                                                             cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                                                                                       <  1:402389/149‑1 (MQ=255)
                                                             cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                                                                                       <  2:669559/149‑1 (MQ=255)
                                                              ccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTg                                                                                      <  4:831038/149‑1 (MQ=255)
                                                                      gtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgc                                                                                                                           >  4:429423/1‑104 (MQ=255)
                                                                      gtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgcgcgc                                                                                                                           <  3:429423/104‑1 (MQ=255)
                                                                      gtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                                                                              <  1:146007/149‑1 (MQ=255)
                                                                      gtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                                                                              <  4:676273/149‑1 (MQ=255)
                                                                       tgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaat                                                                             <  2:142986/149‑1 (MQ=255)
                                                                          cGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatt                                                                          <  4:65208/149‑1 (MQ=255)
                                                                          cGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatt                                                                          >  2:137250/1‑149 (MQ=255)
                                                                           gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatta                                                                         <  4:624194/149‑1 (MQ=255)
                                                                                  tatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTCCGCCGGTTTTTTGCATATTAATATTAGCTTGtt                                                                  >  4:760107/1‑149 (MQ=255)
                                                                                    tttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCa                                                                >  1:329921/1‑149 (MQ=255)
                                                                                           tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGa                                                         <  4:26724/149‑1 (MQ=255)
                                                                                             tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGc                                                       >  3:763151/1‑149 (MQ=255)
                                                                                             tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGc                                                       <  1:811522/149‑1 (MQ=255)
                                                                                             tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCt                                                      >  2:378430/1‑149 (MQ=255)
                                                                                                    ataataATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGa                                                >  3:95530/1‑149 (MQ=255)
                                                                                                     taataaTGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGAc                                               >  4:496219/1‑149 (MQ=255)
                                                                                                        taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACaga                                            <  3:708251/149‑1 (MQ=255)
                                                                                                           tGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGGTGGGCGACAgagaa                                         <  2:387489/149‑2 (MQ=255)
                                                                                                                tcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaa                                    >  1:397082/1‑149 (MQ=255)
                                                                                                                  aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaaaa                                  <  4:54135/149‑1 (MQ=255)
                                                                                                                   cGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAg                                 >  2:527116/1‑149 (MQ=255)
                                                                                                                       aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCa                                                                                    >  2:471282/1‑94 (MQ=255)
                                                                                                                       aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCa                                                                                    <  1:471282/94‑1 (MQ=255)
                                                                                                                       aaGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTc                             >  1:519969/1‑149 (MQ=255)
                                                                                                                        aGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGc                                >  4:164578/1‑145 (MQ=255)
                                                                                                                        aGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGc                                <  3:164578/145‑1 (MQ=255)
                                                                                                                         ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCa                           <  4:68728/149‑1 (MQ=255)
                                                                                                                         ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCa                           <  2:346065/149‑1 (MQ=255)
                                                                                                                          gaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAg                          >  1:37245/6‑149 (MQ=255)
                                                                                                                          gaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAg                          >  4:337668/6‑149 (MQ=255)
                                                                                                                             aCAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCaaa                      <  1:378430/149‑1 (MQ=255)
                                                                                                                                 gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTgg                                                    >  3:514956/1‑116 (MQ=255)
                                                                                                                                 gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTgg                                                    <  4:514956/116‑1 (MQ=255)
                                                                                                                                 gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGc                                <  2:790537/136‑1 (MQ=255)
                                                                                                                                 gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGc                                >  1:790537/1‑136 (MQ=255)
                                                                                                                                 gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagtt                   >  1:537395/1‑149 (MQ=255)
                                                                                                                                   aattCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGt                              >  4:815468/4‑136 (MQ=255)
                                                                                                                                   aattCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGt                              <  3:815468/133‑1 (MQ=255)
                                                                                                                                    aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGc                                                  >  1:759877/1‑115 (MQ=255)
                                                                                                                                    aCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGc                                                  <  2:759877/115‑1 (MQ=255)
                                                                                                                                      tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTg              >  3:760334/1‑149 (MQ=255)
                                                                                                                                      tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTg              >  3:362898/1‑149 (MQ=255)
                                                                                                                                      tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTg              >  1:280005/1‑149 (MQ=255)
                                                                                                                                         gcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTaa           <  4:362898/149‑1 (MQ=255)
                                                                                                                                           gATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAata         >  1:445394/1‑149 (MQ=255)
                                                                                                                                           gATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAata         >  1:585007/1‑149 (MQ=255)
                                                                                                                                           catcgCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAata         >  4:573316/5‑149 (MQ=255)
                                                                                                                                               gCGGAATTACATTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTGGCATATAAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCt     >  4:562480/1‑149 (MQ=255)
                                                                                                                                                cGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCtt    <  3:856284/149‑1 (MQ=255)
                                                                                                                                                 ggAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTg   >  1:364683/1‑149 (MQ=255)
                                                                                                                                                  gAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGc  <  1:800175/149‑1 (MQ=255)
                                                                                                                                                  gAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGc  <  1:514410/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGC  >  NC_000913/4311941‑4312235

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