Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,788,649 T→C N183N (AAT→AAC yibQ → divergent polysaccharide deacetylase domain‑containing protein YibQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,788,6490TC97.7% 145.8 / NA 44N183N (AAT→AACyibQdivergent polysaccharide deacetylase domain‑containing protein YibQ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/20);  minor base A (1/0);  total (24/20)
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.99e-01

GAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTATTCCTCGACGATTCGCAAAATGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGACCTGGCGCGTCGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCG  >  NC_000913/3788410‑3788885
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
gAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACAt                                                                                                                                                                                                                                    <  1:357701/251‑1 (MQ=255)
   cgcATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCg                                                                                                                                                                                                                                 <  1:16452/251‑1 (MQ=255)
      attattCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCgtgt                                                                                                                                                                                                                              <  1:152713/251‑1 (MQ=255)
            cGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAGCCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAAt                                                                                                                                                                                                                         <  1:52138/250‑1 (MQ=255)
            cGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTGACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAAtt                                                                                                                                                                                                                        <  1:275531/251‑1 (MQ=255)
             gTAGTGCGGTCAATAACGTGCACTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGGGGGCAAttt                                                                                                                                                                                                                       >  2:258657/1‑251 (MQ=255)
               agtgcgagcAATAACGTGCCCTATGCCGTGGGGATCAACAGCCACAGGTGTAGCAAGAGGACCTCTAACCTGTGTGGTATGCATAAAGTGATACAGGCACTAGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTGCTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGATCGAAGCGGACCTCCGTGTTCAATTTaa                                                                                                                                                                                                                     <  2:96977/243‑1 (MQ=255)
                gTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAAAGAAGAGGACATCAGTGTGCAATTTAAt                                                                                                                                                                                                                    >  1:364967/1‑251 (MQ=255)
                       cAATAATGTGCCCTATGCCGTGGGGATAAACAAACACATGGGTCGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAAGGCCCAGTCGATGCGCGCCGCGCAAGGTACTGGCGTTAAGGTGATTAATCGGAAGTTGTTCGTCGACTATTCGCAGAACGAAGCTGACAGCCGTGTGCAATTTAATCTCGCaa                                                                                                                                                                                                             <  2:437333/251‑1 (MQ=255)
                              gTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGatct                                                                                                                                                                                                      <  1:258657/251‑4 (MQ=255)
                                 cccTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGgc                                                                                                                                                                                                   >  2:216387/1‑251 (MQ=255)
                                                    acaacCACAAGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGCTCTGGCGCGGCGCAACGGTTCAAcc                                                                                                                                                                                >  2:132737/1‑251 (MQ=255)
                                                              gggtagcaAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAAccattgccattg                                                                                                                                                                      <  2:106921/247‑1 (MQ=255)
                                                                  agcaAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCa                                                                                                                                                                  <  2:45050/251‑1 (MQ=255)
                                                                     aaGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATcc                                                                                                                                                               >  2:277768/1‑251 (MQ=255)
                                                                                 aaCCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCTGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAAc                                                                                                                                                   <  1:250386/251‑1 (MQ=255)
                                                                                      gTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTTGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTAc                                                                                                                                              <  1:366287/251‑1 (MQ=255)
                                                                                        ttGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTAcg                                                                                                                                             <  2:55652/250‑1 (MQ=255)
                                                                                         tGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGAc                                                                                                                                                                                                                                      >  1:270639/1‑160 (MQ=255)
                                                                                         tGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGAc                                                                                                                                                                                                                                      <  2:270639/160‑1 (MQ=255)
                                                                                            tATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGGGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCgtgc                                                                                                                                        >  2:129777/1‑250 (MQ=255)
                                                                                               gCAGAAAGTGAGTCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCATAACCATCGGTAATACCCAGGCGAGGCGCGCCTCGCAAGGTACTGGTGTGAAGGTGATCAAACGGAAGGTGTTCCGCGACGATTCGCAGAACGAAGCGGACAACCGTGTGCAATTTAATCGCGCAATTGATCTGGCTCGACGCAACGGTTCAACCATTGCTATTGGGCATCTTCACCCTTCAACGGTACGCgtgcttc                                                                                                                                     <  2:290352/251‑5 (MQ=255)
                                                                                                   aaaGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCCATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCgtgctacaac                                                                                                                                  >  1:411961/1‑243 (MQ=255)
                                                                                                             aGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTtat                                                                                                                       >  1:151309/1‑251 (MQ=255)
                                                                                                                   tGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCCACGGTACGCGTGCTCCAACAGATGGTTTATAACCTg                                                                                                                 >  1:309026/1‑251 (MQ=255)
                                                                                                                         gTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCCTTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTgccgcc                                                                                                            >  2:319765/1‑250 (MQ=255)
                                                                                                                                   ttACTTCCTCGACAGCGTCACCATCGGAAATACCCAGGCGATGCGCGCCGCGCACGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCACCGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGc                                                                                                 >  2:51842/1‑251 (MQ=255)
                                                                                                                                   ttACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTa                                                                                                    <  1:82789/248‑1 (MQ=255)
                                                                                                                                      cTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGCTGCGCGCCGCGCAAGGTACTGGCGTGACGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGAGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTgg                                                                                              >  2:116614/1‑251 (MQ=255)
                                                                                                                                            cGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGaa                                                                                          >  2:57433/1‑249 (MQ=255)
                                                                                                                                             gACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCACGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGCCATTACGCTGGTGaacgc                                                                                       >  1:225001/1‑248 (MQ=255)
                                                                                                                                              aCAGCGTAACCATCGGTAATACCCAAGCGCTGCGCGCCGCGCACGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGCCGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGGACGCGTGCTACAACAGATGGTTTCTAACCTGCCGCCAGAAATTACGCTGGTGAAAGcc                                                                                      >  2:307353/1‑251 (MQ=255)
                                                                                                                                                         aTCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCt                                                                           >  1:148962/1‑251 (MQ=255)
                                                                                                                                                            ggTAATACCCAGGCGATGAGCGCCGCGCAAGGCACTGGCGTGACGGTGATCAAACGGAAGGTGTTCCTCGACGATTTGCCGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCACCCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTCCAACAGATGGTTTATAACCTGCCGCCCGACATTCCGCTGGTGAAAGCCAGCAGCTTGCtga                                                                         >  2:275752/1‑250 (MQ=255)
                                                                                                                                                              tAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCtgaatg                                                                      >  1:334818/1‑251 (MQ=255)
                                                                                                                                                                                      gcAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTacac                                             >  1:230097/1‑251 (MQ=255)
                                                                                                                                                                                                ggCGTGAAGGTGAACAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCACCCTCCCCCTTCAACGGTACGCGTGCTACAACAGAGGGTTTATAACCGGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCATGCTGAAGGAACCGCAAGGTTATACTTCTACACCGCCGaaaaa                                   <  2:285170/251‑1 (MQ=255)
                                                                                                                                                                                                 gCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAAc                                  <  1:216387/251‑1 (MQ=255)
                                                                                                                                                                                                    ggAAGGGGAACAAACGGAAGGTGTTCCCCGGCGAATCGCAGAACCAAAGGGACCTCCGTGTGGAAATTAAACGCGCAAATGAACTGGCGGGGCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCt                               <  1:260241/250‑1 (MQ=255)
                                                                                                                                                                                                            aTCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTg                        >  1:65446/1‑250 (MQ=255)
                                                                                                                                                                                                                  cGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCAcc                 <  1:57433/251‑1 (MQ=255)
                                                                                                                                                                                                                   ggAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCCGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACcg                >  1:163254/1‑251 (MQ=255)
                                                                                                                                                                                                                             ccTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACACCAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTCCACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgt      >  1:240472/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gacgaTTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTACGCGTGCTACAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCGGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAAtccgttccg  >  1:12818/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
GAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTATTCCTCGACGATTCGCAAAATGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGACCTGGCGCGTCGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACACCGCCGAAAAACGCTGTGCCTGATGCACCGCGTAATCCGTTCCG  >  NC_000913/3788410‑3788885

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