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% 527.6 / 22.1 168intergenic (‑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 (76/80);  new base G (6/6);  total (82/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 170 total aligned reads displayed.

ATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGC  >  NC_000913/4311820‑4312361
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                
                                                                 aGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGCAAAGGAATCAACAATTGAGCCGGTTTTTTGCatattaa                                                                                                                                                                                                            >  4:238620/1‑275 (MQ=255)
                                                                                                                         gagGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGaa                                                                                                                                                                                                                                                                                  <  7:35017/149‑1 (MQ=255)
                                                                                                                           gggTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAAtt                                                                                                                                                                                                                                                                                <  8:568906/149‑1 (MQ=255)
                                                                                                                            ggTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTa                                                                                                                                                                                                                                                                               <  7:1229245/149‑1 (MQ=255)
                                                                                                                            ggTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTTCGCGATGGCGGAATTACCTTTTTCTTATTCCAGCGCGCCTGGTGAGAAGAATTCAACAATTGCGCCGGTTTTTTGc                                                                                                                                                                                                                   >  2:107176/1‑209 (MQ=255)
                                                                                                                              ttttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTAcc                                                                                                                                                                                                                                                                             >  6:334397/1‑149 (MQ=255)
                                                                                                                              ttttCAATACGACAGCATATTCGGCAGCAGCAAAAGCGCTGGTCGACAACATTAAGCCCCCGACTGTGCAGCTGAAATATGTCAGATATTTATCCATAATAATGTTCTCACGGTAAGGGTACAGGTATCTCCGCGATGGCGGAATTTcc                                                                                                                                                                                                                                                                             <  8:905926/149‑1 (MQ=255)
                                                                                                                               tttCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAATTCCGCCATCGCGGAGTTACCt                                                                                                                                                                                                                                                                            <  8:160869/149‑1 (MQ=255)
                                                                                                                                 tCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttt                                                                                                                                                                                                                                                                          >  8:597709/1‑149 (MQ=255)
                                                                                                                                  cAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCtttt                                                                                                                                                                                                                                                                         <  7:1126867/149‑1 (MQ=255)
                                                                                                                                   aaTACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCttttt                                                                                                                                                                                                                                                                        <  7:736327/149‑1 (MQ=255)
                                                                                                                                    aTACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTc                                                                                                                                                                                                                                                                       <  5:1443588/149‑1 (MQ=255)
                                                                                                                                      aCGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCtt                                                                                                                                                                                                                                                                     >  5:1382349/1‑149 (MQ=255)
                                                                                                                                        gACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAt                                                                                                                                                                                                                                                                   >  7:552787/1‑149 (MQ=255)
                                                                                                                                        gACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAt                                                                                                                                                                                                                                                                   >  7:297244/1‑149 (MQ=255)
                                                                                                                                        gACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAt                                                                                                                                                                                                                                                                   >  5:179234/1‑149 (MQ=255)
                                                                                                                                         aCAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAtt                                                                                                                                                                                                                                                                  >  5:843835/1‑149 (MQ=255)
                                                                                                                                          cAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTAttt                                                                                                                                                                                                                                                                 >  8:765080/1‑149 (MQ=255)
                                                                                                                                              atcttCGGCGGCAGTAAAAGCGCTGGTTGACAACCTTAAGCCCACGAGTGTGCCGCAGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgc                                                                                                                                                                                                                                                             <  8:1217357/146‑1 (MQ=255)
                                                                                                                                              atatTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAgc                                                                                                                                                                                                                                                             <  5:430142/149‑1 (MQ=255)
                                                                                                                                                    ggcggcAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCt                                                                                                                                                                                                                                                       <  5:334397/149‑1 (MQ=255)
                                                                                                                                                    ggcggcAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCt                                                                                                                                                                                                                                                       <  6:182604/149‑1 (MQ=255)
                                                                                                                                                     gcggcAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCt                                                                                                                                                                                                                                                       <  5:901913/148‑1 (MQ=255)
                                                                                                                                                       ggcAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGt                                                                                                                                                                                                                                                    <  5:1340434/149‑1 (MQ=255)
                                                                                                                                                         cagcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCCTTTTCTTATTTCAGCGCGCCTGGTga                                                                                                                                                                                                                                                  <  7:765080/149‑1 (MQ=255)
                                                                                                                                                          agcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgag                                                                                                                                                                                                                                                 >  6:516896/1‑149 (MQ=255)
                                                                                                                                                           gcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                                                                                                                                                                                                                >  6:1416736/1‑149 (MQ=255)
                                                                                                                                                           gcaAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTgaga                                                                                                                                                                                                                                                <  5:589872/149‑1 (MQ=255)
                                                                                                                                                            caAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTGAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaa                                                                                                                                                                                                                                               <  8:521638/149‑1 (MQ=255)
                                                                                                                                                            caAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaa                                                                                                                                                                                                                                               <  5:75782/149‑1 (MQ=255)
                                                                                                                                                             aaaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAACTATTTCAGATATTTATTCATAATAATGTTCTCACGGAAAGGGGACAGCTAACTCGGCGATGGCGGAATTACCTTTTTTTTATTTCAGCGCGCCTGGTGagaag                                                                                                                                                                                                                                              >  6:92283/1‑149 (MQ=255)
                                                                                                                                                             aaaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                                                                                                                                                                                                              >  7:1052338/1‑149 (MQ=255)
                                                                                                                                                             aaaaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGagaag                                                                                                                                                                                                                                              >  8:1408155/1‑149 (MQ=255)
                                                                                                                                                               aaGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAgaagaa                                                                                                                                                                                                                                            >  5:658889/1‑149 (MQ=255)
                                                                                                                                                                  cgcTGGTTGACAACATTAAGCCCACGAGCGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGCGAAGAAATc                                                                                                                                                                                                                                         <  8:1256634/149‑1 (MQ=255)
                                                                                                                                                                    cTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaa                                                                                                                                                                                                                                       <  8:935365/149‑1 (MQ=255)
                                                                                                                                                                     tGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATCTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaac                                                                                                                                                                                                                                      >  6:1322080/1‑149 (MQ=255)
                                                                                                                                                                       gTTGACAACATTAAGCCCACGAGTGAGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaacaa                                                                                                                                                                                                                                    >  7:204349/1‑149 (MQ=255)
                                                                                                                                                                         tGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAAtt                                                                                                                                                                                                                                  <  8:297244/149‑1 (MQ=255)
                                                                                                                                                                           acaacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgc                                                                                                                                                                                                                                <  6:73597/149‑1 (MQ=255)
                                                                                                                                                                             aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTTATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCGTTTTCTTATTTCAGCGCGCCTGGTGAGAACAAATCAACAATTgcgc                                                                                                                                                                                                                              >  8:1186885/1‑149 (MQ=255)
                                                                                                                                                                             aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaacaa                                                                                                                                                                                                                                    <  5:544350/143‑1 (MQ=255)
                                                                                                                                                                             aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATcaacaa                                                                                                                                                                                                                                    >  6:544350/1‑143 (MQ=255)
                                                                                                                                                                             aacaTTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                                                                                                                                                                                                              >  8:465784/1‑149 (MQ=255)
                                                                                                                                                                                 ttAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGt                                                                                                                                                                                                                          >  6:823421/1‑149 (MQ=255)
                                                                                                                                                                                     gCCCACGAGTGTGCCACTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACTGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGAttttt                                                                                                                                                                                                                      <  8:401336/149‑1 (MQ=255)
                                                                                                                                                                                         aCGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGACGGTTTTTTGCATATTAATCTTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTCCGACCTACTAACGGCGTTTCCGGGGCTGGTGTGc                                                                                    >  2:348524/1‑275 (MQ=255)
                                                                                                                                                                                          cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatat                                                                                                                                                                                                         >  3:37844/1‑157 (MQ=255)
                                                                                                                                                                                          cGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCat                                                                                                                                                                                                                 <  8:1090254/149‑1 (MQ=255)
                                                                                                                                                                                              tgtgCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc                                                                                                                                                                                                     <  1:37844/157‑1 (MQ=255)
                                                                                                                                                                                                      tGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGc                                                                                                                                                                                                     >  7:1161216/1‑149 (MQ=255)
                                                                                                                                                                                                       gAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCt                                                                                                                                                                                                    <  6:802737/149‑1 (MQ=255)
                                                                                                                                                                                                         aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTgcgc                                                                                                                                                                                                                              <  6:793416/121‑1 (MQ=255)
                                                                                                                                                                                                         aaTATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTg                                                                                                                                                                                                  >  8:636722/1‑149 (MQ=255)
                                                                                                                                                                                                            atTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTCGTTc                                                                                                                                                                                               <  6:438670/149‑1 (MQ=255)
                                                                                                                                                                                                              ttCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCaa                                                                                                                                                                                             <  7:1105532/149‑1 (MQ=255)
                                                                                                                                                                                                               tCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatat                                                                                                                                                                                                         >  5:829868/1‑136 (MQ=255)
                                                                                                                                                                                                               tCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaatat                                                                                                                                                                                                         <  6:829868/136‑1 (MQ=255)
                                                                                                                                                                                                               tCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAAt                                                                                                                                                                                            <  6:845338/149‑1 (MQ=255)
                                                                                                                                                                                                                 cgATTTTTATTCATAATAAGTTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGGATTACCTTTTTCTTGTTTCAGCGCGACTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCg                                                            <  4:15883/274‑1 (MQ=255)
                                                                                                                                                                                                                 aGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGt                                                                                                                                                                                          >  5:315774/1‑149 (MQ=255)
                                                                                                                                                                                                                 aGATATTTATTAATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGt                                                                                                                                                                                          >  5:1097852/1‑149 (MQ=255)
                                                                                                                                                                                                                   atatTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTg                                                                                                                                                                                        <  5:516896/149‑1 (MQ=255)
                                                                                                                                                                                                                      tttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGc                                                                                                                                                                                     >  6:1303904/1‑149 (MQ=255)
                                                                                                                                                                                                                       ttattCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCt                                                                                                                                                                                    >  5:849906/1‑149 (MQ=255)
                                                                                                                                                                                                                         attCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTgg                                                                                                                                                                                  <  7:1408155/149‑1 (MQ=255)
                                                                                                                                                                                                                          ttCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTggg                                                                                                                                                                                 <  7:50145/149‑1 (MQ=255)
                                                                                                                                                                                                                          ttCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTggg                                                                                                                                                                                 <  7:468265/149‑1 (MQ=255)
                                                                                                                                                                                                                                 taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatat                                                                                                                                                                                                               <  8:696367/112‑1 (MQ=255)
                                                                                                                                                                                                                                 taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatat                                                                                                                                                                                                               >  7:696367/1‑112 (MQ=255)
                                                                                                                                                                                                                                 taatGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACaga                                                                                                                                                                          <  7:1122540/149‑1 (MQ=255)
                                                                                                                                                                                                                                     gTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACa                                                                                                                                                                      >  8:817221/1‑149 (MQ=255)
                                                                                                                                                                                                                                     gTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACa                                                                                                                                                                      >  5:233214/1‑149 (MQ=255)
                                                                                                                                                                                                                                      ttCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCatattaat                                                                                                                                                                                                           >  6:1345017/1‑111 (MQ=255)
                                                                                                                                                                                                                                      ttCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaa                                                                                                                                                                     >  7:1379156/1‑149 (MQ=255)
                                                                                                                                                                                                                                      ttCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaa                                                                                                                                                                     <  8:644502/149‑1 (MQ=255)
                                                                                                                                                                                                                                       tctcACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaa                                                                                                                                                                    >  7:1036426/1‑149 (MQ=255)
                                                                                                                                                                                                                                           aCGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACaaaaaaa                                                                                                                                                                <  6:116915/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggggACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTCTTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGTCGTTTCAGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTTAATGGTAGCTTTGGTTTATACaaa                           >  4:32490/1‑275 (MQ=255)
                                                                                                                                                                                                                                                  agaaaaaGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCa                                                                                                                                                         <  5:1019974/143‑1 (MQ=255)
                                                                                                                                                                                                                                                        aGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAgaga                                                                                                                                                                        <  5:40569/128‑1 (MQ=255)
                                                                                                                                                                                                                                                        aGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAgaga                                                                                                                                                                        >  6:40569/1‑128 (MQ=255)
                                                                                                                                                                                                                                                          gTAATTCCGCCATCGCGGAGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagtt                                                                                                                                                 >  6:596361/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagtt                                                                                                                                                 <  8:1161216/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          gTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagtt                                                                                                                                                 >  5:880310/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           tAATTCCGCCATCGCGGCGTTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagtta                                                                                                                                                >  6:1383082/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           tAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagtta                                                                                                                                                <  6:1179079/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           tAACTCCGCGATAGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGc                                                                                                                                                                                                                   <  7:957582/82‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaCTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagttag                                                                                                                                               >  7:94103/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagttagtt                                                                                                                                             <  6:233214/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              cTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAagttagtt                                                                                                                                             <  8:657963/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGtt                                                                                                                                                                                         <  8:596562/104‑1 (MQ=255)
                                                                                                                                                                                                                                                               tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGtt                                                                                                                                                                                         >  7:596562/1‑104 (MQ=255)
                                                                                                                                                                                                                                                               tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACaga                                                                                                                                                                          <  5:551811/119‑1 (MQ=255)
                                                                                                                                                                                                                                                               tCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTTACTGGTAGCTATGGGTTATCACAATAATCTTTGccc               >  2:212867/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                  gcgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTaa                                                                                                                                         <  7:636722/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cgATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACATAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATCTAATCGCCAAGTTTCGCTATCGTTGACTGGTAGCTATGGTTTATACAAATAATCTTTGCCCCGt            >  4:209266/1‑273 (MQ=255)
                                                                                                                                                                                                                                                                         cGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCtt                                                                                                                                  <  5:15826/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                          ggAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCACAGTTAGTTGTAATATTCTTg                                                                                                                                 >  6:397211/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                           gAATTACCTTGTTCGTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATAGTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATAGTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGc  <  2:32490/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                
ATAACTGCAATTGAGATTGAAAATCGCCTTCTGAAGGAGAGGCAAAAATATCAACGCTGACGCCCAGTGTTTTTGCTTCATCTTCAATGCCTTTTTTCATATCTACCCAAAATGGGTTGGAGAGGGTTTTCAATACGACAGCATATTCGGCGGCAGCAAAAGCGCTGGTTGACAACATTAAGCCCACGAGTGTGCCGCTGAAATATTTCAGATATTTATTCATAATAATGTTCTCACGGTAAGGGGACAGGTAACTCCGCGATGGCGGAATTACCTTTTTCTTATTTCAGCGCGCCTGGTGAGAAGAAATCAACAATTGCGCCGGTTTTTTGCATATTAATATTAGCTTGTTCAATGTTGAGCTGGGCGACAGAGACAAAAAAAGCGTCCAGCAAAGTTAGTTGTAATATTCTTGCCGAGGCATTACGACCTAATAACGGCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGC  >  NC_000913/4311820‑4312361

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