Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,788,598 C→T G166G (GGC→GGT 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,5980CT98.1% 169.8 / NA 53G166G (GGC→GGTyibQdivergent polysaccharide deacetylase domain‑containing protein YibQ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/25);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCGCTGGAGAAAAATACGCTACGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTATTCCTCGACGATTCGCAAAATGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGACCTGGCGCGTCGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAC  >  NC_000913/3788359‑3788842
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
cAGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGa                                                                                                                                                                                                                                            >  1:397569/1‑250 (MQ=255)
 aGCCGCTGGAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGa                                                                                                                                                                                                                                            >  2:171257/1‑249 (MQ=255)
        ggAGAAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGAGCACCAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCGTTACTTCCTCGACAGCGGAACCATCGGTAATACCCAGGCGATGAGCGCCGCGCTAGGTACTGGCGTGaattggatca                                                                                                                                                                                                                                   <  2:385750/251‑9 (MQ=255)
          agaAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCaa                                                                                                                                                                                                                                  >  1:55652/1‑250 (MQ=255)
           gaAAAATACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAAc                                                                                                                                                                                                                                <  2:335228/251‑1 (MQ=255)
             aaaaTACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACgg                                                                                                                                                                                                                              <  2:63326/251‑1 (MQ=255)
                aTACGCTGCGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAg                                                                                                                                                                                                                           <  1:274655/251‑1 (MQ=255)
                aTACGCTGCGCACGGCGATGAGCAGCGACGAAATTGCGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGACGGTGATCAAACGGa                                                                                                                                                                                                                             >  2:275531/1‑249 (MQ=255)
                   cgctgcgcCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGtg                                                                                                                                                                                                                        >  2:82789/1‑250 (MQ=255)
                     ctgcgcCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGGGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTACCCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTAATGGCGTGAAGGTGATCACACGGCAggggtt                                                                                                                                                                                                                      >  2:11314/4‑247 (MQ=255)
                        cgcCCGGCGATGAGAAGCGACGAAATTGAGCTCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGGGGGGATCAACAACCACAGTGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATGCCCAGGCGATGCGCGCCGCGCAAGGTACGGGCGTGAAGGTTATCAACCGGAAGTTGTTCCt                                                                                                                                                                                                                   <  2:210554/251‑1 (MQ=255)
                        cgcCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCt                                                                                                                                                                                                                   >  2:312149/1‑251 (MQ=255)
                                 attaggagcGGCGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGGGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTc                                                                                                                                                                                                          <  2:384996/245‑1 (MQ=255)
                                   gAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGc                                                                                                                                                                                                        <  1:312149/251‑1 (MQ=255)
                                           acgaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCagaacgaa                                                                                                                                                                                                <  1:171257/251‑8 (MQ=255)
                                                   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)
                                                                                                                             gACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGAtt                                                                                                                                                                                                           <  1:149463/158‑1 (MQ=255)
                                                                                                                             gACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGAtt                                                                                                                                                                                                           >  2:149463/1‑158 (MQ=255)
                                                                                                                                    aaCCTGTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCTGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAAc                                                                                                       <  1:250386/251‑1 (MQ=255)
                                                                                                                                         gTTTGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTTGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTAc                                                                                                  <  1:366287/251‑1 (MQ=255)
                                                                                                                                           ttGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGATCTGGCGCGTCGCAACGGTTCAACCATTGCCATTGGGCATCCTCACCCTTCAACGGTAcg                                                                                                 <  2:55652/250‑1 (MQ=255)
                                                                                                                                            tGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGAc                                                                                                                                                                                          <  2:270639/160‑1 (MQ=255)
                                                                                                                                            tGGTATGCAGAAAGTGATGCAGGCACTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGTACTGGCGTGAAGGTGATCAAACGGAAGGTGTTCCTCGACGATTCGCAGAACGAAGCGGAc                                                                                                                                                                                          >  1:270639/1‑160 (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)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
CAGCCGCTGGAGAAAAATACGCTACGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTAACCATCGGTAATACCCAGGCGATGCGCGCCGCGCAAGGCACTGGCGTGAAGGTGATCAAACGGAAGGTATTCCTCGACGATTCGCAAAATGAAGCGGACATCCGTGTGCAATTTAATCGCGCAATTGACCTGGCGCGTCGCAACGGTTCGACCATTGCCATTGGGCATCCTCACCCTGCAACGGTACGCGTGTTGCAACAGATGGTTTATAACCTGCCGCCAGACATTACGCTGGTGAAAGCCAGCAGCTTGCTGAATGAACCGCAGGTTGATACTTCTACAC  >  NC_000913/3788359‑3788842

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: