Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 71,170 T→C S262P (TCG→CCG)  araC → DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,1700TC94.5% 152.1 / ‑1.6 55S262P (TCG→CCG) araCDNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (30/22);  total (32/23)
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.35e-01

ACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTG  >  NC_000913/70940‑71410
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
aCATCAGCGATCACCTGGCAGACAGCAATTTTGATTTTCCCTGCGTCTCCCAGCATGTTTGATTGTCGCCGTCGCGTCTTTCCCATCTTATCCGCCAGCAGTGAGGGATTAGCGTCTGTAGCTGGCGCTAGGACAATCTCATTTGTCAGGCGAAGCTCCTTTTGAGCACTACCCTGATGCCTATCGCCACCGTCGTTCGCAAAGATGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATg                                                                                                                                                                                                                               <  2:66302/250‑1 (MQ=255)
  aTCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAACGTCTTAAGCTGGCGCGAGGACCAACGCCTTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGGTCAACTCTATTTCCCGCGCGTATTTAAAAAATGCAc                                                                                                                                                                                                                            >  2:449235/1‑251 (MQ=255)
         aTCAACTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCGGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCACCGCTTTAGTCAGGAGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGcc                                                                                                                                                                                                                     >  2:63486/1‑251 (MQ=255)
            aCCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGGGTATTTAAAAAATGCACCGGGGcccgc                                                                                                                                                                                                                  >  2:113048/1‑248 (MQ=255)
                     aCAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcgag                                                                                                                                                                                                         <  1:376450/251‑1 (MQ=255)
                     aCAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgagcg                                                                                                                                                                                                           >  1:650206/1‑249 (MQ=255)
                         cAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTc                                                                                                                                                                                                     <  1:601565/251‑1 (MQ=255)
                                  tatCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGTTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCAGGATGCCTATCGCCCCCGTCGGTCGCAATGTTGGCTTAAACTATCAAATCTATTTCCCGCGAATATTTAAAAAATGCAACGGGGCCAGCCCGGGCGCGTTTCGTGcggggt                                                                                                                                                                                            >  2:186674/1‑246 (MQ=255)
                                   atCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGAACTGTTGGTTTTGCAGACCAACTCTCTTTCCCGAGGGGATTGAAAAACTGACCCGGGGACAGCCCGACCGAGTTTAGTGCCGGTtg                                                                                                                                                                                           >  2:182516/1‑251 (MQ=255)
                                        cAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTga                                                                                                                                                                                        >  1:232920/1‑249 (MQ=255)
                                           cGTCGCACAGCATTTTTGCTTGTCGACGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGTGAGGACCAACGCATTAGGCAGGCGAAGCTGCTTTTTAGCCATACACGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGATTTGACGATCAAATCTATTTCCCGCGAGGATTTACAAAATGCCACGGGGCCAGCCCTAGCGAGTTTCGTGACGGTTGTgaggaa                                                                                                                                                                                    >  2:527924/1‑250 (MQ=255)
                                             tCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGCGCGAGTTTCGTGCCGGTTGTGAAGaaaa                                                                                                                                                                                  >  2:501893/1‑250 (MQ=255)
                                                cacaGCCTGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGCTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGc                                                                                                                                                                                                                      >  1:93428/1‑211 (MQ=255)
                                                cacaGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTg                                                                                                                                                                              >  1:21702/1‑251 (MQ=255)
                                                   aGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAat                                                                                                                                                                           >  2:136248/1‑251 (MQ=255)
                                                         tttGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGaaaaa                                                                                                                                                                                 <  1:363942/239‑1 (MQ=255)
                                                         tttGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGGTGTGAAGaaaaa                                                                                                                                                                                 >  2:363942/1‑239 (MQ=255)
                                                          ttGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGTAGAAAAAGAGAATGATGt                                                                                                                                                                      >  1:537762/1‑249 (MQ=255)
                                                          ttGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTa                                                                                                                                                                     <  1:87326/250‑1 (MQ=255)
                                                            gCTTTTCGCCGTCGCGGCTGGCACATCTTTACCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGTATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGc                                                                                                                                                                   <  2:130696/250‑1 (MQ=255)
                                                                      gTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGt                                                                                                                                                        >  1:678381/1‑251 (MQ=255)
                                                                          cgTCTGTCACATCTGTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATa                                                                                                                                                    <  1:70152/251‑1 (MQ=255)
                                                                                       ttttCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAAt                                                                                                                                       >  1:529467/1‑251 (MQ=255)
                                                                                               cagcagTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAAt                                                                                                                               >  1:155294/1‑251 (MQ=255)
                                                                                                       atggATTAGTTTCTTTAGCTGGCGCGCGAGTCAACGCATGAGTCAGGCCAAGTTGCTCTTGAGCACTACCCGGAAGCCTATCGCAACCGTCGGTTGAAATGTTGGTTTTGACGATCAACCCTATTTCCCGCGAGTAGTTAAAAAATGCACCGTGGCCAGCCCGAGCGAGTTGCGTTCCTGTTGTGAAGAAAAAGTGAATGATGTAGCCGCCAAGGTGTCATAATTGGTAACGAATCAGACAATTGACGGCt                                                                                                                       <  2:151733/249‑1 (MQ=255)
                                                                                                          taTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAACTGCACCGGGGCCAGCCCGAGCGCGTTTCGTTCCGGTTGTGAAGAAAAAGTGACTGATGTAGCCGTCCAGTTGTCATAATTGGTAACGAATCAgcaaatttcaggcttg                                                                                                                     >  2:595249/2‑235 (MQ=255)
                                                                                                            ttAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACg                                                                                                                  >  1:267813/1‑251 (MQ=255)
                                                                                                                        gCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAgg                                                                                                      <  1:285517/251‑1 (MQ=255)
                                                                                                                           ggCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACGACCCGGATGCATATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCTGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGGAGCCGTCAAGTTTTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAggg                                                                                                     <  2:657461/249‑1 (MQ=255)
                                                                                                                             agcgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTg                                                                                                 <  1:135833/250‑1 (MQ=255)
                                                                                                                              gcgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGc                                                                                                <  1:557316/251‑1 (MQ=255)
                                                                                                                               cgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCCATGTTGGTTTTGACGATCAACTCTATTTCCCGCGCGTATTTAAAAACTGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGGGAAGAAAAAGGGAATGATGTAGCCGTCAAGTTGTCCAAATTGGTAAAGAAACAGACacatg                                                                                                                             >  2:216722/1‑217 (MQ=255)
                                                                                                                               cgAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTg                                                                                                                             <  1:216722/221‑1 (MQ=255)
                                                                                                                                     ccAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGGGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGCCAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGACTcc                                                                                         >  2:500393/1‑251 (MQ=255)
                                                                                                                                       aaCGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTTTCATAATTGGTAACGACTCAGACAATTGACGGCTTGACGTAGAAGCATAGGGTTTGCAGAATCCCt                                                                                       >  1:34925/1‑251 (MQ=255)
                                                                                                                                        aCGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTg                                                                                      <  1:417994/251‑1 (MQ=255)
                                                                                                                                          gCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCt                                                                                    <  1:63486/251‑1 (MQ=255)
                                                                                                                                          gCATTAGCCAGGCGAATCTGCTTTTGAGCACTCCCCGGATGCCTATCGCCACCGGCGGTCGCAATGTTGGATTTGACAATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCTGGTTTCGTGCCGGTTGTGAAGAACAAGTGAATGATGTAGCGGTCAAGTTGTCCTATTTGGTAACGAATCAGCCAATTGACGGCTTGAAGGAGTAGCAGAGGGTTTgaaaaattcatgct                                                                                    >  2:529905/1‑238 (MQ=255)
                                                                                                                                                gTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTcc                                                                              <  1:136248/251‑1 (MQ=255)
                                                                                                                                                      cGAAGCTGCTTTCGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTg                                                                         >  1:218901/1‑250 (MQ=255)
                                                                                                                                                            tgctTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACa                                                                                                                                 <  1:470654/188‑1 (MQ=255)
                                                                                                                                                            tgctTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACg                                                                                                                          <  1:966/195‑1 (MQ=255)
                                                                                                                                                            tgctTTTGAGCACTACCCGGATGCCTATCGCCACAGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACa                                                                                                                                 >  2:470654/1‑188 (MQ=255)
                                                                                                                                                                  tGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGTGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATACCTGCTTCGTCCATTTGACAGGCACAttc                                                             >  2:42607/1‑249 (MQ=255)
                                                                                                                                                                      cACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCAttt                                                                          >  1:40352/1‑233 (MQ=255)
                                                                                                                                                                      cACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTAGTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCAttt                                                                          <  2:40352/233‑1 (MQ=255)
                                                                                                                                                                       aCTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACACTTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCa                                                         >  2:92053/1‑249 (MQ=255)
                                                                                                                                                                          aCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGTAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGAGGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCCAGCAt                                                    >  2:359402/1‑251 (MQ=255)
                                                                                                                                                                            ccGGATGCCTATCGCCACCGTCTGTCGCAATGTTGGTTTTGACTATCATCTCTATTTCCCGCGAGTTTTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGGTTCTGGCCTGTTGGGAATAAAAAGTGAATGATGTAGCCGTCAAGTTGTAATAATTGGTAACGAATCAGACACTTGACCGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCTTTTGACAGGCACATTATGCAAGCATTg                                                  <  2:610878/251‑1 (MQ=255)
                                                                                                                                                                             cGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGc                                                 >  1:513337/1‑251 (MQ=255)
                                                                                                                                                                                              cGTCGGGCGCAATGGTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAAGGCGCCGGGGCCAGCCCGAGCGAGTTGCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCGTGACGGAGTAGCATAGGGTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTGccc                               <  2:238956/251‑1 (MQ=255)
                                                                                                                                                                                               gacggtcgCAATGTTGGTTTTGACGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCAGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTAccc                               >  1:491782/3‑251 (MQ=255)
                                                                                                                                                                                                             tGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTa                  <  1:92053/250‑1 (MQ=255)
                                                                                                                                                                                                                      cGATCAACTCTATTTCCCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATg        >  1:502078/1‑251 (MQ=255)
                                                                                                                                                                                                                            acgctATTTCCCGCGAGTATTTAAAAAAAGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCggtg  <  2:502078/248‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
ACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTG  >  NC_000913/70940‑71410

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