Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 440,662 C→A 6.2% S147I (AGC→ATC) ‡ ispA ← geranyl diphosphate/farnesyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913440,6620CA 6.2% 275.1 / 6.1 81S147I (AGC→ATC) ‡ispAgeranyl diphosphate/farnesyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (35/41);  new base A (2/3);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

GCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCA  >  NC_000913/440514‑440809
                                                                                                                                                    |                                                                                                                                                   
gcTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGgcgcg                                                                                                                                                    >  2:1375963/1‑149 (MQ=255)
  ttAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCg                                                                                                                                                  <  2:1218374/150‑1 (MQ=255)
   tAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGc                                                                                                                                                 >  2:415564/1‑150 (MQ=255)
    aaTGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGcc                                                                                                                                                <  2:715123/150‑1 (MQ=255)
     aTGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCa                                                                                                                                               >  2:1297642/1‑150 (MQ=255)
      tGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAg                                                                                                                                              >  1:236320/1‑150 (MQ=255)
         aCCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTc                                                                                                                                           >  2:1104317/1‑150 (MQ=255)
          ccAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCa                                                                                                                                          >  1:1062445/1‑150 (MQ=255)
          ccAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggt                                                                                                                                          <  1:915796/150‑4 (MQ=255)
           cAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggta                                                                                                                                         <  1:158420/150‑5 (MQ=255)
            aaGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGa                                                                                                                                        <  2:1100715/150‑1 (MQ=255)
              ggCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGaaa                                                                                                                                      <  2:694863/150‑1 (MQ=255)
                cGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggta                                                                                                                                         <  2:970379/145‑5 (MQ=255)
                cGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggta                                                                                                                                         >  1:970379/1‑141 (MQ=255)
                 gAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAg                                                                                                                                         <  2:176860/144‑1 (MQ=255)
                 gAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAg                                                                                                                                         >  1:176860/1‑144 (MQ=255)
                 gAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatca                                                                                                                                   >  2:250494/1‑150 (MQ=255)
                   aCGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatcatc                                                                                                                                 <  2:1284077/150‑1 (MQ=255)
                    cggcggCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCg                                                                                                                                <  1:40746/150‑1 (MQ=255)
                     ggcggcGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGa                                                                                                                               <  2:946624/150‑1 (MQ=255)
                       cggcgCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGaaa                                                                                                                             >  1:828736/1‑150 (MQ=255)
                        ggcgCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAAt                                                                                                                            <  2:1313824/150‑1 (MQ=255)
                          cgcgAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTc                                                                                                                          >  2:877541/1‑150 (MQ=255)
                            cgAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTg                                                                                                                        <  1:1300510/150‑1 (MQ=255)
                                    tGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTc                                                                                                                <  2:1158523/150‑1 (MQ=255)
                                       gcCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGGCCGa                                                                                                             <  2:1170626/150‑1 (MQ=255)
                                        ccGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGac                                                                                                            <  1:1314846/150‑1 (MQ=255)
                                        ccGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGac                                                                                                            <  1:13327/150‑1 (MQ=255)
                                        ccGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTACGGCAATACCACTGGCGCTCGACAGTTCAGAAATCATCGAAATTCTGTCTCGGTCCGac                                                                                                            <  2:499570/150‑1 (MQ=255)
                                          ggTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGacac                                                                                                          >  2:652053/1‑150 (MQ=255)
                                                aagACGATGAATACGCTCAAGCGCGACCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCgg                                                                                                    >  2:949951/3‑150 (MQ=255)
                                                aTGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCgg                                                                                                    >  2:1257472/1‑150 (MQ=255)
                                                 tGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGc                                                                                                   <  1:1104317/150‑1 (MQ=255)
                                                       tGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCg                                                                                             >  1:1097836/1‑150 (MQ=255)
                                                        gAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCgg                                                                                            <  2:927126/150‑1 (MQ=255)
                                                        gAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCgg                                                                                            <  1:1367241/150‑1 (MQ=255)
                                                         aaTACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGc                                                                                           <  1:250494/150‑1 (MQ=255)
                                                             cGCTCAAGAGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCg                                                                                       >  2:552586/1‑150 (MQ=255)
                                                                  aaGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTaa                                                                                  <  2:925449/150‑1 (MQ=255)
                                                                     cgcgTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAAt                                                                               >  1:283045/1‑150 (MQ=255)
                                                                         tCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCgag                                                                           >  1:466370/1‑150 (MQ=255)
                                                                         tCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCgag                                                                           >  2:1091612/1‑150 (MQ=255)
                                                                          ccAGAGGTACGTGTTTGCCTTCCGCCTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCTTCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCgaga                                                                          <  1:422192/150‑1 (MQ=255)
                                                                           cAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGaa                                                                         >  2:1043256/1‑150 (MQ=255)
                                                                             gagGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACg                                                                       <  1:218039/150‑1 (MQ=255)
                                                                              agGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGc                                                                      <  1:1389370/150‑1 (MQ=255)
                                                                                     tgtTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGtt                                                               >  2:1133781/1‑150 (MQ=255)
                                                                                     tgtTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGtt                                                               <  1:1050204/150‑1 (MQ=255)
                                                                                     tgtTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGtt                                                               >  2:94318/1‑150 (MQ=255)
                                                                                     tgtTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGtt                                                               >  1:535555/1‑150 (MQ=255)
                                                                                       tttGCCTTCCGCGTCTAAATCTAATGCCTGACCTCCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTg                                                             <  1:716962/150‑1 (MQ=255)
                                                                                       tttGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTg                                                             <  2:1062445/150‑1 (MQ=255)
                                                                                        ttGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGt                                                            >  1:685911/1‑150 (MQ=255)
                                                                                             ttCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGc                                                       <  1:1257472/150‑1 (MQ=255)
                                                                                              tCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCg                                                      <  2:1097836/150‑1 (MQ=255)
                                                                                              tCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCg                                                      >  1:400398/1‑150 (MQ=255)
                                                                                                 gcgTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCg                                                   >  1:1127362/1‑150 (MQ=255)
                                                                                                     cTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAg                                               <  1:52300/150‑1 (MQ=255)
                                                                                                         aTCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGTGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCgag                                           >  1:169200/1‑150 (MQ=255)
                                                                                                         aTCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCgag                                           >  1:227343/1‑150 (MQ=255)
                                                                                                          tCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCgaga                                          <  1:1385952/150‑1 (MQ=255)
                                                                                                             aaTGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGTCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATc                                       <  2:304127/150‑1 (MQ=255)
                                                                                                              aTGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATcg                                      <  1:360735/150‑1 (MQ=255)
                                                                                                                gCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATcgcg                                    >  1:894917/1‑150 (MQ=255)
                                                                                                                  cTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGtt                                  <  2:588510/150‑1 (MQ=255)
                                                                                                                   tGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGttt                                 <  1:652053/150‑1 (MQ=255)
                                                                                                                    gACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTg                                <  1:415564/150‑1 (MQ=255)
                                                                                                                    gACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTg                                <  1:202661/150‑1 (MQ=255)
                                                                                                                    gACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTg                                <  1:877541/150‑1 (MQ=255)
                                                                                                                      ccaccGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCt                              >  1:929733/1‑150 (MQ=255)
                                                                                                                        accGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTc                            >  1:1325348/1‑150 (MQ=255)
                                                                                                                        accGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTc                            >  1:1119124/1‑150 (MQ=255)
                                                                                                                        accGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTc                            >  2:1072099/1‑150 (MQ=255)
                                                                                                                          cGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGc                          >  2:115704/1‑150 (MQ=255)
                                                                                                                           gCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGcc                         <  1:1091612/150‑1 (MQ=255)
                                                                                                                               aTTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAAc                     <  1:94318/150‑1 (MQ=255)
                                                                                                                                 tCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACtt                   >  1:1360506/1‑150 (MQ=255)
                                                                                                                                   cGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTca                 >  1:1186822/1‑150 (MQ=255)
                                                                                                                                     gCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTcaca               <  2:283045/150‑1 (MQ=255)
                                                                                                                                     gCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTcaca               <  1:949951/150‑1 (MQ=255)
                                                                                                                                        aTACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATgg            >  2:71048/1‑150 (MQ=255)
                                                                                                                                         ctgaaCTGGCGAGCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGc           >  1:109721/5‑150 (MQ=255)
                                                                                                                                                 gcggTCGCCAGTTCAGAAATCATCGAAATTCTGGGGCGGTCCGACAGTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGc   <  2:373213/150‑1 (MQ=255)
                                                                                                                                                 gcgcTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGc   <  1:1372601/150‑1 (MQ=255)
                                                                                                                                                  cgcTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCa  <  2:236320/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCA  >  NC_000913/440514‑440809

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: