Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 441,071 C→A 5.0% V11F (GTT→TTT) ‡ ispA ← geranyltranstransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913441,0710CA 5.0% 413.7 / 6.3 131V11F (GTT→TTT) ‡ispAgeranyltranstransferase
Reads supporting (aligned to +/- strand):  ref base C (61/63);  new base A (4/3);  total (65/66)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

GCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGC  >  NC_000913/440922‑441215
                                                                                                                                                     |                                                                                                                                                
gCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGTTAGGACTTCTTTAc                                                                                                                                                  >  2:240485/1‑150 (MQ=255)
 ccGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCTCGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACg                                                                                                                                                 >  2:1603144/1‑150 (MQ=255)
 ccGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGAACACGGGAGTGTTCTGAAAGGGCAGTGGGGAGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACg                                                                                                                                                 >  1:920161/1‑150 (MQ=255)
   gAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCa                                                                                                                                               <  2:800935/150‑1 (MQ=255)
      cATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGc                                                                                                                                            >  2:46889/1‑150 (MQ=255)
        tatGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCtt                                                                                                                                          >  2:1537234/1‑150 (MQ=255)
        tatGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTCCTTAACGCAGGCtt                                                                                                                                          >  1:80573/1‑150 (MQ=255)
        aatGACCGGTGGCATAAACCAGGAAAGGGCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCGTTAAGCAGGcca                                                                                                                                          <  1:246571/149‑3 (MQ=255)
          tGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTTACGCAGGCTTCg                                                                                                                                        >  1:1616600/1‑150 (MQ=255)
          tGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCg                                                                                                                                        <  1:41210/150‑1 (MQ=255)
          tGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCg                                                                                                                                        <  2:554235/150‑1 (MQ=255)
          tGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCg                                                                                                                                        <  1:1690484/150‑1 (MQ=255)
          tGACCGGTGGCATAAACCAGGAAAGGGCGCAGGCGCTTACCACCTAATAATGCGCCATACAGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCg                                                                                                                                        <  1:931169/150‑1 (MQ=255)
           gACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGa                                                                                                                                       <  1:693352/150‑1 (MQ=255)
            aCCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATACTGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAg                                                                                                                                      <  1:636087/150‑1 (MQ=255)
             ccGGGGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGt                                                                                                                                     <  2:1061315/150‑1 (MQ=255)
               ggTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtg                                                                                                                                   >  2:655098/1‑150 (MQ=255)
                gTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTGCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgc                                                                                                                                  <  1:586824/150‑1 (MQ=255)
                  ggCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgctg                                                                                                                                <  2:798330/150‑1 (MQ=255)
                   gCATAAACCAGGAAAGGTCGAAGGCGCTTACCACCTAATAATCAGCCATACTGCATGGTTTCGACCACGGGAGTCTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCAGCTTAACGCAGGCTTCGAGTtgctgg                                                                                                                               >  1:1270169/1‑149 (MQ=255)
                    cATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAACCGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCg                                                                                                                              <  1:717636/150‑1 (MQ=255)
                    cATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCg                                                                                                                              <  1:576486/150‑1 (MQ=255)
                    cATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCg                                                                                                                              >  1:1819142/1‑150 (MQ=255)
                        aaCCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCCCGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGaaa                                                                                                                          >  1:1634127/1‑150 (MQ=255)
                          ccAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGt                                                                                                                        >  2:759024/1‑150 (MQ=255)
                          ccAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGt                                                                                                                        <  1:1183485/150‑1 (MQ=255)
                           cAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTc                                                                                                                       <  1:385943/150‑1 (MQ=255)
                             ggAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCa                                                                                                                     >  1:1546956/1‑150 (MQ=255)
                             ggAAAGGTCGCAGGCGCTTACCACCAAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCa                                                                                                                     <  1:316653/150‑1 (MQ=255)
                               aaaGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCCCGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCAtt                                                                                                                   >  1:509137/1‑150 (MQ=255)
                                 aGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTAc                                                                                                                 <  2:741444/150‑1 (MQ=255)
                                   gTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTc                                                                                                               >  2:1629501/1‑150 (MQ=255)
                                     cGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCAt                                                                                                             >  2:1121327/1‑150 (MQ=255)
                                         ggCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGCCCACGGGAGTGTTCTGCAAGGGCAGTGGGGCGATAAACCGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAACGTCCCTTACTCCTTGTc                                                                                                         >  2:278662/1‑150 (MQ=255)
                                         ggCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGGTGCTGCGGAAAGTCCATTACTCATTGTc                                                                                                         >  1:861771/1‑150 (MQ=255)
                                          gcgcTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTcc                                                                                                        <  1:1000046/150‑1 (MQ=255)
                                          gcgcTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTAGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTcc                                                                                                        <  1:93466/150‑1 (MQ=255)
                                           cgcTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCg                                                                                                       <  1:1022407/150‑1 (MQ=255)
                                           cgcTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCg                                                                                                       <  1:1433539/150‑1 (MQ=255)
                                               tACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGgtgt                                                                                                   >  2:853681/1‑150 (MQ=255)
                                               tACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGgtgt                                                                                                   >  2:1266811/1‑150 (MQ=255)
                                                 ccaccTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTaa                                                                                                 >  1:458355/1‑150 (MQ=255)
                                                       aataatGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAgggg                                                                                           >  1:53559/1‑150 (MQ=255)
                                                        ataatGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGt                                                                                          <  2:1660293/150‑1 (MQ=255)
                                                          attGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGTAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCGATTACTCATTGTCCGGTGTAAAAGGGGTTa                                                                                        <  1:1306499/148‑1 (MQ=255)
                                                           atGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTag                                                                                       >  2:1241207/1‑150 (MQ=255)
                                                              cgcCATACTGCATGGTTTCGACCACGGGAGTGATCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGTCCTCCTAAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCAGGTGTAAAAGGGGTTagata                                                                                    >  2:372635/1‑148 (MQ=255)
                                                                ccATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAgg                                                                                  >  2:42628/1‑150 (MQ=255)
                                                                ccATACTGCATGGTTACGACCCCGGGAGTGTGCTGAAAGGGCTGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAgg                                                                                  <  2:137170/150‑1 (MQ=255)
                                                                  aTACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGGTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCg                                                                                <  2:651183/150‑1 (MQ=255)
                                                                  aTACTGCATGGTTTCGACCACGGGAGTGGTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCg                                                                                <  2:626126/150‑1 (MQ=255)
                                                                     cTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTcg                                                                             >  2:1530081/1‑149 (MQ=255)
                                                                       gCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTc                                                                           >  1:988010/1‑150 (MQ=255)
                                                                        cATGGTTTCGACCACGGGAGGGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGTCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGGGTAAAAGGGGTTAGAGAGGCGTCTTCa                                                                          <  1:743752/150‑1 (MQ=255)
                                                                               tCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAg                                                                   <  1:1629501/150‑1 (MQ=255)
                                                                                 gACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGAc                                                                 >  1:1225475/1‑150 (MQ=255)
                                                                                  aCCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGAca                                                                <  1:245623/150‑1 (MQ=255)
                                                                                  aCCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGAca                                                                >  2:1475334/1‑150 (MQ=255)
                                                                                    cACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCAGTACTCATTGTACGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATAGTCATAcagc                                                              >  1:1096423/1‑150 (MQ=255)
                                                                                      cgtgAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTATAAGGGGTTAGAGAGGCGTCTTCATTGTCAGAcagcag                                                            <  2:1633446/147‑1 (MQ=255)
                                                                                      cGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTTCTGCGGAAAGTCCATTACTCATTGACCGGTGTAAAAGGGGTTAGAGAGGCGTCGTCATTGTCAGAcagaag                                                            >  1:1101323/1‑150 (MQ=255)
                                                                                      cGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGAaagcag                                                            >  2:1201069/1‑150 (MQ=255)
                                                                                      cGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCGCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGAcagcag                                                            <  2:1640772/150‑1 (MQ=255)
                                                                                        ggAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGaa                                                          >  2:1024259/1‑150 (MQ=255)
                                                                                            tgtTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTg                                                      <  1:1537234/150‑1 (MQ=255)
                                                                                            tgtTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTg                                                      <  1:1530090/150‑1 (MQ=255)
                                                                                             gtTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTTTCAGACAGCAGAATTTGt                                                     >  2:382820/1‑150 (MQ=255)
                                                                                              ttCTGAAAGGGCAGTGGGGCGATAAAACTGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTa                                                    <  2:1241140/150‑1 (MQ=255)
                                                                                              ttCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTa                                                    <  2:1605987/150‑1 (MQ=255)
                                                                                               tCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTGCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTAc                                                   <  2:402831/150‑1 (MQ=255)
                                                                                                cTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTAcg                                                  <  1:655098/150‑1 (MQ=255)
                                                                                                 tGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCTAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGAAcgc                                                 >  2:1041221/1‑150 (MQ=255)
                                                                                                 tGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTAcgc                                                 >  1:440745/1‑150 (MQ=255)
                                                                                                 tGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTAcgc                                                 >  2:204361/1‑150 (MQ=255)
                                                                                                  gAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTAcgcg                                                >  2:1733130/1‑150 (MQ=255)
                                                                                                   aaaGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACgcgc                                               <  2:521463/150‑1 (MQ=255)
                                                                                                        gCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTc                                          <  2:564357/150‑1 (MQ=255)
                                                                                                            gggggCGATAAAACGGCTCAGCGCCTGGTGGGCCTGCTAAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCt                                      <  1:889971/149‑1 (MQ=255)
                                                                                                             ggggCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCtt                                     >  2:942475/1‑150 (MQ=255)
                                                                                                                   aTAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTa                               >  1:1007364/1‑150 (MQ=255)
                                                                                                                       aaCGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAAttt                           >  2:794719/1‑150 (MQ=255)
                                                                                                                            cTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCt                      >  1:758081/1‑150 (MQ=255)
                                                                                                                            cTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCt                      >  2:725816/1‑150 (MQ=255)
                                                                                                                              cAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGTAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGGCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGc                    <  1:964128/150‑1 (MQ=255)
                                                                                                                              cAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGCCAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGc                    >  2:452697/1‑150 (MQ=255)
                                                                                                                                 cgcCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGcccc                 <  2:177792/150‑1 (MQ=255)
                                                                                                                                 cgcCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGcccc                 >  1:105449/1‑150 (MQ=255)
                                                                                                                                 cgcCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGcccc                 >  1:1367679/1‑150 (MQ=255)
                                                                                                                                   ccTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTg               >  1:685509/1‑150 (MQ=255)
                                                                                                                                   ccTGGTTGGCCTGCTTAACGAAGGCATCGAGTTGCTGCTGAAAGTCCATTAATCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCATg               <  2:864892/150‑1 (MQ=255)
                                                                                                                                     tGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGAc             >  2:1047209/1‑150 (MQ=255)
                                                                                                                                     ctcgaagCCTGCGTTAAGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGAc             >  2:1552254/7‑150 (MQ=255)
                                                                                                                                       cgaagCCTGCGTTAAGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGt           >  2:1211068/5‑150 (MQ=255)
                                                                                                                                       cgaagCCTGCGTTAAGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGt           >  1:2075/5‑150 (MQ=255)
                                                                                                                                        ttGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTg          <  2:909316/150‑1 (MQ=255)
                                                                                                                                        gaagCCTGCGTTAAGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTg          >  2:98400/4‑150 (MQ=255)
                                                                                                                                          ggCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGcc        <  1:1121327/150‑1 (MQ=255)
                                                                                                                                                cTTACCGCAGGCGTCGAGTTGCTGCGGAAAGGCCATTACTCATTGTCCGGCGTAAAAGGGGTTAGAGAGGCGTCATCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGc  <  1:546601/150‑1 (MQ=255)
                                                                                                                                                cTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGc  <  1:853681/150‑1 (MQ=255)
                                                                                                                                                cTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGGCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGc  <  2:963340/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTGC  >  NC_000913/440922‑441215

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