Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,311,337 C→T A236T (GCA→ACA)  alsB ← D‑allose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,311,3370CT98.5% 234.3 / NA 68A236T (GCA→ACA) alsBD‑allose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (37/30);  minor base G (0/1);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGCGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTCACCACCTT  >  NC_000913/4311093‑4311579
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
gggatcatTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGc                                                                                                                                                                                                                                              <  1:102294/248‑1 (MQ=255)
     catTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTg                                                                                                                                                                                                                                           >  1:544794/1‑249 (MQ=255)
                  tGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTGGTGACTGCCTGAGCAACACCCATTg                                                                                                                                                                                                                             <  2:292350/250‑1 (MQ=255)
                       gACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCg                                                                                                                                                                                                                       <  1:282115/251‑1 (MQ=255)
                             gATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCAt                                                                                                                                                                                                                 <  1:106784/251‑1 (MQ=255)
                             gATTGAATCGACCACTATAAATTCAGGTGCTTTCTCCACCGGCATAACATTGCCGGATTTCTCAGCGTCAACCCTCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTACGGCGTTTGTGACTGCCTGAGCCACACCCATTGCCATAGTGTcct                                                                                                                                                                                                                 >  2:402815/1‑249 (MQ=255)
                              aTTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCAtt                                                                                                                                                                                                                <  1:323839/251‑1 (MQ=255)
                                tGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTc                                                                                                                                                                                                               <  1:574834/250‑1 (MQ=255)
                                 gAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCAGCGTGGCATTcg                                                                                                                                                                                                              <  2:511756/250‑1 (MQ=255)
                                      gACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAAt                                                                                                                                                                                                        <  1:109471/251‑1 (MQ=255)
                                      gACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAAt                                                                                                                                                                                                        >  1:277803/1‑251 (MQ=255)
                                       aCCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATa                                                                                                                                                                                                       <  1:445858/251‑1 (MQ=255)
                                        ccAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACTCCCATTGCCATCGTGTCATTCGCGCAATAg                                                                                                                                                                                                      <  1:183741/251‑1 (MQ=255)
                                        ccAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAAt                                                                                                                                                                                                        >  2:492823/1‑249 (MQ=255)
                                                aaTTCCGGTGCTTTATCCAGCGGTATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGTGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTTCCTGAGCAACACCCATTGCCATCGTGTCAGTCGCGCAATAGATCGCttt                                                                                                                                                                                              <  1:234471/251‑1 (MQ=255)
                                                  ttcaggTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACACGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCAGGCGTTTGTGACTGCATGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTa                                                                                                                                                                                             >  2:465138/5‑250 (MQ=255)
                                                                      ggATAACCTTGCCGGATTTCTCAGAGTCACCCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGCATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAac                                                                                                                                                                        >  2:116554/1‑251 (MQ=255)
                                                                       gATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTcg                                                                                                                                                                                                              >  1:258091/1‑212 (MQ=255)
                                                                       gATAACCTTGCCGGAGTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGGGCCGATATCCGCCGGGTTCTGGGCAAACGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCTGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTcg                                                                                                                                                                                                              <  2:258091/212‑1 (MQ=255)
                                                                           aCCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGtt                                                                                                                                                                   >  1:403845/1‑251 (MQ=255)
                                                                             cTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAacac                                                                                                                                                                      >  2:221752/1‑246 (MQ=255)
                                                                             cTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAacac                                                                                                                                                                      <  1:221752/246‑1 (MQ=255)
                                                                              ttGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAgg                                                                                                                                                                >  1:30620/1‑250 (MQ=255)
                                                                               tGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTAATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCCTTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTg                                                                                                                                                               >  2:363104/1‑251 (MQ=255)
                                                                                        tctcAGCGTCAACCAACAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCTCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGGTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCAGTCGCGCAATAGATCGCTGTAATATTCGGATTACGTTGCAACACGTTAGTGGCGTCATc                                                                                                                                                       <  2:494820/250‑1 (MQ=255)
                                                                                            aGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGCTATCCGCCGGGATCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATCTTCGGATTACGTTGCAACACGTTAGCGGCGACATCCatgg                                                                                                                                                  >  1:228374/1‑248 (MQ=255)
                                                                                                  aaCCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCCTTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTa                                                                                                                                             >  1:570091/1‑250 (MQ=255)
                                                                                                    ccctCAGCGTCAGGCCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATGTGTCCGGCTTCCTCCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATg                                                                                                                                          <  1:24525/248‑1 (MQ=255)
                                                                                                    ccATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATg                                                                                                                                          >  1:191278/1‑251 (MQ=255)
                                                                                                     cATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGc                                                                                                                                         >  1:287025/1‑251 (MQ=255)
                                                                                                        cAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCACAACGTTAGTGGCGACATCCAGTGCTTTAATGCGGt                                                                                                                                      >  2:219660/1‑251 (MQ=255)
                                                                                                         aGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTc                                                                                                                                     >  1:500767/1‑251 (MQ=255)
                                                                                                            ttCAGCCCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCCCCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAata                                                                                                                                                                                          <  2:415594/195‑1 (MQ=255)
                                                                                                            ttCAGACCCGATGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAata                                                                                                                                                                                          >  1:415594/1‑195 (MQ=255)
                                                                                                                gACCCGTTGCGCCGCTATCCCCCGGGTTCTGCGCAACCGTCGCGGACATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTACTATTCGGATTACGTTGCAACACGTTAGTGGCGCCATCCAGTGCTTTAATGCGGTCCCAGTCg                                                                                                                              >  1:482255/1‑251 (MQ=255)
                                                                                                                  cccGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGc                                                                                                                            >  1:128733/1‑251 (MQ=255)
                                                                                                                   ccGTTGCGCCGATATCCTCCGGGTTCTGGGTATCCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCGGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCATCACTTTAGTGGCGACATCCAGTGCTTTAATGCGGACCCAGTCGGcg                                                                                                                           <  2:536991/251‑2 (MQ=255)
                                                                                                                   ccGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCATCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCACCACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGGGGCGCAATCCAGTGCTTTAATGCGGTCCCAGTCGGc                                                                                                                            >  2:48455/1‑250 (MQ=255)
                                                                                                                    cGTTGCGCCGTTATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGCGCAACACCCATTGCCCTCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGTTTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTCAATGCGGTCCCAGTCGGc                                                                                                                            >  1:404702/1‑249 (MQ=255)
                                                                                                                        gcgcCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAggctg                                                                                                                      <  1:174431/251‑1 (MQ=255)
                                                                                                                          gcCGATATCCGCCGAGTTCAGGGCAACCGTACCGGTCCTTTGACCGGCTTCCACAATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCCTTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATCTTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGTCAggctggc                                                                                                                    >  2:386467/1‑251 (MQ=255)
                                                                                                                             gATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGctggctgg                                                                                                                 >  1:295634/1‑251 (MQ=255)
                                                                                                                                         gTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTg                                                                                                      >  1:216200/1‑250 (MQ=255)
                                                                                                                                               ggCAACCGTCGCGGCAATTTGTCCGGCTTCCCCCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCCTTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGgc                                                                                               >  2:76694/1‑251 (MQ=255)
                                                                                                                                                      gTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCt                                                                                                                                                                                                >  1:255189/1‑147 (MQ=255)
                                                                                                                                                      gTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAAGATTCGGATTACGTTGCAACACGTTAGt                                                                                                                                                                <  2:482800/179‑1 (MQ=255)
                                                                                                                                                      gTCGCGGTCATTTGTCCGGCTTCCACCATTTGGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCAGTCGCGCAATAGATCGCt                                                                                                                                                                                                <  2:255189/147‑1 (MQ=255)
                                                                                                                                                      gTCGCGGTCATTTGTCCGGATTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGt                                                                                                                                                                >  1:482800/1‑179 (MQ=255)
                                                                                                                                                          cgGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTg                                                                                    >  1:204692/1‑251 (MQ=255)
                                                                                                                                                                   gTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTTCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGAGGCGACATCCAGTGCTTTAATGCGATCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCg                                                                           >  2:160138/1‑251 (MQ=255)
                                                                                                                                                                   gTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCg                                                                           >  2:114146/1‑251 (MQ=255)
                                                                                                                                                                            ccaccaTTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCGTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCa                                                                  <  2:102603/251‑1 (MQ=255)
                                                                                                                                                                                    tttcGGGCTTCCGGAATGCCACCTGTACCGACGACCTGCTCTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCGTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTGTAATGCTGTACCAGGCTGCGGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTacgac                                                           <  2:570091/247‑1 (MQ=255)
                                                                                                                                                                                         ggCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCtt                                                     <  1:253664/251‑1 (MQ=255)
                                                                                                                                                                                          gcgtcaggAATGCCATCTGTACCGACGACCAGCACTGTTCCCTGTTAGCCGGCGTGTGGGACTGCCTGAGCAACCCCCATTGCCATCGTGTCATTCGCTCAATAGATCGCTTTAATATTCTGATTACGTTGCAACGCGTTAGTGGCGTCATCCAGTGCGTTAATGCGTTCCCAGTCGGCAGGGTGGCTGGCTACAAGCTTGATCTGGCTTGGTTTTCTGAAGGCTTCGGTTTTACCATTACGACGCGCTTc                                                    <  2:223209/245‑1 (MQ=255)
                                                                                                                                                                                                  aaTGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAgg                                            <  1:114146/251‑1 (MQ=255)
                                                                                                                                                                                                      ccATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGtt                                        <  2:204692/251‑1 (MQ=255)
                                                                                                                                                                                                            gTACCGACGACCAGCACTTTTCCCGTTTGTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGGCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTAcc                                     <  1:358100/248‑1 (MQ=255)
                                                                                                                                                                                                              aCCGCCGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTTACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCGTCGGTGGCTCCATTACGACGCGCTTCACCGGAGGCGTTACCGGCt                                 <  2:294602/250‑1 (MQ=255)
                                                                                                                                                                                                              aCCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCtt                                <  2:287025/251‑1 (MQ=255)
                                                                                                                                                                                                                 gacgacCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAAGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTAc                             >  2:300118/1‑251 (MQ=255)
                                                                                                                                                                                                                   cgacCAGCACTTTTCCCGTTTTTCAGGCGTTTGTGCATGCTTGAGCACCCCCCCTTGCCATCGTGTCATTCGAGCAATAGATCGCTTTAATATTCGGATTCCGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTAccc                           >  2:284862/1‑251 (MQ=255)
                                                                                                                                                                                                                    gacCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCt                          >  1:51041/1‑251 (MQ=255)
                                                                                                                                                                                                                     acCAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGCGGCGTTACCGGCTTTACCCTc                         >  2:167718/1‑251 (MQ=255)
                                                                                                                                                                                                                       cAGCACTTTTCCCGTTTTTCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCa                        >  2:228728/1‑250 (MQ=255)
                                                                                                                                                                                                                                tCCCTTTTTTCCGTCGTTTGTGACTGCCTGAGCAACACTCGTGGCCATCGTGTCATTCGCTCAATCGATCGCTTTTATATTCGGATTACGTTGCGACACGTTAGTGGCGACATCCCGTCCTTTAATGCGGTCCCAGTCGGCAGGCTGGGTGGCTACAAGCTTGATCTGGCTTGCTTTCTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTTCCGGCTTTACCCTCAATGATTGCg               <  2:24906/250‑1 (MQ=255)
                                                                                                                                                                                                                                       tttCCGGCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTcac        <  1:116554/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ggCGTTTGTGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTCACCACCtt  >  1:196489/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
GGCATCATTCAGTTTTATTGAGTGACCAGGATTGAATCGACCAGTTTAAATTCCGGTGCTTTATCCAGCGGGATAACCTTGCCGGATTTCTCAGCGTCAACCATCAGCTTCAGACCCGTTGCGCCGATATCCGCCGGGTTCTGGGCAACCGTCGCGGTCATTTGTCCGGCTTCCACCATTTTGCGGGCTTCCGGAATGCCATCTGTACCGACGACCAGCACTTTTCCCGTTTTTCCGGCGTTTGCGACTGCCTGAGCAACACCCATTGCCATCGTGTCATTCGCGCAATAGATCGCTTTAATATTCGGATTACGTTGCAACACGTTAGTGGCGACATCCAGTGCTTTAATGCGGTCCCAGTCGGCAGGCTGGCTGGCGACAAGCTTGATCTGGCTTGCTTTTTTGAAGGCTTCGGTGGCACCATTACGACGCGCTTCACCGGAGGCGTTACCGGCTTTACCCTCAATGATTGCGACTTCACCACCTT  >  NC_000913/4311093‑4311579

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: