Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,089 A→G 6.0% 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 6.0% 321.4 / 8.2 100intergenic (‑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 (42/52);  new base G (3/3);  total (45/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 102 total aligned reads displayed.

AGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGC  >  NC_000913/4311942‑4312235
                                                                                                                                                   |                                                                                                                                                  
agGGTTTTCAATACGACAGCATTTTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAACGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAt                                                                                                                                                   <  3:566718/149‑1 (MQ=255)
   gTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGaa                                                                                                                                                    >  4:239738/1‑145 (MQ=255)
   gTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGaa                                                                                                                                                    <  3:239738/145‑1 (MQ=255)
   gTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAc                                                                                                                                                <  4:513048/149‑1 (MQ=255)
   gTTTTCAATACAACAGCATATTCGGCGGCTGCAAAAGCGCTGGTTGACAACATTAAGCCGACGAGTGTGCCGCTAAAATATTTCAGATATTTATTCATAATAATGTTCTCACGTTAAGGGGACAGGTAACTGCTCGATGGCGGAATTAc                                                                                                                                                <  4:572056/149‑1 (MQ=255)
    ttttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                                                                               >  4:626881/1‑149 (MQ=255)
    ttttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                                                                               >  4:211720/1‑149 (MQ=255)
    ttttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                                                                               >  3:830175/1‑149 (MQ=255)
     tttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCt                                                                                                                                              <  2:530936/149‑1 (MQ=255)
     tttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCACGGAGTTACCt                                                                                                                                              <  1:489773/149‑1 (MQ=255)
      ttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCtt                                                                                                                                             <  3:626881/149‑1 (MQ=255)
         aaTACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                                                                                          <  2:256237/149‑1 (MQ=255)
          aTACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGTTGGCGGAATTACCTTTTTc                                                                                                                                         >  2:456219/1‑149 (MQ=255)
          aTACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                                                                         <  1:504621/149‑1 (MQ=255)
             cGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                                                                                      <  2:789038/149‑1 (MQ=255)
             cGACAGCATATTCGGCGGCAGCAAAAGCGCTGCTTGACAACCTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                                                                                      <  4:754259/149‑1 (MQ=255)
                cAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAttt                                                                                                                                   <  1:299062/149‑1 (MQ=255)
                 aGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTATCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGCATAACCTTTTTCTTATTTc                                                                                                                                  >  3:184179/1‑149 (MQ=255)
                  gCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGGAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCa                                                                                                                                 >  3:165398/1‑149 (MQ=255)
                   cATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAg                                                                                                                                >  3:507891/1‑149 (MQ=255)
                    atatTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgc                                                                                                                               <  3:402030/149‑1 (MQ=255)
                           gcggcAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTg                                                                                                                        <  3:118186/149‑1 (MQ=255)
                             ggcAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGt                                                                                                                      <  2:451145/149‑1 (MQ=255)
                              gcagcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTg                                                                                                                     >  1:211719/1‑149 (MQ=255)
                                  caAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagna                                                                                                                 <  3:180290/149‑1 (MQ=255)
                                   aaaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                                                                                >  2:540694/1‑149 (MQ=255)
                                   aaaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                                                                                >  2:693886/1‑149 (MQ=255)
                                       gcgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAt                                                                                                            >  2:931482/1‑149 (MQ=255)
                                           tGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaac                                                                                                        >  1:778801/1‑149 (MQ=255)
                                            ggTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaaca                                                                                                       >  3:415733/1‑149 (MQ=255)
                                               tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAtt                                                                                                    <  3:322645/149‑1 (MQ=255)
                                                 acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGGCAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                                                                                      >  4:198741/1‑113 (MQ=255)
                                                 acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGGCAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTa                                                                                                                                      <  3:198741/113‑1 (MQ=255)
                                                 acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAt                                                                                                            <  3:140413/139‑1 (MQ=255)
                                                 acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAAt                                                                                                            >  4:140413/1‑139 (MQ=255)
                                                 acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                                                                                  >  3:406861/1‑149 (MQ=255)
                                                 acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                                                                                  >  1:765103/1‑149 (MQ=255)
                                                   aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                                                                                >  3:751384/1‑149 (MQ=255)
                                                      aTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                                                                                             >  1:146377/1‑149 (MQ=255)
                                                           gCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttttt                                                                                        <  2:211719/149‑1 (MQ=255)
                                                            cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                                                                                       <  4:489528/149‑1 (MQ=255)
                                                                 gAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCata                                                                                  <  2:756562/149‑1 (MQ=255)
                                                                          gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatta                                                                         <  4:577266/149‑1 (MQ=255)
                                                                           cTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAg                                                                        >  2:211907/1‑149 (MQ=255)
                                                                           cTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTACTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAg                                                                        <  4:881954/149‑1 (MQ=255)
                                                                              aaaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCtt                                                                     >  4:357493/1‑149 (MQ=255)
                                                                                atatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                                                                                     >  2:377206/1‑131 (MQ=255)
                                                                                atatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                                                                                     <  1:377206/131‑1 (MQ=255)
                                                                                atatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGt                                                                   <  2:759330/149‑1 (MQ=255)
                                                                                 tatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGtt                                                                  <  1:540694/149‑1 (MQ=255)
                                                                                 tatTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGtt                                                                  <  4:406861/149‑1 (MQ=255)
                                                                                     tCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt                                                              <  1:898371/149‑1 (MQ=255)
                                                                                          tatttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGa                                                         <  4:830175/149‑1 (MQ=255)
                                                                                           atttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAg                                                        >  2:861895/1‑149 (MQ=255)
                                                                                            tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGc                                                       >  2:66397/1‑149 (MQ=255)
                                                                                            tttattCATAATAATGTTCTCACGGCAAGGGGACAGGTAACTCCGCGATGTCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGc                                                       >  4:196944/1‑149 (MQ=255)
                                                                                                ttCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTggg                                                   >  3:746064/1‑149 (MQ=255)
                                                                                                  cATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCg                                                 >  3:313061/1‑149 (MQ=255)
                                                                                                   ataataATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTTAGCTGGGCGa                                                >  2:123805/1‑149 (MQ=255)
                                                                                                   ataataATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGa                                                <  4:507891/149‑1 (MQ=255)
                                                                                                     aataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACa                                              >  1:37031/1‑149 (MQ=255)
                                                                                                     aataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACa                                              >  4:115064/1‑149 (MQ=255)
                                                                                                      ataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACag                                             >  2:262450/1‑149 (MQ=255)
                                                                                                       taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACaga                                            <  2:219356/149‑1 (MQ=255)
                                                                                                        aatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGCCagag                                           >  1:916788/1‑149 (MQ=255)
                                                                                                        aatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACagag                                           <  3:357493/149‑1 (MQ=255)
                                                                                                         atGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAgaga                                          >  2:351485/1‑149 (MQ=255)
                                                                                                               tcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTga                                                                                                                    <  4:465200/69‑1 (MQ=255)
                                                                                                               tcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTga                                                                                                                    >  3:465200/1‑69 (MQ=255)
                                                                                                                 aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaaaa                                  <  1:861895/149‑1 (MQ=255)
                                                                                                                 aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaaaa                                  <  2:37031/149‑1 (MQ=255)
                                                                                                                 aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaaaa                                  <  4:521353/149‑1 (MQ=255)
                                                                                                                  cGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAg                                 <  3:591583/149‑1 (MQ=255)
                                                                                                                   ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGc                                <  3:138244/149‑1 (MQ=255)
                                                                                                                   ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGc                                <  1:618134/149‑1 (MQ=255)
                                                                                                                       aGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTcc                            <  2:146377/149‑1 (MQ=255)
                                                                                                                        agaaaaaGGTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCa                           >  4:768558/7‑149 (MQ=255)
                                                                                                                           gACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                                                                                     >  1:245013/1‑72 (MQ=255)
                                                                                                                           gACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                                                                                     <  2:245013/72‑1 (MQ=255)
                                                                                                                             cAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAATCGTCCAGCaaa                      <  1:693886/149‑1 (MQ=255)
                                                                                                                             cAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCaaa                      <  4:184129/149‑1 (MQ=255)
                                                                                                                             cAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCaaa                      <  2:765103/149‑1 (MQ=255)
                                                                                                                              aGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAag                     >  1:102977/1‑149 (MQ=255)
                                                                                                                              aGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAag                     >  2:126668/1‑149 (MQ=255)
                                                                                                                                gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagtt                   >  2:714231/1‑149 (MQ=255)
                                                                                                                                 tAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagtta                  >  2:333499/1‑149 (MQ=255)
                                                                                                                                 tAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagtta                  >  2:121380/1‑149 (MQ=255)
                                                                                                                                  aaCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagttag                 <  4:415733/149‑1 (MQ=255)
                                                                                                                                  aaCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagttag                 >  4:575511/1‑149 (MQ=255)
                                                                                                                                  aaCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagttag                 <  3:115064/149‑1 (MQ=255)
                                                                                                                                  aaCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagttag                 <  4:751384/149‑1 (MQ=255)
                                                                                                                                    cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagttagtt               <  3:211720/149‑1 (MQ=255)
                                                                                                                                    cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagttagtt               <  1:931482/149‑1 (MQ=255)
                                                                                                                                     tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTg              >  2:234861/1‑149 (MQ=255)
                                                                                                                                        gcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTaa           <  2:102977/149‑1 (MQ=255)
                                                                                                                                        gccatcgCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTaa           <  1:333499/143‑1 (MQ=255)
                                                                                                                                              gCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCt     >  3:508503/1‑149 (MQ=255)
                                                                                                                                               cGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCtt    <  1:593342/149‑1 (MQ=255)
                                                                                                                                               cGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCtt    <  2:452995/149‑1 (MQ=255)
                                                                                                                                                 gAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGc  <  2:740271/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGC  >  NC_000913/4311942‑4312235

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