Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,634,882 C→T Q143Q (CAG→CAA rob ← DNA‑binding transcriptional dual regulator Rob

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,634,8820CT97.8% 132.2 / ‑3.9 46Q143Q (CAG→CAArobDNA‑binding transcriptional dual regulator Rob
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (19/25);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

CACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCC  >  NC_000913/4634636‑4635125
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
cacaTATTCGCCGCCCTGCAGCATCACCGGAGGCCCCGTCAGTCCATAGCCATCTGCCTGATCCTGGGCTAACTCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTgg                                                                                                                                                                                                                                                  <  2:114069/250‑1 (MQ=255)
 acaTATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATCGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGCACTGATAACGCATTTCATGGCGGAAATCAGGGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGt                                                                                                                                                                                                                                                >  1:64262/1‑251 (MQ=255)
     atTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCATTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAac                                                                                                                                                                                                                                             <  1:470586/250‑1 (MQ=255)
        cgccgcCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTTGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACAcc                                                                                                                                                                                                                                          <  1:259485/250‑1 (MQ=255)
          ccgccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCCTCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGCCTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAAACCGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAAAACCac                                                                                                                                                                                                                                        >  2:286515/1‑249 (MQ=255)
          ccgccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAAt                                                                                                                                                                                                                                       <  2:494739/251‑1 (MQ=255)
            gccCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGCAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATc                                                                                                                                                                                                                                      >  1:69808/1‑250 (MQ=255)
                     caTCACCGGATGCCCCGACAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGTTCCAGCGACAAGGAGTAGCTTTGGGTAacacacaacagcagcgtac                                                                                                                                                                                                                            >  1:267015/1‑236 (MQ=255)
                          ccGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTcc                                                                                                                                                                                                                       <  1:476748/251‑1 (MQ=255)
                            ggATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCCGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTcc                                                                                                                                                                                                                       >  1:492223/1‑249 (MQ=255)
                             gATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGGTAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAACGAACAAGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAgg                                                                                                                                                                                                                    >  1:351146/1‑251 (MQ=255)
                                        aGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCGAGAAATGAGAACGCATTGCAGGGCGGAAATCAAAAATTGGCTCCAGCGAACAGGAGTAGAATTGGGTAACACCAAGCAGCGGCGTATATGCCAGGGTGACaaagtt                                                                                                                                                                                                         >  1:335526/1‑248 (MQ=255)
                                        aGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAAttt                                                                                                                                                                                                         <  1:69820/251‑1 (MQ=255)
                                           aCATAGCCTTCTGCCGTATCCTGGGCTAACGCGGTGGTATAGAAGTCCTCTTGCTCGTCGTCTTTATCCTGTCTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACGGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTGTGGGTAACACCAATCAGCGGCTTATCTGCCAGGGTTACACATTtgtg                                                                                                                                                                                                      <  2:460878/251‑1 (MQ=255)
                                           aCATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCGTGCTCGTCGTCTTTATCCGGACTCGGACGCGTTTCATTCAGGACGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATGTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTtgtg                                                                                                                                                                                                      <  1:455667/251‑1 (MQ=255)
                                                gCCATCTGCCTGATCCTCGACTAACGCGGTAGTATAGAATACATCATGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGACGAGTAAATCGTGCCAGCACTAATAACGCCTTTCAGGGCGGGAATCAGAGATTTGGTCCAGCGCACAGGAGTAACTTTGGGTAACACCAATCAGCGGGGTATCTTCCTGGGAGGCAAAATTGTGATCTg                                                                                                                                                                                                 >  2:454985/1‑251 (MQ=255)
                                                   aTCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGc                                                                                                                                                                                               >  1:491028/1‑250 (MQ=255)
                                                     cTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACCCCAAACAACGGCGTATCTTCCCGGGTGACAAATTTGTGCTCTGGcct                                                                                                                                                                                             >  2:364138/1‑248 (MQ=255)
                                                       gCCTGATCCTGGGCTAACGCGGTGGTATAGAATTCCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCAt                                                                                                                                                                                             <  1:48268/248‑1 (MQ=255)
                                                          tGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGa                                                                                                                                                                                        >  1:185638/1‑250 (MQ=255)
                                                                 gggCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCAccc                                                                                                                                                                                <  1:497157/251‑1 (MQ=255)
                                                                 gggCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCAcc                                                                                                                                                                                 <  1:443417/250‑1 (MQ=255)
                                                                      aaCGCGGTGGTATAGCATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGCATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGTTAAAGCATTTCATGGCGGAAATCAGAGATTTGCTCCTACGAACAGGAGTAGCTTTTGGTAAAACCAAACAGCGGGGTTTCTTCCCGGGTGACAAATTTGTGCTATGGCATAGTGAAGTCAccccagcg                                                                                                                                                                           >  2:217581/1‑246 (MQ=255)
                                                                       aCGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGGGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGGAACACCAATCAGCGGCGTATCTTCCAGGGTGCCAAATTTTTGCCCTGGCCTAGTGAGTTCAccccagcgg                                                                                                                                                                          >  2:321309/1‑245 (MQ=255)
                                                                       aCGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCCGCGAACAGGAGGAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGAGAATTCACCCAGGCGt                                                                                                                                                                          >  1:266630/1‑251 (MQ=255)
                                                                            gtggtCTAGAGTACCTCGTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCAGTAGAGCACCGTCGGAATGGTCGGCGCGTTGACGAGCAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATAAGAGATTTGGTCCCACGAACAAGAGTAGCTTTGGGTAACACCAATCAGCGGCGGATATTCCAGGGTGACAAATTTTTGCTCTGGCATAGTGAATTCACCCAGGCgaggagg                                                                                                                                                                     >  2:454957/1‑245 (MQ=255)
                                                                               ttataGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATGCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCAc                                                                                                                                                                                  <  2:475517/234‑1 (MQ=255)
                                                                               gtATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCAc                                                                                                                                                                                  >  1:475517/1‑235 (MQ=255)
                                                                                   aGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGGCGCGTTTCATTCAGGCCGAAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGAAAACGCATTTCATGGAGGAAATCAGAGATTTGCTCCAGCAAAAAGGAGAAACTTTGGGTAAAACCAAACAGCGGCGTATCTTACAAGGTGAAAAATTTGTGCTCTAGCATAGTAGATTCACCCAGGCgacagcgcgggg                                                                                                                                                                >  2:446682/1‑238 (MQ=255)
                                                                                   aGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGCCAACGCATTTCATGGCGGCAATCCGAGCTTTGCTCCAGCGAACAGGGGTAGCTTTGGGTAACACCAATCAACGGCCTATCTTCCCAGGTGGCCAATTTTTGCTCTGGCATAGTGAATTCACCCAGGCggttcggcgggggg                                                                                                                                                              >  2:446677/1‑238 (MQ=255)
                                                                                       tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCa                                                                                                                                                                                                                      <  1:97499/191‑1 (MQ=255)
                                                                                       tACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGGAACACCAATCAGCGGGGTATCTTCCa                                                                                                                                                                                                                      >  2:97499/1‑191 (MQ=255)
                                                                                        cccTCTTGCGTGTCGTCTTTTGCCTTACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATAcc                                                                                                                                                         <  2:454163/250‑1 (MQ=255)
                                                                                               gCTCGTCGTCTTTATCCTGACTCGTCCGCGTTTCAATCAGGCCGTAGAGCACCGGCGGAATGGTCGTCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAAGCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAAGTCACCCAGGCGGAGCGGCGGGCGAATACCAAAGgcg                                                                                                                                                  <  2:492223/251‑1 (MQ=255)
                                                                                                            aTCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCAt                                                                                                                                            <  1:432990/244‑1 (MQ=255)
                                                                                                                         cgcgTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATGTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAg                                                                                                                        <  1:341967/251‑1 (MQ=255)
                                                                                                                                    caGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTGGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCt                                                                                                             <  2:484107/251‑1 (MQ=255)
                                                                                                                                            agagCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCaaa                                                                                                      <  2:391019/250‑1 (MQ=255)
                                                                                                                                                    ggcggAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGAGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGcttctt                                                                                             <  1:402601/251‑1 (MQ=255)
                                                                                                                                                                         ccGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCGTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAAGTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACGGCTTTTTGAATGCGCGGGTAAATGTCTg                                                                        <  2:187933/251‑1 (MQ=255)
                                                                                                                                                                                      gCCAGAACTGATAGAGCATTTCATGGCTGAAATCAGAGATTGGCTCCAGCGAACAGTAGTAGGTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATGTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCGCCATTCAGGAGAACGGCGGGAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGaa                                                            <  2:498157/250‑1 (MQ=255)
                                                                                                                                                                                      gCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAAGTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAg                                                           <  1:46668/251‑1 (MQ=255)
                                                                                                                                                                                         aGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCgg                                                        <  1:30875/251‑1 (MQ=255)
                                                                                                                                                                                                              ggCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGa                                   <  2:17788/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAAcg                                                                                                                                 <  2:50084/135‑1 (MQ=255)
                                                                                                                                                                                                                                               aGTAGCTTTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGAccgcc  <  1:217581/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   
CACATATTCGCCGCCCTGCAGCATCACCGGATGCCCCGTCAGTACATAGCCATCTGCCTGATCCTGGGCTAACGCGGTGGTATAGAATACCTCTTGCTCGTCGTCTTTATCCTGACTCGGACGCGTTTCATTCAGGCCGTAGAGCACCGGCGGAATGGTCGGCGCGTTGCCGAGAAAATCGTGCCAGAACTGATAACGCATTTCATGGCGGAAATCAGAGATTTGCTCCAGCGAACAGGAGTAGCTCTGGGTAACACCAATCAGCGGCGTATCTTCCAGGGTGACAAATTTGTGCTCTGGCATAGTGAATTCACCCAGGCGTAGCGGCGGGCGAATACCAAAGGCGCTCCATTCAGGAGAACGGCGGTAAAGTGCAGGAGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCC  >  NC_000913/4634636‑4635125

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