Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,785,280 G→A G611G (GGC→GGT ppsA ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,785,2800GA97.4% 270.0 / ‑2.2 76G611G (GGC→GGTppsAphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (40/34);  total (42/34)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

TCAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGC  >  NC_000913/1785137‑1785423
                                                                                                                                               |                                                                                                                                               
tcAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCcg                                                                                                                                             <  1:488374/148‑1 (MQ=255)
 cAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGcc                                                                                                                                         >  2:458385/1‑151 (MQ=255)
       ccGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGa                                                                                                                                    >  2:954573/1‑150 (MQ=255)
        cGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAg                                                                                                                                  <  2:598986/151‑1 (MQ=255)
         gCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGc                                                                                                                                 >  2:323742/1‑151 (MQ=255)
            ttCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCg                                                                                                                               >  1:1139649/1‑150 (MQ=255)
            ttCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCg                                                                                                                               >  1:893089/1‑150 (MQ=255)
                 cTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCAt                                                                                                                          >  1:225851/1‑150 (MQ=255)
                   gATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCg                                                                                                                        <  1:3806/150‑1 (MQ=255)
                       tACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTtc                                                                                                                   >  2:623232/1‑151 (MQ=255)
                        aCGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGt                                                                                                                     <  1:288777/148‑1 (MQ=255)
                        aCGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTtct                                                                                                                  <  2:1109704/151‑1 (MQ=255)
                        aCGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTtc                                                                                                                   <  1:733106/150‑1 (MQ=255)
                           gTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCtt                                                                                                               >  1:1042120/1‑151 (MQ=255)
                             aCGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCa                                                                                                             <  1:188479/151‑1 (MQ=255)
                              cgaaCGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCAt                                                                                                            >  1:388668/4‑151 (MQ=255)
                              cGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCAt                                                                                                            <  2:1026991/151‑1 (MQ=255)
                               gCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATc                                                                                                           >  1:323586/1‑151 (MQ=255)
                                  cGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTg                                                                                                         <  1:727823/150‑1 (MQ=255)
                                     aCGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTc                                                                                                     >  1:500960/1‑151 (MQ=255)
                                           tCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAAc                                                                                                >  1:615120/1‑150 (MQ=255)
                                           tCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACg                                                                                               >  2:1408/1‑151 (MQ=255)
                                           tCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACg                                                                                               >  1:1187566/1‑151 (MQ=255)
                                             aTGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGc                                                                                              >  1:37467/1‑150 (MQ=255)
                                             aTGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCt                                                                                             >  1:423091/1‑151 (MQ=255)
                                                 tctcAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTcacc                                                                                         >  2:1038823/1‑151 (MQ=255)
                                                 tctcAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTcacc                                                                                         >  1:683730/1‑151 (MQ=255)
                                                    cAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTcac                                                                                          <  2:855683/147‑1 (MQ=255)
                                                       cGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAc                                                                                   >  1:290871/1‑151 (MQ=255)
                                                        gTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAcc                                                                                  >  1:149296/1‑151 (MQ=255)
                                                         ttGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGAcc                                                                                  >  2:151900/1‑150 (MQ=255)
                                                           ggTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAgg                                                                               >  1:603603/1‑151 (MQ=255)
                                                               aGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTgg                                                                           <  2:145184/151‑1 (MQ=255)
                                                                  cccATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCa                                                                         <  2:581303/150‑1 (MQ=255)
                                                                     aTGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCTGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGATCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGGTGGCATAtt                                                                     >  2:158414/1‑151 (MQ=255)
                                                                     aTGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAtt                                                                     <  1:973233/151‑1 (MQ=255)
                                                                     aTGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAtt                                                                     <  2:487342/151‑1 (MQ=255)
                                                                     aTGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAtt                                                                     <  1:1058816/151‑1 (MQ=255)
                                                                      tGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttc                                                                    >  1:356999/1‑151 (MQ=255)
                                                                      tGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttc                                                                    >  2:697198/1‑151 (MQ=255)
                                                                       gTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCGCCACCGACCAGGTTGGCATAttc                                                                    <  2:256202/150‑1 (MQ=255)
                                                                         cGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgtt                                                                 >  2:937839/1‑151 (MQ=255)
                                                                         cGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgtt                                                                 >  1:1117498/1‑151 (MQ=255)
                                                                         cGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgt                                                                  <  2:318414/150‑1 (MQ=255)
                                                                          gTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttc                                                                <  2:340868/151‑1 (MQ=255)
                                                                            tGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATAttcgttcg                                                               <  1:473305/150‑1 (MQ=255)
                                                                             gcgcACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGAt                                                             <  2:399697/151‑1 (MQ=255)
                                                                              cgcACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGAtt                                                            >  1:587646/1‑151 (MQ=255)
                                                                                  cacGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTaaa                                                        <  2:895629/151‑1 (MQ=255)
                                                                                   acGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTaaa                                                        >  2:51649/1‑150 (MQ=255)
                                                                                    cGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAAt                                                      >  1:722475/1‑151 (MQ=255)
                                                                                     gTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATc                                                     >  2:690817/1‑151 (MQ=255)
                                                                                               ttCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACgg                                           >  1:1246292/1‑151 (MQ=255)
                                                                                                   cacTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACaa                                       <  1:760840/151‑1 (MQ=255)
                                                                                                   cacTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACa                                        <  2:837881/150‑1 (MQ=255)
                                                                                                    acTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAAt                                      >  1:256517/1‑151 (MQ=255)
                                                                                                     cTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATg                                     >  1:377367/1‑151 (MQ=255)
                                                                                                       ccAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATg                                     <  1:1186358/149‑1 (MQ=255)
                                                                                                       ccAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGa                                    <  1:787384/150‑1 (MQ=255)
                                                                                                        cAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGAcg                                  <  1:1121185/151‑1 (MQ=255)
                                                                                                            gcgAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACgcgc                               <  2:505740/150‑1 (MQ=255)
                                                                                                             cgAAACAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACgcgc                               <  2:174770/149‑1 (MQ=255)
                                                                                                                aaCAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCgg                          >  2:672920/1‑151 (MQ=255)
                                                                                                                 aCAGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCgg                          >  1:1016533/1‑150 (MQ=255)
                                                                                                                   aGTCGCGGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGAt                        >  2:360470/1‑150 (MQ=255)
                                                                                                                       gcgGAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAAc                   <  2:256521/151‑1 (MQ=255)
                                                                                                                         ggAAGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAAcgc                 <  2:50376/151‑1 (MQ=255)
                                                                                                                            aGCTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACAACCGACCAGGATGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTAGGATAAAACGCgg               >  2:1101678/1‑150 (MQ=255)
                                                                                                                              cTGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGCGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAAcg                  <  2:801526/145‑1 (MQ=255)
                                                                                                                               tGTCGGAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGACTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCa             >  1:1079812/1‑149 (MQ=255)
                                                                                                                                   ggAAACATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGc       >  2:424375/1‑151 (MQ=255)
                                                                                                                                     aaaCATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGt     >  2:1060653/1‑151 (MQ=255)
                                                                                                                                      aaCATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTc    <  2:588860/151‑1 (MQ=255)
                                                                                                                                      aaCATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTc    <  1:882096/151‑1 (MQ=255)
                                                                                                                                        cATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTcgc  >  2:923946/1‑151 (MQ=255)
                                                                                                                                        cATAGCGACCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTcgc  <  2:269803/151‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TCAACCACCGCTTTCGCCTGATCTACGGTACGCACGAACGGGATCATGATCTCAACGTTGGTCAGTCCCATGTCGTTGCGCACACGTTTCACTGCTTCACACTCCAGCGCGAAACAGTCGCGGAAGCTGTCGGAAACATAGCGGCCCGCGCCACGGAAGCCGAGCATCGGGTTCTCTTCATCTGGCTCGTAACGCTCACCACCGACCAGGTTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGC  >  NC_000913/1785137‑1785423

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