Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,667,460 G→A Q369* (CAG→TAG)  mlc ← DNA‑binding transcriptional repressor Mlc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,667,4600GA98.0% 166.0 / ‑3.8 51coding (1105/1221 nt)mlcDNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/19);  minor base . (1/0);  total (32/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTGAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAAC  >  NC_000913/1667216‑1667698
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
gCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCACAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑tata                                                                                                                                                                                                                                           <  1:568234/251‑1 (MQ=255)
       aaaCTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGa                                                                                                                                                                                                                                    <  2:317231/250‑1 (MQ=255)
         aCTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGAGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGa                                                                                                                                                                                                                                    >  2:147209/1‑249 (MQ=255)
               cctttATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACTCTGATGTGCTAAC‑TATACGCAGGAAGGGCctg                                                                                                                                                                                                                            <  1:275941/248‑1 (MQ=255)
                    aTGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGAc                                                                                                                                                                                                                       <  1:594499/251‑1 (MQ=255)
                        ttCACAAATATAAACCAGGAAAATAATTAACCTTGCAAGTCTAAGTTATGATTTCCTGGCCCAAATTAAGATAGCGCAAATTTTGGTAGAACAGTTCAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGACAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACggct                                                                                                                                                                                                                   >  2:109266/1‑249 (MQ=255)
                          cacaAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTCGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGc                                                                                                                                                                                                                 >  2:474023/1‑251 (MQ=255)
                                   aaaCCAGGAAAATAATTAACCTTGAAAGTCTCAGTTATGCTTTCCTGGCCCAAAATGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTACCCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCCGTGCAGCGCCTGCCATCGTGCACTGGTTAGAAAAC‑TGAGTACTCTCAACGATGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGCCTg                                                                                                                                                                                                           >  2:637122/1‑248 (MQ=255)
                                    aaCCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACATTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTGTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACGCTCAACGCTGAGTTGCTAAC‑TATACGCAGGAAGGGCCTTCTGACGGATGCTGTCTGAGAt                                                                                                                                                                                                       <  2:38605/251‑1 (MQ=255)
                                                aatTAACCTTGAAAGTCTAAGTTATGCTTTCATGGCCCAAATTGAGATATCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAaga                                                                                                                                                                                            <  1:186429/250‑1 (MQ=255)
                                                 atTAACCTTGAAAGACTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAACC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGCATGCTGTCTGAGATGACCGGGAAGAgg                                                                                                                                                                                          >  2:562504/1‑251 (MQ=255)
                                                     acattGAAAGTCTAAGTTATGATTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAATATGTTAACCATGCTACAGACGAATCAACGAAGAACCGTTATACATCGCGTCTTTTACCAGTGAAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCCACGCTGATGTGCTAAA‑TATACGCAGGAAGGGCCTGATGACGGATGCTGTCTGAGAAGCCCGGGaaagggaaat                                                                                                                                                                                      >  2:344853/4‑243 (MQ=255)
                                                        ttGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAAAAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTAATCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTg                                                                                                                                                                                  >  1:265540/1‑251 (MQ=255)
                                                           aaaGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCg                                                                                                                                                                               <  1:95697/251‑1 (MQ=255)
                                                             aGTCTAAATTAGGCTTTCCTGGCCCAAAGTGATATAGCGCAAATTTTGGTAGAACATTAAAAAAATTTTATCCCTTCATCAGACGAATCACCAAAGAACCGTTAAACTTCGCGTCTTTTACCAGTGCAGCGCCTTCCAGCGTGCCCTGTTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGTGCCTGCTGTCGGATGCTGTCTGAGATGACCGGGAAGAGTA‑TTTCTGCCGCt                                                                                                                                                                             <  2:265540/251‑1 (MQ=255)
                                                                  aaGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCttt                                                                                                                                                                           <  1:405903/248‑1 (MQ=255)
                                                                     ttATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGCTGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTa                                                                                                                                                                      >  1:415134/1‑250 (MQ=255)
                                                                      tATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAAc                                                                                                                                                                    >  1:795368/1‑251 (MQ=255)
                                                                       aTGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAAc                                                                                                                                                                    <  1:562504/250‑1 (MQ=255)
                                                                            ttCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTg                                                                                                                                                                >  1:166370/1‑249 (MQ=255)
                                                                                  gCCCAAATTGAGATAGCGCAACTTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGAAAGGGCATGCTGACGGATGCTGTCTGAGATGACAGGGAAGAGGA‑TATCTGCAGCTTTACTTAACGGGGAGCCAATc                                                                                                                                                        >  2:810633/1‑251 (MQ=255)
                                                                                  gCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATc                                                                                                                                                        <  1:774639/251‑1 (MQ=255)
                                                                                       aaTTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTAt                                                                                                                                                   >  2:210695/1‑251 (MQ=255)
                                                                                          tGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTAtttt                                                                                                                                                >  1:598063/1‑251 (MQ=255)
                                                                                          tGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCGCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTAttt                                                                                                                                                 <  2:264554/250‑1 (MQ=255)
                                                                                                            ggTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTaaataaat                                                                                                                               <  1:63260/250‑1 (MQ=255)
                                                                                                                 aaCAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCAcc                                                                                                                          >  1:672651/1‑250 (MQ=255)
                                                                                                                                  aacacTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATgacc                                                                                                        >  2:9824/5‑248 (MQ=255)
                                                                                                                                      cagCAACAGACGAATCAACCAGGAACCGTTATACATCGCGTCTTTTACAAGTGCAGCGCCTGCCATCGGGCCCTGGTTAGAAAAC‑TGAGTATTCTCAACGCTGATGTGCTACC‑TATACGCCGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAGGAGGAAGAT‑TGCCGCTTTACTTAACGGTGAGCCAATCAGGATTTTTTGTGGGTTAAATACATTCACCatatagtaaaaaaagcgcacg                                                                                                     >  2:51989/3‑231 (MQ=255)
                                                                                                                                              gACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATgcgcg                                                                                            >  2:766909/1‑251 (MQ=255)
                                                                                                                                               aCGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGCGCCAATCAGTATTTTTTGTGGGTTAAATAAATTAACCATGATGGCAAGAATGCGCCCGACATgcgcgc                                                                                           >  1:336989/1‑251 (MQ=255)
                                                                                                                                                   aTCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACCTGCGCGCCCAc                                                                                       >  2:596352/1‑251 (MQ=255)
                                                                                                                                                                cGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCACGACATGCGCGCCCACACCGGTAatgatg                                                                          >  2:133974/1‑251 (MQ=255)
                                                                                                                                                                 gTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGt                                                                         >  1:600101/1‑251 (MQ=255)
                                                                                                                                                                     taCATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCGGGAAGGGCCTGCTGCCGGATGCTGTCTGAGCTGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCCATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGACAGAATGCGCCCGACATGCGCGCCCACCCCGGTCATGATGTCttt                                                                     >  2:452027/1‑251 (MQ=255)
                                                                                                                                                                          cgagTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGcc                                                                 <  1:147209/247‑1 (MQ=255)
                                                                                                                                                                            cgTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGt                                                              >  2:196561/1‑251 (MQ=255)
                                                                                                                                                                                      cAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATcgccgc                                                    >  2:255202/1‑251 (MQ=255)
                                                                                                                                                                                           cAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAAt                                               <  1:224597/251‑1 (MQ=255)
                                                                                                                                                                                              cgcCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATgcc                                            >  1:161775/1‑251 (MQ=255)
                                                                                                                                                                                               gaatggcATCGTGGCCTGGTTA‑ACAACATGGGTACTCTCCACGTTGCTGCGCTACC‑TATCCGCAGGAAGGGCCTGCTGACGGATGATGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGCGCCAATCAGTATTTTTTGTGGGCTAAATCACTTCACCATGCTGGCAAGAATGCGACCGCCATGCAAGCCACCCCCGGTAATGATGTCTTTTGCCCGTAGATCGCCgcgccctgccc                                           >  2:518135/7‑244 (MQ=255)
                                                                                                                                                                                                ccTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATgccgc                                          >  2:514484/1‑251 (MQ=255)
                                                                                                                                                                                                  tGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCt                                        >  1:170865/1‑251 (MQ=255)
                                                                                                                                                                                                   gCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCt                                        <  1:596352/250‑1 (MQ=255)
                                                                                                                                                                                                   gCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTg                                       <  1:717021/251‑1 (MQ=255)
                                                                                                                                                                                                                ggTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTc                           <  1:804155/250‑1 (MQ=255)
                                                                                                                                                                                                                  ttAGAAAAC‑TGAGTACTCTAAACGCTGATGTGCT‑ACATATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATAAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGAAAGAATGCGCCCGACATGCGCACCCACCCCGGAAAAGATGTCTTTTGCCAGTCGCTCGCCGCTCAATGCCGCCTGACCCAAAGAATCCAc                        >  2:550844/1‑251 (MQ=255)
                                                                                                                                                                                                                   tAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCAc                        <  1:210695/250‑1 (MQ=255)
                                                                                                                                                                                                                                   tctcAACGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGCTGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGAAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTcc         >  1:758313/1‑250 (MQ=255)
                                                                                                                                                                                                                                         cGCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAAc  >  1:739122/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gCTGATGTGCTAAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAAc  >  2:132082/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                
GCATCACTAACTCCACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAAC‑TGAGTACTCTCAACGCTGATGTGCTGAC‑TATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGA‑TATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTTGTCCATGTAAC  >  NC_000913/1667216‑1667698

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