Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,785,500 +T coding (1613/2379 nt) ppsA ← phosphoenolpyruvate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,785,4961.T95.5% 227.4 / ‑2.6 70N539K (AAC→AAAppsAphosphoenolpyruvate synthase
Reads supporting (aligned to +/- strand):  ref base . (1/2);  major base T (31/33);  minor base C (3/0);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTCCGG  >  NC_000913/1785347‑1785636
                                                                                                                                                      |                                                                                                                                             
ttGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGt                                                                                                                                              >  1:199719/1‑151 (MQ=255)
 tGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGtt                                                                                                                                             >  2:288067/1‑151 (MQ=255)
 tGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGtt                                                                                                                                             >  2:663297/1‑151 (MQ=255)
  ggCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGtt                                                                                                                                             <  1:108758/150‑1 (MQ=255)
  ggCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGtt                                                                                                                                             <  1:502085/150‑1 (MQ=255)
   gCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGtttt                                                                                                                                          <  2:172545/151‑1 (MQ=255)
    cATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGttttt                                                                                                                                         <  1:1111778/151‑2 (MQ=255)
    cATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGttttt                                                                                                                                         >  1:656749/1‑150 (MQ=255)
    cATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGttttt                                                                                                                                         <  1:358163/151‑2 (MQ=255)
     atatTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGtttttg                                                                                                                                        >  1:1157398/1‑149 (MQ=255)
        ttcgttcgATTTAACATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGtttttgca                                                                                                                                      >  1:757365/1‑146 (MQ=255)
        ttcgttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGtttttgca                                                                                                                                      <  1:1107943/150‑5 (MQ=255)
          cgttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACt                                                                                                                                    <  2:1127985/151‑1 (MQ=255)
          cgttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACt                                                                                                                                    <  2:848070/151‑1 (MQ=255)
          cgttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGACCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACt                                                                                                                                    >  2:921713/1‑151 (MQ=255)
           gttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTg                                                                                                                                   <  1:890281/151‑1 (MQ=255)
           gttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTg                                                                                                                                   >  2:810882/1‑151 (MQ=255)
           gttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTg                                                                                                                                   <  1:495118/151‑1 (MQ=255)
            ttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGc                                                                                                                                  >  1:1254135/1‑151 (MQ=255)
              cgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCgg                                                                                                                                <  2:239801/151‑1 (MQ=255)
               gATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGt                                                                                                                               <  1:593820/151‑1 (MQ=255)
                    aaaaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCt                                                                                                                           <  1:408298/150‑1 (MQ=255)
                    aaaaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTg                                                                                                                          <  1:142414/151‑1 (MQ=255)
                    aaaaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTg                                                                                                                          <  2:101181/151‑1 (MQ=255)
                      aaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCGGACGACCATCGTAAAATTCTCGCGGAGAGTCGAAACCTTTCATCATCTCGCTGATTTCGTTTTTGCATCTGCGGTTCCTGAt                                                                                                                        <  1:757752/151‑1 (MQ=255)
                      aaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGAt                                                                                                                        >  2:1042138/1‑151 (MQ=255)
                         cAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGt                                                                                                                     >  1:1229802/1‑151 (MQ=255)
                              agaCGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGCTTTTGCAACTGCGGTTCCTGATCGTCaaa                                                                                                                 >  2:779126/1‑150 (MQ=255)
                                      aaTGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCa                                                                                                         <  1:111418/150‑1 (MQ=255)
                                      aaTGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAg                                                                                                        <  2:648662/151‑1 (MQ=255)
                                        tGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTg                                                                                                      >  2:209742/1‑151 (MQ=255)
                                            gcgcTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGcgcg                                                                                                  >  2:534174/1‑151 (MQ=255)
                                               cTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCgtgg                                                                                               >  2:690474/1‑151 (MQ=255)
                                               cTTCGGATAAAACGCGGCACCCAGCGAAGCAATCCCTTCAGCGAGACCACCAACGGCGAATCCACCCGGAGAATCCAAACCGATCATCATCTCACGGATTTCGCATTTCCAACTGCAGTTCCCGTTCGTCAAACTCCAGCAGTGCACgt                                                                                                 >  1:588469/1‑149 (MQ=255)
                                                 tCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCgtgg                                                                                               <  2:745930/149‑1 (MQ=255)
                                                  cGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCgtgggtg                                                                                            <  2:553700/151‑1 (MQ=255)
                                                   ggATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCgtgggtgg                                                                                           >  1:1194992/1‑151 (MQ=255)
                                                       aaaaCGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGc                                                                                       >  1:306055/1‑151 (MQ=255)
                                                                  cccAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACgg                                                                            <  1:1122977/151‑1 (MQ=255)
                                                                     aGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACggat                                                                          >  2:412112/1‑148 (MQ=255)
                                                                          cgcgATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTtgatga                                                                     <  2:644607/150‑1 (MQ=255)
                                                                           gcgATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATa                                                                   <  1:644325/151‑1 (MQ=255)
                                                                            cgATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATaa                                                                  <  1:84601/151‑1 (MQ=255)
                                                                             gATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAACCCTTTCAACATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATaaa                                                                 <  1:747425/150‑1 (MQ=255)
                                                                              aTCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAAt                                                                >  1:413592/1‑151 (MQ=255)
                                                                                      agtcagACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGAcg                                                        <  1:1369761/151‑1 (MQ=255)
                                                                                      agtcagACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGAcg                                                        >  2:884565/1‑151 (MQ=255)
                                                                                         cagACGACCAACGTAAAATTCTCGCGGATAATCAAAACCTGTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCATTCATACGGTTGATGATAAATTCCAg                                                           <  1:1219848/145‑1 (MQ=255)
                                                                                          agACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACgcgc                                                     >  2:1280597/1‑150 (MQ=255)
                                                                                            acgacCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAg                                                  <  2:708749/151‑1 (MQ=255)
                                                                                             cgacCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAg                                                  <  1:907292/150‑1 (MQ=255)
                                                                                              gacCAACGTAAAATTCACGCGGAGAATCAAAACATTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGCCAAACTCAAGCATTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAgg                                                 >  1:1280828/1‑150 (MQ=255)
                                                                                                ccAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGccc                                              >  2:275911/1‑151 (MQ=255)
                                                                                                ccAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGccc                                              <  2:1111110/151‑1 (MQ=255)
                                                                                                    cGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGc                                          >  1:1179329/1‑151 (MQ=255)
                                                                                                       aaaaTTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCtt                                       <  2:922812/151‑1 (MQ=255)
                                                                                                       aaaaTTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCtt                                       >  1:603511/1‑151 (MQ=255)
                                                                                                       aaaaTTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCt                                        <  2:1026591/150‑1 (MQ=255)
                                                                                                           ttCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCttcgt                                    >  1:488323/1‑150 (MQ=255)
                                                                                                            tCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCttcgtt                                   <  1:1100142/150‑1 (MQ=255)
                                                                                                             cACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCttcgttcg                                 >  1:603397/1‑151 (MQ=255)
                                                                                                                gcgGAGAATCAAAACCTTTCATCATCTCGCTGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTa                              <  2:790699/151‑1 (MQ=255)
                                                                                                                gcgGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTa                              <  1:689345/151‑1 (MQ=255)
                                                                                                                    agaATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTaggca                          <  1:298511/151‑1 (MQ=255)
                                                                                                                      aaTCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTTATCGTCATACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGCCGCGCAAGGCCCACGCCTTCGTTCGGTaggcagg                        <  1:1422660/151‑1 (MQ=255)
                                                                                                                      aaTCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTaggcagg                        <  2:662582/151‑1 (MQ=255)
                                                                                                                       aTCAAAACCTTTCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTaggcagg                        >  1:218139/1‑150 (MQ=255)
                                                                                                                               cTTTCATCATCTCGCCGATTTCGCTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTTCGCGTGGGTGGACGCCAATCATCCGCTTGATGATAAATTCCAGACGCGGAAGGCCCACGCCTTCGTTCGGTAGCCGGGCGACGTCg                >  1:529888/1‑150 (MQ=255)
                                                                                                                                tttCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGa               >  1:562241/1‑150 (MQ=255)
                                                                                                                                tttCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGa               >  2:24632/1‑150 (MQ=255)
                                                                                                                                tttCATCATCTCGCGGATTTCGTTTATGGAACTGCGGGTCCTGATCGTCAAACTCAAGAATTGCGTGTGGGTGGACGCCAATAATATGATTGATGATAGAGTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGTTGGTGAAGt                  >  1:471838/1‑147 (MQ=255)
                                                                                                                                 ttCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGaa              >  2:690256/1‑150 (MQ=255)
                                                                                                                                 ttCATCATCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGACGaa              >  2:186442/1‑150 (MQ=255)
                                                                                                                                  tcatcaTCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGTCAATCATACGGTTGATGATAAATTCCAGACGCGCAAAGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGa               >  1:516119/1‑148 (MQ=255)
                                                                                                                                  tcatcaTCTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGaa              >  2:388054/1‑149 (MQ=255)
                                                                                                                                    atcatcTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAACGCCCACGCCTTCGTTAGGTAGGCAGGCGAAGTCGAAAg            >  1:106149/1‑149 (MQ=255)
                                                                                                                                      catcTCGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACg        >  2:869774/1‑151 (MQ=255)
                                                                                                                                          tcGCGGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTcc    >  1:1058306/1‑151 (MQ=255)
                                                                                                                                           cgcgGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTcc    <  2:199728/150‑1 (MQ=255)
                                                                                                                                           cgcgGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTc     <  2:217526/149‑1 (MQ=255)
                                                                                                                                           cgcgGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTc     <  2:907761/148‑1 (MQ=255)
                                                                                                                                            gcgGATTTCGTTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTCCgg  >  2:1395907/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
TTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTTGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTCCGG  >  NC_000913/1785347‑1785636

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