Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 70,371 T→C intergenic (‑323/‑16) araB ← / → araC ribulokinase/DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091370,3710TC97.2% 235.9 / NA 73intergenic (‑323/‑16)araB/araCribulokinase/DNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base C (37/33);  minor base A (2/0);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGT  >  NC_000913/70128‑70606
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
tCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGt                                                                                                                                                                                                                                            >  1:945614/1‑250 (MQ=255)
   tATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTa                                                                                                                                                                                                                                           <  1:48662/248‑1 (MQ=255)
     tGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATATACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCCGTAAAAGACGCAGTGACGGCATTGTCTGATGCAATATGGCCCATTGGTTTCTTCTCTGACTGGCGGGAGTATGac                                                                                                                                                                                                                                       >  2:16050/1‑249 (MQ=255)
       gATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCTGGAGTATGAAAAg                                                                                                                                                                                                                                    <  1:780488/251‑1 (MQ=255)
           aaaaGGCTATGGACTAGCAAAGTGTGACGACGTGCAATTAATCAATGTGGACTTTTCTGCCGTGATTATAGACACCTTTGTTTCGCGTTTATGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCAGTGAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGACAGTAAAAGACGCAGTGACGGAAATGTCTGATGCAAGATGGGCAATTGGTGTCTTCTTGGAATGGCGGGAGAATGAAAAGTATg                                                                                                                                                                                                                                <  1:795587/251‑1 (MQ=255)
                     ggCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGTCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAgcg                                                                                                                                                                                                                       <  1:16050/250‑1 (MQ=255)
                     ggCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAgcgc                                                                                                                                                                                                                      >  1:600877/1‑251 (MQ=255)
                        aTAGCAAAGTGGGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGTGTTTTTGGCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGTGGTTACCGGTTTGGTTAGCGAGAAGAGTCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCaaa                                                                                                                                                                                                                   <  1:920834/251‑1 (MQ=255)
                             aaaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGAt                                                                                                                                                                                                              <  1:984034/251‑1 (MQ=255)
                              aaGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATc                                                                                                                                                                                                             <  1:53138/251‑1 (MQ=255)
                                    gtagactagcaatgacacAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGTTTTTTTGTCATGGCTTGGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTTACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGTTAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCctg                                                                                                                                                                                                        <  2:386886/233‑1 (MQ=255)
                                    gACGCCGTGCAAATAATCTATGTGGACTTTTCTGCCGTGATTATCGTCACCTTTGTTACGCTTTTGTGTCATGTCTTTGGTCCCGCTTTGTTACAGAATTATTTTAATAAGCGGTGTTACCGGTTTGGTTAGCGAGAATAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATTGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCctgc                                                                                                                                                                                                       <  1:566832/251‑1 (MQ=255)
                                    gACGCCGTGCAAATAATCAATGTGGACTTTTCTGACGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGATTTGGTTCCACTTTGTTACAGAATGCTTTTAATAAGTGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGCCGGCAGTGTTTGAGGCAAAATGGACAATTGGTTTATGTTCTGGTTGGCGGGTGTATGAAAAGGATGGCTgtgtctatgactgaatacctgc                                                                                                                                                                                                       >  2:198652/1‑230 (MQ=255)
                                            gCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGa                                                                                                                                                                                               <  1:667881/251‑1 (MQ=255)
                                                         gtgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGAGACTCGTTTaatg                                                                                                                                                                                  <  1:732763/251‑3 (MQ=255)
                                                              cTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCAt                                                                                                                                                                             >  1:176331/1‑251 (MQ=255)
                                                                  tCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCtg                                                                                                                                                                          <  1:380897/250‑1 (MQ=255)
                                                                        gTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGAGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCCTCTGGAGGCgg                                                                                                                                                                   >  1:920822/1‑251 (MQ=255)
                                                                          gATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCACGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGAGGGAGGATGAAAAGTATGGCTGTAGCGCACAATGATCCCCCGCTGCCGGGGTACTCGGTTattacccattcgtgggcgggg                                                                                                                                                                 >  2:43314/1‑231 (MQ=255)
                                                                          gATTATAGACACTTTTGTTACGCGTATTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGTATGAAAAGTATGGCTGAAGCGCACAATGATCCCCTGCTGCCGGGATACTCGATTAATGCCCATCTGGTGGCgggg                                                                                                                                                                 >  1:707550/1‑250 (MQ=255)
                                                                              ataGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGGAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTa                                                                                                                                                             >  1:441589/1‑250 (MQ=255)
                                                                                  acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGcc                                                                                                                                                       <  1:504924/251‑1 (MQ=255)
                                                                                   cacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCg                                                                                                                                                      <  1:125225/251‑1 (MQ=255)
                                                                                                ggTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACg                                                                                                                                        <  2:258970/250‑1 (MQ=255)
                                                                                                 gTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGGATGAAAAGTATTGGTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTTATGCCCATCTGGTGGCGGGttttaacccgatttggggccactg                                                                                                                                         >  2:163689/1‑230 (MQ=255)
                                                                                                 gTTTTTGTCATGGCTTTGGGCCCGTTTTGTTAAAGAATGCGTTTAATAAGCGGGGTTTCCGGTTTGGGTAGCGAGAATAGCCAGGAAAAGACGCAGTTACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGACAACgg                                                                                                                                       <  2:461978/251‑1 (MQ=255)
                                                                                                   ttttGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGGATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGACGGTTGAGGCCAACGGtt                                                                                                                                     >  1:235800/1‑251 (MQ=255)
                                                                                                     ttGTCATGGCTTTGGTCCCGCTTTGTTACAGCATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACCATTGGTTTCTTCTCTGAATGGCGGGAGTATGGAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTTCCGGGGAACTCGCTTAATGCCCATCTGGGGGGGGGTTTACACCCGATTTTGGCCAACGGTTAt                                                                                                                                   >  2:193338/1‑251 (MQ=255)
                                                                                                      tGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGGTACTCGTTTAATGCCCATCTGGTGGCGGGGTTAAAGCCGATTGAGGcccaagggtatc                                                                                                                                  >  1:990605/1‑241 (MQ=255)
                                                                                                      tGTCATGGCTTTGGGCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTAt                                                                                                                                   <  1:896872/250‑1 (MQ=255)
                                                                                                          aTGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAACGCCGGTTGAGGCCAAAGGGTATCTCgg                                                                                                                              >  1:879126/1‑250 (MQ=255)
                                                                                                                 tGGTCCCGCTTTGTGACAGAATGCTTTAAATAAGCGGGGTTACCGAGTTGGTTAGCGAGAAGAGCCAGTAAAAGGCGCAGTGACGGCAATGTCTGATGCAATATGGACAAGTGGTTTCTTCTCTGAATGGCGGAAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCCTCTGGTGGCGGATTTAACGCCGATTGAGGCCAACGGTTATCTCGAttgttt                                                                                                                        <  2:1203/250‑1 (MQ=255)
                                                                                                                      ccGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACTGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGGGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAATCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATGTTTTTATcga                                                                                                                  <  2:945456/251‑1 (MQ=255)
                                                                                                                      ccGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGGTTTAGGCCAACGGTTTTCTCGGTTTTTTTTTcgg                                                                                                                  >  1:749589/1‑250 (MQ=255)
                                                                                                                      ccGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGGATGGCTGTAGCGCAAAATGATCCCCCGCTGCCGGGGTACTCGGTTAATGCCCATCTGGTGGCGGGGTTTACCCCCGTTTAAGGCCAACGGTTTCCCGGttttttttatcg                                                                                                                  >  2:332508/1‑246 (MQ=255)
                                                                                                                        gCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGCGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTCTGAAAAGAAAGGCTGAAGCGCAAAATGTTCCCCTGCTGCCGGGGTAACAGTTTAATCCCCATTTGGTGGCGGGTTTAACGACGATGGAGGCCAAAGGTTATCTCGGtttttttttcctcc                                                                                                                >  2:307309/1‑244 (MQ=255)
                                                                                                                         cTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGAGCAAAATGATCCCCTGCTGCAGGGATTCTCGTTTAATGCCCATCTGGTGGCGGGGTTAACCCCGATTGAGGCCAACGGGTTTTTCGAttttttttttgcccg                                                                                                               >  1:783817/1‑243 (MQ=255)
                                                                                                                            tGTTACAGAATGCTTTTAATAAGCGGGGATACAGGTTTGATTAGCGAGAAGAGCCAGTAAAAGAAGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGAGGGAGGATGAAAAGTATGACTGAAGAGCAAAATGATCCCCTGCTGCCGGGATAATCGGTTAATGCCCATATGGTGGCGGGCTTAACGCAGGTTGAGGCCAAAGGTTATATCGGGTTTTTTATCgacagacc                                                                                                           >  1:850027/1‑251 (MQ=255)
                                                                                                                                 cAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGGGAATATCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTGCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGGGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTgg                                                                                                      <  2:920822/251‑1 (MQ=255)
                                                                                                                                 cAGAATGCTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGGGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAACGCCGATTGGGGCCCACGGTTATCTCGGTTTTTTTTTCGACCCACCCCTgg                                                                                                      >  1:122616/1‑251 (MQ=255)
                                                                                                                                        cTTTTAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCCAAATGATCCCCTGCTGCCGGGATACTCGTTTAATTCCCATCTGGGGGCGGGTTTAACGCCGATTGAGGCCAACGGGTATCTCGAtttttttttccgccggacgatggggaata                                                                                                 >  1:594829/1‑228 (MQ=255)
                                                                                                                                          gttAATAAGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGaaa                                                                                              <  2:176331/249‑1 (MQ=255)
                                                                                                                                            taataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGGCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTTAATGGCGGGAGTATTAAAAGTATGGCTGTAGCGCAAAATGATCCCCTGCTGCCGGGATACTTGTTTAATGCCCATCTTGTGGCGGGTTTAACGCCGATTGAGGCCAACCGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAgg                                                                                            <  1:751589/250‑1 (MQ=255)
                                                                                                                                            taataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAAAGCCGGTTGGGGCCAAAGGTTATCTCGGTTTTTTTATCGACCGGCCGCTGGGGATGGAAggg                                                                                           >  1:704292/1‑250 (MQ=255)
                                                                                                                                              ataaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGGTATCTCGATTTTTTTTTCGCCCGACCGCTGGGGATGAAAgggtt                                                                                         >  1:451860/1‑248 (MQ=255)
                                                                                                                                                aaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGGTACTCGTTTAATGCCCATCTGGTTGCGGGTTTTACCCCGATTTAGGCCCACGGGTTTCTCCCTTTTTTTTTTCACCGACCGCTTGGGATGGAAGGTtatt                                                                                       >  2:324811/1‑250 (MQ=255)
                                                                                                                                                aaGCGGGGTTACCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAAAGCCGATTTAGGCCAACGGTTATCTCGATTTTTTTATCGACCGGCCGCTGGGAATGAAAGGGtat                                                                                        >  1:166861/1‑250 (MQ=255)
                                                                                                                                                          aCCGGTTTGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAAtct                                                                            <  1:745035/251‑1 (MQ=255)
                                                                                                                                                                           agaagaGCCAGGAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAAGTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAg                                                            <  2:707550/251‑1 (MQ=255)
                                                                                                                                                                             aagaGACAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGGTTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAggg                                                          >  1:799465/1‑251 (MQ=255)
                                                                                                                                                                                agCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACAATTCGCGGTCAGGGggt                                                       <  1:318253/251‑1 (MQ=255)
                                                                                                                                                                                 gCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGGATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGGTACTCGTTTAATGCCCCTCTGGTGGCGGGGTTTACGCCGCTTGAGGCCAAAGGTTATCTCTGTTTTTTTTTCGGCCCACCGCTGGGGATGAAAGGTTATTTTTTCAATTCTCCCATTTGCGGGCAggggggg                                                      >  2:62768/1‑249 (MQ=255)
                                                                                                                                                                                   cAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGGAAGGGTATATTTTCAATCTCACCATTCGCGGGCAggggggggg                                                    >  1:124726/1‑247 (MQ=255)
                                                                                                                                                                                    aGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGgtggtg                                                   <  2:451860/251‑1 (MQ=255)
                                                                                                                                                                                         aaGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATg                                                                                                                                                                                                                <  1:76390/94‑1 (MQ=255)
                                                                                                                                                                                          aGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTTAGGCCAACGGTTATCTCGATTTTTTTATCGACCGGCCGCTGGGGATGAAAGGGTATATTCTCCATCTACCCATTCGCGGTCAGGGGGTGGTGAAAAAt                                             >  1:317152/1‑251 (MQ=255)
                                                                                                                                                                                              gCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCa                                                             <  2:64796/231‑1 (MQ=255)
                                                                                                                                                                                              gCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGGGGGGAAAAATCAgg                                         >  1:329620/1‑251 (MQ=255)
                                                                                                                                                                                              gCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAAAGCCGGTTTGGGCCCAAGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCCCCCTTCGCgggca                                                             >  1:64796/1‑228 (MQ=255)
                                                                                                                                                                                                 gTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGAc                                      <  2:525064/251‑1 (MQ=255)
                                                                                                                                                                                                   gACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCCTCTGGTGGCGGGTTTAACCCCGGATGGAGCCAACGGTTATTTCGATTTTTTTATCACACTACCGCTGGGAGTGAAAGGTTATATTCTCAATCTCCCCATTCTCGGGCAgggggggggggaaaaacaggtacg                                     >  2:649635/1‑231 (MQ=255)
                                                                                                                                                                                                    aCGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTTAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGGATGAAAGGTTATATTCTCAATCTCACCATTCGCGGGCAGGGGGTGGGGAAAAATCAGGGACg                                     >  1:285767/1‑249 (MQ=255)
                                                                                                                                                                                                       gCAATGTCTTATGCAAGATGGACATTTGGTTTCTTCTCTGAATGGCGGGAGGTTGAAAATTATAGCTTAAGCGCAAAAAGAGCCCCTGCTTCCGGGATACTCTTTTAATGCCCATCTGGAGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATAGTCTCAATCGCAACATTCGCGCTCAGGGGGTGGTGAAAAATCAGGTACGAGAAt                                <  2:745034/251‑1 (MQ=255)
                                                                                                                                                                                                       gCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGGGGCGGGGTTAACGCCGGTTGAAGGCAACGGGTATCTCCGTTTTTTTTTCGCCCGACCGCTGGGGATGAAAGGTTTTATTCTCAATTCTCCCAATCGCGGTCCGGGGGGGGGGTAAAAAAAAGGACGAGAAt                                >  2:537717/1‑251 (MQ=255)
                                                                                                                                                                                                         aaTGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAAttt                              <  2:553346/251‑1 (MQ=255)
                                                                                                                                                                                                           tGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCCAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGGGGCGGGGTTTACGCCGATTGAGGGCAACGGTTATCTCGGTTTTTTTATCGACCGACCGCTGGGGATGAAAGGTTATATTCTCAATCTCAACATTCGCGGTCAGGGGGTGGTGAAAAAACAGGGGAGAGGATTTGt                            >  2:774079/1‑251 (MQ=255)
                                                                                                                                                                                                                gATGCAATATGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGTTTGAGGCCAACGGTTATTTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCTGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGcc                       >  1:224758/1‑251 (MQ=255)
                                                                                                                                                                                                                         tGGACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGGAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGGGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCggggg              >  1:241931/1‑249 (MQ=255)
                                                                                                                                                                                                                          agACAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGGTTAAAGCCGATTTAGGCCCACGGGTTTTTCGGTTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTTTCAATCTCACCATTCGCGGGCAGGGGGTGGGGAAAAATCAGGGACGAGAATTTTTTTGCCGACCgggggg             >  1:45825/2‑248 (MQ=255)
                                                                                                                                                                                                                            aCAATTGGTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGata           <  2:945614/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gTTTCTTCTCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGCAAGGGTATATTCTCAATCTCCCCATTCGCGGGCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGAAATTTTGCt    >  1:976612/1‑251 (MQ=255)
                                                                                                                                                                                                                                      tcttctCTGAATGGCGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAATGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGTGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTTTGCCGACCGGGTGATATTTTGCTGt  <  1:749800/250‑1 (MQ=255)
                                                                                                                                                                                                                                             tGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGaa                                 <  1:22106/212‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                
TCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAGAAGAGCCAGTAAAAGACGCAGTGACGGCAATGTCTGATGCAATATGGACAATTGGTTTCTTCTCTGAATGGTGGGAGTATGAAAAGTATGGCTGAAGCGCAAAATGATCCCCTGCTGCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGT  >  NC_000913/70128‑70606

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