Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,158,519 T→G 6.8% P217P (CCT→CCG ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,158,5190TG 6.8% 188.5 / 2.6 69P217P (CCT→CCGptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base T (38/25);  new base G (2/3);  total (41/28)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGC  >  NC_000913/1158371‑1158668
                                                                                                                                                    |                                                                                                                                                     
tGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCtt                                                                                                                                                      <  2:232226/150‑1 (MQ=255)
tGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCtt                                                                                                                                                      >  2:49017/1‑150 (MQ=255)
  tCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTc                                                                                                                                                    <  1:15344/150‑1 (MQ=255)
     ttCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAg                                                                                                                                                 >  1:858063/1‑150 (MQ=255)
       cATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTGTGGGATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGAt                                                                                                                                               >  1:773548/1‑150 (MQ=255)
         tttGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGc                                                                                                                                             <  2:832791/150‑1 (MQ=255)
         tttGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAAAGGTACGTTCCAGATgt                                                                                                                                             <  1:361702/150‑2 (MQ=255)
           tGGCCGCCGATTGGTTCTGCAATCCATACCTTCTCTCAGGGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAAAGGTACGTTCCAGATgtgg                                                                                                                                           <  1:1191576/150‑4 (MQ=255)
           tGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAAAGGTACGTTCCAGATgtgg                                                                                                                                           <  1:1223815/150‑4 (MQ=255)
            ggCCGCCGATTGGTTCTGCAATCCAGCCCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGGAGTTGCGGTTGGCATTTACGGTTTCATCGAACGTGGCGTGGTACCGTTTGGTCTGCACCACATCTGGGCCGTCCCTTTCCAGATGCAgg                                                                                                                                          >  2:81138/1‑149 (MQ=255)
               cgccgATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTg                                                                                                                                       >  1:312717/1‑150 (MQ=255)
                gccgATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTgg                                                                                                                                      <  1:713878/150‑1 (MQ=255)
                   gATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGa                                                                                                                                   <  2:265962/150‑1 (MQ=255)
                      tGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATa                                                                                                                                <  2:790704/150‑1 (MQ=255)
                          tCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACAcc                                                                                                                            >  1:406839/1‑150 (MQ=255)
                           cTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCa                                                                                                                           >  2:1010880/1‑150 (MQ=255)
                           cTGCAAACCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGACCCTTGCCTGGTACCGGTTTGGATGCACCACATCTGGAACGTACCATTCCAGATGCAGATTGGTGAATACAccc                                                                                                                           >  2:352207/1‑149 (MQ=255)
                            tGCAATCCAGACCTTCTATCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTATGGCATTTACGGTTTCATCGAACGTTGCATGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCaa                                                                                                                          <  1:219537/150‑1 (MQ=255)
                               aaTCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACgc                                                                                                                       >  1:1044977/1‑150 (MQ=255)
                                aTCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGTATTCACCaatctg                                                                                                                      <  2:958966/150‑5 (MQ=255)
                                   cAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACgcagca                                                                                                                   <  1:1012386/150‑1 (MQ=255)
                                       ccTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGcaggtc                                                                                                               >  2:60555/1‑150 (MQ=255)
                                          tctctcAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGcaggtcagg                                                                                                            <  2:521411/150‑1 (MQ=255)
                                             ctcAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGATCGTTGCCGGGTACCGTTTGGTCTGCACCACATCTGGAACGTCCCGTTAACGATGCAGAGTGGTGAATACACCAACGCAGcaggtctggtgt                                                                                                         >  1:207027/1‑148 (MQ=255)
                                             ctcAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCTGATGCAGATTGGTGAATACACCACCGCAGCAGGTCAGGttt                                                                                                         >  1:1198345/1‑150 (MQ=255)
                                             ctcAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGttt                                                                                                         >  2:1297080/1‑150 (MQ=255)
                                                aGTGGGCTGCTTACCAGACCCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGAAGGTCAGGTTTTcc                                                                                                      >  2:1069438/1‑150 (MQ=255)
                                                    ggCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACgg                                                                                                  >  1:799512/1‑150 (MQ=255)
                                                     gctgctTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGc                                                                                                 >  2:398411/1‑150 (MQ=255)
                                                      ctgctTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACGAACGCAGCAGGTCAGGTTTTCCACGGCg                                                                                                >  1:323897/1‑150 (MQ=255)
                                                      ctgctTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCg                                                                                                <  2:312717/150‑1 (MQ=255)
                                                        gctTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGAc                                                                                              <  2:994979/150‑1 (MQ=255)
                                                          ttACCAGACCCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACATACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACAt                                                                                            >  1:602360/1‑150 (MQ=255)
                                                           tACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTGAGGTTTTCCCCGGCGCCAtt                                                                                           >  2:859521/1‑150 (MQ=255)
                                                            aCCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACAAAATCCGCTGGTCAGGTTTTCCAAGGCGCCATTc                                                                                          >  2:82500/1‑150 (MQ=255)
                                                                  aCCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGtt                                                                                    >  2:1308053/1‑150 (MQ=255)
                                                                   cccGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTta                                                                                   <  2:773571/150‑1 (MQ=255)
                                                                    ccGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTATGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGCCATTCCGCGTtat                                                                                  >  1:626456/1‑150 (MQ=255)
                                                                          gTTGCGTTTGGCATTTAAGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCTACGCAGCAGGTCAGGTGTTCCACGGCGACATTCCCCGTTATATGGCg                                                                            >  2:657531/1‑150 (MQ=255)
                                                                              cGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTg                                                                        >  1:560690/1‑150 (MQ=255)
                                                                                 ttGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGAcc                                                                     <  2:750015/150‑1 (MQ=255)
                                                                                  tGGCAGTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGGCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGAccc                                                                    <  2:162741/150‑1 (MQ=255)
                                                                                   ggCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCg                                                                   >  2:1000738/1‑150 (MQ=255)
                                                                                   ggCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCg                                                                   >  2:12020/1‑150 (MQ=255)
                                                                                    gCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGa                                                                  >  2:1263307/1‑150 (MQ=255)
                                                                                      aTTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACt                                                                >  1:415139/1‑150 (MQ=255)
                                                                                                tCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAAc                                                      >  1:412605/1‑150 (MQ=255)
                                                                                                 cATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAAcg                                                     >  2:1003809/1‑149 (MQ=255)
                                                                                                    cGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTc                                                  >  2:874907/1‑150 (MQ=255)
                                                                                                             ccTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCtt                                         >  1:529393/1‑150 (MQ=255)
                                                                                                             ccTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCtt                                         >  1:472882/1‑150 (MQ=255)
                                                                                                               tGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTcc                                       <  2:1044977/150‑1 (MQ=255)
                                                                                                               tGGTACCGATTGGTCGGCACCACCTCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTGTGCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTcc                                       <  2:836196/150‑1 (MQ=255)
                                                                                                                 gTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATCCACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTg                                     <  1:49017/150‑1 (MQ=255)
                                                                                                                 gTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTg                                     <  2:272179/150‑1 (MQ=255)
                                                                                                                   acgctttGGTCTGCTCCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGCTCAGGTTATCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGtt                                   <  1:1171145/146‑1 (MQ=255)
                                                                                                                     cGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCa                                 >  1:964015/1‑150 (MQ=255)
                                                                                                                      gTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCaa                                <  1:412986/150‑1 (MQ=255)
                                                                                                                          ggTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCACCGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGGTCAAAATg                            >  2:233041/1‑150 (MQ=255)
                                                                                                                            tCTGCACCTCATCTGGAACGTACCTTTCCAGATGTAGATTGGTGAATACTCCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCTTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTa                          <  1:554058/150‑1 (MQ=255)
                                                                                                                              tGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACg                        >  1:413334/1‑150 (MQ=255)
                                                                                                                               gaaccacATCTGGAACGAACCTTACCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACgg                       <  2:214855/148‑1 (MQ=255)
                                                                                                                               gCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCATCGCAGCAGGTCAGGTTTTCCCCGGCGACATTCCGCGGTATAGGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCATCCTGGTCAAAATGTACgg                       >  1:1061863/1‑150 (MQ=255)
                                                                                                                                  ctgcaTCTGGAAAGGTACGTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCt                    >  2:562148/4‑150 (MQ=255)
                                                                                                                                    ccaTCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCTTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGc                  <  2:653045/149‑1 (MQ=255)
                                                                                                                                       tCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACTGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTTATCTGTCTCGTGGCTTCCGGTTCAAAATGTACGGTCTGCCAg               >  2:794298/1‑150 (MQ=255)
                                                                                                                                           gAACGTACCTTTCAGATGCATATTGGTGAATACTCCAACGCAGCAGTTCAGGATTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGcc          <  2:606021/150‑1 (MQ=255)
                                                                                                                                            aaCGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGcc          <  2:1043610/150‑1 (MQ=255)
                                                                                                                                                   cTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTg   >  2:326421/1‑150 (MQ=255)
                                                                                                                                                   cTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGACTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTg   >  1:1298011/1‑150 (MQ=255)
                                                                                                                                                   cTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGATTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAACCTGTGTGGTGGCATCCTGTTCTAAATGTACGGTCTGCCAGCAGCCGCAATTg   >  1:252392/1‑150 (MQ=255)
                                                                                                                                                    tttCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCCGATGCCGCAATTGc  >  2:84344/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTCCTTCATTTGGCCGCCGATTGGTTCTGCAATCCAGACCTTCTCTCAGTGGGCTGCTTACCAGAACCCGGTAGTTGCGTTTGGCATTTACGGTTTCATCGAACGTTGCCTGGTACCGTTTGGTCTGCACCACATCTGGAACGTACCTTTCCAGATGCAGATTGGTGAATACACCAACGCAGCAGGTCAGGTTTTCCACGGCGACATTCCGCGTTATATGGCGGGTGACCCGACTGCGGGTAAACTGTCTGGTGGCTTCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGC  >  NC_000913/1158371‑1158668

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