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→AAA) ppsAphosphoenolpyruvate 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

TTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑TTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTCCGG  >  NC_000913/1785347‑1785636
                                                                                                                                                      |                                                                                                                                             
ttGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑t                                                                                                                                              >  1:199719/1‑151 (MQ=255)
 tGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑tt                                                                                                                                             >  2:288067/1‑151 (MQ=255)
 tGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑tt                                                                                                                                             >  2:663297/1‑151 (MQ=255)
  ggCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑tt                                                                                                                                             <  1:108758/150‑1 (MQ=255)
  ggCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑tt                                                                                                                                             <  1:502085/150‑1 (MQ=255)
   gCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑ttt‑t                                                                                                                                          <  2:172545/151‑1 (MQ=255)
    cATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑ttt‑tt                                                                                                                                         <  1:1111778/151‑2 (MQ=255)
    cATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑ttt‑tt                                                                                                                                         >  1:656749/1‑150 (MQ=255)
    cATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑ttt‑tt                                                                                                                                         <  1:358163/151‑2 (MQ=255)
     atatTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑ttt‑ttg                                                                                                                                        >  1:1157398/1‑149 (MQ=255)
        ttcgttcgATTTAACATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑ttt‑ttgca                                                                                                                                      >  1:757365/1‑146 (MQ=255)
        ttcgttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑ttt‑ttgca                                                                                                                                      <  1:1107943/150‑5 (MQ=255)
          cgttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACt                                                                                                                                    <  2:1127985/151‑1 (MQ=255)
          cgttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACt                                                                                                                                    <  2:848070/151‑1 (MQ=255)
          cgttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGACCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACt                                                                                                                                    >  2:921713/1‑151 (MQ=255)
           gttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTg                                                                                                                                   <  1:890281/151‑1 (MQ=255)
           gttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTg                                                                                                                                   >  2:810882/1‑151 (MQ=255)
           gttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTg                                                                                                                                   <  1:495118/151‑1 (MQ=255)
            ttcgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGc                                                                                                                                  >  1:1254135/1‑151 (MQ=255)
              cgATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCgg                                                                                                                                <  2:239801/151‑1 (MQ=255)
               gATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGt                                                                                                                               <  1:593820/151‑1 (MQ=255)
                    aaaaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCt                                                                                                                           <  1:408298/150‑1 (MQ=255)
                    aaaaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTg                                                                                                                          <  1:142414/151‑1 (MQ=255)
                    aaaaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTg                                                                                                                          <  2:101181/151‑1 (MQ=255)
                      aaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCGGACGACCATCGTAAAATTCTCGCGGAGAGTCGAAACCTTTCATCATCTCGCTGATTTCGTTTT‑TGCATCTGCGGTTCCTGAt                                                                                                                        <  1:757752/151‑1 (MQ=255)
                      aaTCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGAt                                                                                                                        >  2:1042138/1‑151 (MQ=255)
                         cAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGt                                                                                                                     >  1:1229802/1‑151 (MQ=255)
                              agaCGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGCTTT‑TGCAACTGCGGTTCCTGATCGTCaaa                                                                                                                 >  2:779126/1‑150 (MQ=255)
                                      aaTGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCa                                                                                                         <  1:111418/150‑1 (MQ=255)
                                      aaTGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAg                                                                                                        <  2:648662/151‑1 (MQ=255)
                                        tGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTg                                                                                                      >  2:209742/1‑151 (MQ=255)
                                            gcgcTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGcgcg                                                                                                  >  2:534174/1‑151 (MQ=255)
                                               cTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCgtgg                                                                                               >  2:690474/1‑151 (MQ=255)
                                               cTTCGGATAAAACGCGGCACCCAGCGAAGCAATCCCTTCAGCGAGACCACCAACGGCGAATCCACCCGGAGAATCCAAACCGATCATCATCTCACGGATTTCGCATT‑TCCAACTGCAGTTCCCGTTCGTCAAACTCCAGCAGTGCACgt                                                                                                 >  1:588469/1‑149 (MQ=255)
                                                 tCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCgtgg                                                                                               <  2:745930/149‑1 (MQ=255)
                                                  cGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCgtgggtg                                                                                            <  2:553700/151‑1 (MQ=255)
                                                   ggATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCgtgggtgg                                                                                           >  1:1194992/1‑151 (MQ=255)
                                                       aaaaCGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGc                                                                                       >  1:306055/1‑151 (MQ=255)
                                                                  cccAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACgg                                                                            <  1:1122977/151‑1 (MQ=255)
                                                                     aGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACggat                                                                          >  2:412112/1‑148 (MQ=255)
                                                                          cgcgATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTtgatga                                                                     <  2:644607/150‑1 (MQ=255)
                                                                           gcgATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATa                                                                   <  1:644325/151‑1 (MQ=255)
                                                                            cgATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATaa                                                                  <  1:84601/151‑1 (MQ=255)
                                                                             gATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAACCCTTTCAACATCTCGCGGATTTCG‑TTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATaaa                                                                 <  1:747425/150‑1 (MQ=255)
                                                                              aTCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAAt                                                                >  1:413592/1‑151 (MQ=255)
                                                                                      agtcagACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGAcg                                                        <  1:1369761/151‑1 (MQ=255)
                                                                                      agtcagACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGAcg                                                        >  2:884565/1‑151 (MQ=255)
                                                                                         cagACGACCAACGTAAAATTCTCGCGGATAATCAAAACCTGTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCATTCATACGGTTGATGATAAATTCCAg                                                           <  1:1219848/145‑1 (MQ=255)
                                                                                          agACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACgcgc                                                     >  2:1280597/1‑150 (MQ=255)
                                                                                            acgacCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAg                                                  <  2:708749/151‑1 (MQ=255)
                                                                                             cgacCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAg                                                  <  1:907292/150‑1 (MQ=255)
                                                                                              gacCAACGTAAAATTCACGCGGAGAATCAAAACATTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGCCAAACTCAAGCATTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAgg                                                 >  1:1280828/1‑150 (MQ=255)
                                                                                                ccAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGccc                                              >  2:275911/1‑151 (MQ=255)
                                                                                                ccAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGccc                                              <  2:1111110/151‑1 (MQ=255)
                                                                                                    cGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGc                                          >  1:1179329/1‑151 (MQ=255)
                                                                                                       aaaaTTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCtt                                       <  2:922812/151‑1 (MQ=255)
                                                                                                       aaaaTTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCtt                                       >  1:603511/1‑151 (MQ=255)
                                                                                                       aaaaTTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCt                                        <  2:1026591/150‑1 (MQ=255)
                                                                                                           ttCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCttcgt                                    >  1:488323/1‑150 (MQ=255)
                                                                                                            tCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCttcgtt                                   <  1:1100142/150‑1 (MQ=255)
                                                                                                             cACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCttcgttcg                                 >  1:603397/1‑151 (MQ=255)
                                                                                                                gcgGAGAATCAAAACCTTTCATCATCTCGCTGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTa                              <  2:790699/151‑1 (MQ=255)
                                                                                                                gcgGAGAATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTa                              <  1:689345/151‑1 (MQ=255)
                                                                                                                    agaATCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTaggca                          <  1:298511/151‑1 (MQ=255)
                                                                                                                      aaTCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTTATCGTCATACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGCCGCGCAAGGCCCACGCCTTCGTTCGGTaggcagg                        <  1:1422660/151‑1 (MQ=255)
                                                                                                                      aaTCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTaggcagg                        <  2:662582/151‑1 (MQ=255)
                                                                                                                       aTCAAAACCTTTCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTaggcagg                        >  1:218139/1‑150 (MQ=255)
                                                                                                                               cTTTCATCATCTCGCCGATTTCGCTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTTCGCGTGGGTGGACGCCAATCATCCGCTTGATGATAAATTCCAGACGCGGAAGGCCCACGCCTTCGTTCGGTAGCCGGGCGACGTCg                >  1:529888/1‑150 (MQ=255)
                                                                                                                                tttCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGa               >  1:562241/1‑150 (MQ=255)
                                                                                                                                tttCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGa               >  2:24632/1‑150 (MQ=255)
                                                                                                                                tttCATCATCTCGCGGATTTCG‑TTTATGGAACTGCGGGTCCTGATCGTCAAACTCAAGAATTGCGTGTGGGTGGACGCCAATAATATGATTGATGATAGAGTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGTTGGTGAAGt                  >  1:471838/1‑147 (MQ=255)
                                                                                                                                 ttCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGaa              >  2:690256/1‑150 (MQ=255)
                                                                                                                                 ttCATCATCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGACGaa              >  2:186442/1‑150 (MQ=255)
                                                                                                                                  tcatcaTCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGTCAATCATACGGTTGATGATAAATTCCAGACGCGCAAAGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGa               >  1:516119/1‑148 (MQ=255)
                                                                                                                                  tcatcaTCTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGaa              >  2:388054/1‑149 (MQ=255)
                                                                                                                                    atcatcTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAACGCCCACGCCTTCGTTAGGTAGGCAGGCGAAGTCGAAAg            >  1:106149/1‑149 (MQ=255)
                                                                                                                                      catcTCGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACg        >  2:869774/1‑151 (MQ=255)
                                                                                                                                          tcGCGGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTcc    >  1:1058306/1‑151 (MQ=255)
                                                                                                                                           cgcgGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTcc    <  2:199728/150‑1 (MQ=255)
                                                                                                                                           cgcgGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTc     <  2:217526/149‑1 (MQ=255)
                                                                                                                                           cgcgGATTTCG‑TTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTc     <  2:907761/148‑1 (MQ=255)
                                                                                                                                            gcgGATTTCGTTTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTCCgg  >  2:1395907/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                             
TTGGCATATTCGTTCGATTTAAAATCAGAGAGACGGACAATGACGCGCTTCGGATAAAACGCGGCACCCAGCGTCGCGATCCCTTCAGTCAGACGACCAACGTAAAATTCACGCGGAGAATCAAAACCTTTCATCATCTCGCGGATTTCG‑TTT‑TGCAACTGCGGTTCCTGATCGTCAAACTCAAGCAGTGCGCGTGGGTGGACGCCAATCATACGGTTGATGATAAATTCCAGACGCGCAAGGCCCACGCCTTCGTTCGGTAGGCAGGCGAAGTCGAAAGCACGGTCCGG  >  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: ‑