Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,965,907 C→T 6.4% L95L (CTG→CTA flhB ← flagellar biosynthesis protein FlhB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,965,9070CT 6.4% 242.8 / 4.5 78L95L (CTG→CTAflhBflagellar biosynthesis protein FlhB
Reads supporting (aligned to +/- strand):  ref base C (34/39);  new base T (2/3);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

ACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACC  >  NC_000913/1965759‑1966055
                                                                                                                                                    |                                                                                                                                                    
aCTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCa                                                                                                                                                     <  1:396632/150‑1 (MQ=255)
   ccGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCg                                                                                                                                                  >  1:283406/1‑150 (MQ=255)
   ccGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCg                                                                                                                                                  <  1:1234832/150‑1 (MQ=255)
      cgcCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGgca                                                                                                                                               >  2:444383/1‑150 (MQ=255)
        ccAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTGATTAGCGgcgtg                                                                                                                                             <  2:734205/150‑4 (MQ=255)
         cAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTGATTAGCGgcgtgg                                                                                                                                            <  1:122083/150‑5 (MQ=255)
           gTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCagcagc                                                                                                                                          >  2:250923/1‑150 (MQ=255)
            tCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTGATTAGCGgcgtggtgc                                                                                                                                         <  1:316456/150‑8 (MQ=255)
              tGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGcc                                                                                                                                       <  2:283406/150‑1 (MQ=255)
              tGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCAAATCAGCGGCAGCAGCGcc                                                                                                                                       <  2:118638/150‑1 (MQ=255)
                 gCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGc                                                                                                                                    >  1:1199947/1‑150 (MQ=255)
                 gCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGc                                                                                                                                    >  1:363766/1‑150 (MQ=255)
                 gCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGc                                                                                                                                    <  2:575532/150‑1 (MQ=255)
                          aTCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTtct                                                                                                                           >  1:181211/1‑150 (MQ=255)
                            ccGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTctct                                                                                                                         <  2:1190407/150‑1 (MQ=255)
                             cGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTg                                                                                                                        <  1:1135214/150‑1 (MQ=255)
                             cGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTg                                                                                                                        >  2:876069/1‑150 (MQ=255)
                                ttAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATc                                                                                                                     <  1:133141/150‑1 (MQ=255)
                                  taTGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATcag                                                                                                                   <  1:830445/149‑1 (MQ=255)
                                    tGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATcagca                                                                                                                 <  1:10340/150‑1 (MQ=255)
                                       ccGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGaa                                                                                                              <  2:242567/150‑1 (MQ=255)
                                             gCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTg                                                                                                                        >  1:244135/1‑134 (MQ=255)
                                             gCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTg                                                                                                                        <  2:244135/134‑1 (MQ=255)
                                                     aGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAgg                                                                                                >  1:60554/1‑150 (MQ=255)
                                                     aGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAgg                                                                                                >  1:864724/1‑150 (MQ=255)
                                                     aGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCACAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAgg                                                                                                >  2:1008599/1‑150 (MQ=255)
                                                      gTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGa                                                                                               >  1:1219750/1‑150 (MQ=255)
                                                        ttGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATc                                                                                             <  2:181211/150‑1 (MQ=255)
                                                        ttGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATc                                                                                             >  2:1006493/1‑150 (MQ=255)
                                                         tGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCa                                                                                            >  2:1073490/1‑150 (MQ=255)
                                                           gAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGa                                                                                          >  1:163904/1‑150 (MQ=255)
                                                             aaaCTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGAtt                                                                                        <  1:567246/150‑1 (MQ=255)
                                                               aCTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCg                                                                                      <  1:387211/150‑1 (MQ=255)
                                                                cTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCgg                                                                                     <  1:960528/150‑1 (MQ=255)
                                                                cTTCGGCTGCAAGGATTTGCCGCTAAATACCACCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCgg                                                                                     <  2:523893/150‑1 (MQ=255)
                                                                 ttCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCggg                                                                                    <  1:1226441/150‑1 (MQ=255)
                                                                      cTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCAtt                                                                               <  1:1011849/150‑1 (MQ=255)
                                                                      cTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGAGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCAtt                                                                               <  2:869474/150‑1 (MQ=255)
                                                                        gCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCa                                                                                            <  1:842834/135‑1 (MQ=255)
                                                                        gCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCa                                                                                            >  2:842834/1‑135 (MQ=255)
                                                                                  gCCGATAAATACCAGCCCTCCAATCATGACCGGATAAATGAGCGCCACCAGAACAACGCCGCTAATCAGCGGCAGCAGAGCCATAATGGCTTCTCTGATCAGCAGAATATTCTGCCCGTGGATCAGATTCGGGTCATTGATAATACtttg                                                                   <  2:326358/150‑1 (MQ=255)
                                                                                   ccGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGa                                                                  >  2:1138729/1‑150 (MQ=255)
                                                                                    cGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGAt                                                                 >  1:237133/1‑150 (MQ=255)
                                                                                    cGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGAt                                                                 <  1:979930/150‑1 (MQ=255)
                                                                                      cTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCa                                                               <  2:1199947/150‑1 (MQ=255)
                                                                                      cTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCa                                                               <  1:386920/150‑1 (MQ=255)
                                                                                         aaTACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCaaaa                                                            <  2:1017591/150‑1 (MQ=255)
                                                                                           tACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATg                                                          <  1:1096328/150‑1 (MQ=255)
                                                                                            aCCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGc                                                         <  2:60554/150‑1 (MQ=255)
                                                                                            aCCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGc                                                         >  1:625626/1‑150 (MQ=255)
                                                                                             ccAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCTTTGAAAATACTCTGAACAAAATGCa                                                        >  1:624204/1‑150 (MQ=255)
                                                                                             ccAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCa                                                        >  2:657310/1‑150 (MQ=255)
                                                                                               aGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGc                                                      <  1:1358210/150‑1 (MQ=255)
                                                                                                   cTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTtctc                                                                                                                          <  2:985065/78‑1 (MQ=255)
                                                                                                   cTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTtctc                                                                                                                          >  1:985065/1‑78 (MQ=255)
                                                                                                   cTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCTTTGATAATACTGTGATCAAAATGCAGCCCAg                                                  >  1:795606/1‑150 (MQ=255)
                                                                                                       cAGCATCACAGGAGAAATAAACACAACAAGCACCACGACCATAATCAGCGGCACCAGCGCCAGAATAGCATCTCTGATAAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGga                                              >  2:1413281/1‑150 (MQ=255)
                                                                                                           aTGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCATAATTATTTGCCTGAGTATGAGATTAGGGTGATTGATAATAGTGTTATCAATATGCATCCCAGAGgagatg                                          <  2:589324/150‑3 (MQ=255)
                                                                                                               ccGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGc                                      >  1:274246/1‑150 (MQ=255)
                                                                                                                   agaAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGa                                  <  1:1181950/150‑1 (MQ=255)
                                                                                                                    gaAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGAc                                 >  1:344268/1‑150 (MQ=255)
                                                                                                                      aaTGAGCGCCACCAGCACCACCCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCAGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACaa                               >  2:494016/1‑150 (MQ=255)
                                                                                                                          aGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGa                           >  1:1082682/1‑150 (MQ=255)
                                                                                                                              ccaccaGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACggg                       <  1:40096/150‑1 (MQ=255)
                                                                                                                               caccaGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGTCGATTGTCGGGc                      >  1:608777/1‑150 (MQ=255)
                                                                                                                                accaGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGTCGATTGTCGGGca                     <  2:1219921/150‑2 (MQ=255)
                                                                                                                                accaGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGTCGATTGTCGGGca                     <  2:608777/150‑2 (MQ=255)
                                                                                                                                   aGCACCACTCCTCTAATCAGCGGCAGCAGAGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGACAGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGAAGACCAGCGGAGAGCATTCCATACAATCGAAGGGCCAGc                  >  2:642721/1‑150 (MQ=255)
                                                                                                                                   aGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGc                  >  1:333234/1‑150 (MQ=255)
                                                                                                                                    gCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGTCGATTGTCGGGcatgct                 <  2:625626/150‑6 (MQ=255)
                                                                                                                                      ctgctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGac               >  2:1277152/6‑150 (MQ=255)
                                                                                                                                        gctgccgcTGATTAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGacac             >  2:1143489/4‑150 (MQ=255)
                                                                                                                                        cacGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGccac             >  1:1309790/1‑150 (MQ=255)
                                                                                                                                          cgccgcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGacacac           >  1:963308/1‑150 (MQ=255)
                                                                                                                                             cgcTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGc        <  1:1006493/150‑1 (MQ=255)
                                                                                                                                               cTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCa      <  1:250923/150‑1 (MQ=255)
                                                                                                                                                  aTCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAAc   <  1:76043/150‑1 (MQ=255)
                                                                                                                                                   tCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAAcc  >  1:964801/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGGGTTGAGTTTGGAAAACTTCGGCTGCAAGGATTTGCCGCTAAATACCAGCCCTCCCAGCATGACCGGAGAAATGAGCGCCACCAGCACCACGCCGCTAATCAGCGGCAGCAGCGCCAGCATGGCTTCTCTGATCAGCAGAATAATCTGCCCGAGGATCAGATTCGGGTCATTGATAATACTGTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACC  >  NC_000913/1965759‑1966055

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