Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,312,089 A→G 7.1% 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 7.1% 222.0 / 6.4 70intergenic (‑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 (28/37);  new base G (3/2);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

TTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTA  >  NC_000913/4311946‑4312225
                                                                                                                                               |                                                                                                                                        
ttttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTa                                                                                                                                       >  4:100545/1‑147 (MQ=255)
ttttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGACCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTa                                                                                                                                       <  3:100545/147‑1 (MQ=255)
  ttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCt                                                                                                                                    <  2:894018/148‑1 (MQ=255)
    cAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCtttt                                                                                                                                 <  4:109396/149‑1 (MQ=255)
    cAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttt                                                                                                                                  >  1:295128/1‑148 (MQ=255)
      aTACGACAGCATATTCGGCTGCAGCAAAAGCGCTGGTTGACAGCATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                                                               >  3:532210/1‑149 (MQ=255)
      aTACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                                                                                >  2:141198/1‑148 (MQ=255)
           aCAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAt                                                                                                                           >  2:497529/1‑148 (MQ=255)
              gCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCa                                                                                                                       <  4:321967/149‑1 (MQ=255)
                 tatTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgc                                                                                                                     >  1:857904/1‑148 (MQ=255)
                      ggcggcAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCACTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGcc                                                                                                                <  1:753324/148‑1 (MQ=255)
                              caAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                                                                        >  1:859483/1‑148 (MQ=255)
                              caAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                                                                        >  2:112553/1‑148 (MQ=255)
                              caAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                                                                        <  2:785373/148‑1 (MQ=255)
                                 aaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAgaagaa                                                                                                    <  3:202446/149‑1 (MQ=255)
                                   gcgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGaaa                                                                                                   <  1:478879/148‑1 (MQ=255)
                                   gcgcTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGaaa                                                                                                   <  2:857904/148‑1 (MQ=255)
                                    cgcTGGTTGACAACATTAAGCCCACGAGAGTGCCGCTGAAATATTTCAGCTATTTATTCATAATAATGTTCTCACCGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATc                                                                                                 <  3:765692/149‑1 (MQ=255)
                                      cTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATca                                                                                                <  1:399560/148‑1 (MQ=255)
                                          ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaacaa                                                                                            <  2:333497/148‑1 (MQ=255)
                                          ttGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAt                                                                                           <  4:752636/149‑1 (MQ=255)
                                             acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                                                                        >  4:639330/1‑149 (MQ=255)
                                              caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                                                                               >  3:41750/1‑141 (MQ=255)
                                              caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                                                                               <  4:41750/141‑1 (MQ=255)
                                              caacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                                                                        >  2:801556/1‑148 (MQ=255)
                                               aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcg                                                                                       <  1:141198/148‑1 (MQ=255)
                                                acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCt                                                                                                               <  1:907146/123‑1 (MQ=255)
                                                acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCt                                                                                                               >  2:907146/1‑123 (MQ=255)
                                                acaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGcc                                                                                     <  3:691472/149‑1 (MQ=255)
                                                   ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCgg                                                                                   >  1:168705/1‑148 (MQ=255)
                                                    tAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGt                                                                                  >  2:772464/1‑148 (MQ=255)
                                                        cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                                                                             <  4:247302/149‑1 (MQ=255)
                                                        cccACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttttt                                                                              <  1:801556/148‑1 (MQ=255)
                                                          cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                                                                             >  4:329197/1‑147 (MQ=255)
                                                          cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGtttttt                                                                             <  3:329197/147‑1 (MQ=255)
                                                          cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                                                                           >  4:468109/1‑149 (MQ=255)
                                                          cACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                                                                           >  3:105309/1‑149 (MQ=255)
                                                            cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCa                                                                          >  2:645589/1‑148 (MQ=255)
                                                             gAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCat                                                                         <  1:772464/148‑1 (MQ=255)
                                                                  tgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaa                                                                    >  1:421511/1‑148 (MQ=255)
                                                                      gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatta                                                               <  4:105309/149‑1 (MQ=255)
                                                                      gCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatatt                                                                <  1:112553/148‑1 (MQ=255)
                                                                        tGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc                                                             <  3:639330/149‑1 (MQ=255)
                                                                          aaaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt                                                            <  1:645589/148‑1 (MQ=255)
                                                                                       atttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAg                                              <  3:1250/149‑1 (MQ=255)
                                                                                           attCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTg                                           <  1:497529/148‑1 (MQ=255)
                                                                                            ttCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTggg                                         >  4:480149/1‑149 (MQ=255)
                                                                                                 aataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttttt                                                                              >  3:162539/1‑107 (MQ=255)
                                                                                                 aataatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGttttt                                                                              <  4:162539/107‑1 (MQ=255)
                                                                                                   taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACaga                                  >  3:827560/1‑149 (MQ=255)
                                                                                                     atGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACagag                                 <  2:168705/148‑1 (MQ=255)
                                                                                                            cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaa                          >  2:385182/1‑148 (MQ=255)
                                                                                                            cACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaa                          >  1:377042/1‑148 (MQ=255)
                                                                                                             aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaaaa                        <  4:532210/149‑1 (MQ=255)
                                                                                                              cGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAACGTTGAGCTGGGCGACAGAGACaaaaaaa                        >  2:416697/1‑148 (MQ=255)
                                                                                                               ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTTAATGTTGAGCTGGGCGACAGAGACAAAAAAAGc                      <  3:572319/149‑1 (MQ=255)
                                                                                                               ggTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAg                       <  2:859483/148‑1 (MQ=255)
                                                                                                                         cAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCaa             <  1:385182/148‑1 (MQ=255)
                                                                                                                           ggTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAGGCGTCCAGCAAagt          >  3:404121/1‑149 (MQ=255)
                                                                                                                           ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt                                                    >  4:466932/1‑107 (MQ=255)
                                                                                                                           ggTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt                                                    <  3:466932/107‑1 (MQ=255)
                                                                                                                            gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaaaa                        >  3:322098/1‑134 (MQ=255)
                                                                                                                            gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaaaa                        <  4:322098/134‑1 (MQ=255)
                                                                                                                             tAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAg                       <  1:669407/134‑1 (MQ=255)
                                                                                                                             tAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAg                       >  2:669407/1‑134 (MQ=255)
                                                                                                                              aattCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagtta        >  1:209672/4‑148 (MQ=255)
                                                                                                                               attCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAGGCGTCCAGCAAagttagt      <  4:404121/147‑1 (MQ=255)
                                                                                                                               aCTCCGCGATGGCGGAATTAACTTTTTTTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGGCGGTTTTTTGCATATTTATATTTGCCTGTTCCATGGTGAAGCGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagttagt      <  3:480149/149‑1 (MQ=255)
                                                                                                                                 tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTg    >  3:68998/1‑149 (MQ=255)
                                                                                                                                    gcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTa  <  2:295128/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                        
TTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTA  >  NC_000913/4311946‑4312225

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