Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 440,663 T→G 7.2% S147R (AGC→CGC) ‡ ispA ← geranyl diphosphate/farnesyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913440,6630TG 7.2% 292.5 / 13.4 98S147R (AGC→CGC) ‡ispAgeranyl diphosphate/farnesyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (50/40);  new base G (3/4);  total (54/44)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCA  >  NC_000913/440515‑440809
                                                                                                                                                    |                                                                                                                                                  
cTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTc                                                                                                                                                   >  1:1165939/1‑150 (MQ=255)
  tAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGc                                                                                                                                                 >  1:63357/1‑150 (MQ=255)
      gCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGt                                                                                                                                             >  1:185132/1‑150 (MQ=255)
       cACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGtt                                                                                                                                            >  2:1609742/1‑150 (MQ=255)
       cACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGtt                                                                                                                                            >  2:1542687/1‑150 (MQ=255)
        aCCAAGGCGAACGGCGGCGCGTATCAATGCGCGGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTGCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCATTACCACTGGCGAGCGCCAGtgg                                                                                                                                           <  1:1166056/150‑3 (MQ=255)
         ccAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCa                                                                                                                                          <  2:469110/150‑1 (MQ=255)
           aaGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGa                                                                                                                                        <  1:728510/150‑1 (MQ=255)
           aaGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggtat                                                                                                                                        <  2:404249/150‑6 (MQ=255)
           aaGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggtat                                                                                                                                        <  2:1713970/150‑6 (MQ=255)
            aGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGaa                                                                                                                                       >  2:454430/1‑150 (MQ=255)
               cGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggta                                                                                                                                         <  2:1835826/145‑5 (MQ=255)
               cGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGAGCGCCAGtggta                                                                                                                                         >  1:1835826/1‑141 (MQ=255)
                   cggcggCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCg                                                                                                                                <  1:94212/150‑1 (MQ=255)
                    ggcggcGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGa                                                                                                                               >  2:697780/1‑150 (MQ=255)
                       ggcgCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAAt                                                                                                                            <  2:1824506/150‑1 (MQ=255)
                           cgAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTg                                                                                                                        <  1:1731886/150‑1 (MQ=255)
                              aTCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTcg                                                                                                                     <  1:985379/150‑1 (MQ=255)
                               tCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTcgc                                                                                                                    <  2:63357/150‑1 (MQ=255)
                               tCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTcgc                                                                                                                    <  1:1631202/150‑1 (MQ=255)
                                cAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTcgcg                                                                                                                   >  2:1702619/1‑150 (MQ=255)
                                  aTGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGt                                                                                                                 <  1:201990/150‑1 (MQ=255)
                                    gcgcCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTg                                                                                                                        <  2:1274582/141‑1 (MQ=255)
                                    gcgcCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTg                                                                                                                        >  1:1274582/1‑141 (MQ=255)
                                     cgcCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCg                                                                                                              >  2:466291/1‑150 (MQ=255)
                                     cgcCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCg                                                                                                              >  1:398095/1‑150 (MQ=255)
                                      gcCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGa                                                                                                             >  2:1724391/1‑150 (MQ=255)
                                      gcCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGa                                                                                                             >  2:748887/1‑150 (MQ=255)
                                         ggTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGacac                                                                                                          >  2:534507/1‑150 (MQ=255)
                                             ttATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTcc                                                                                                      >  1:297418/1‑150 (MQ=255)
                                               ttGACGATGAATACGGTCAAGCGCGTCCAGAGGTTCGTGTTTGCCTTCCGCGTCTAAATCTATTGGCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCgg                                                                                                    <  2:1984068/149‑1 (MQ=255)
                                                tGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGc                                                                                                   <  1:746442/150‑1 (MQ=255)
                                                tGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGc                                                                                                   <  2:185132/150‑1 (MQ=255)
                                                 gACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatcatc                                                                                                                                 >  1:657266/1‑119 (MQ=255)
                                                 gACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAatcatc                                                                                                                                 <  2:657266/119‑1 (MQ=255)
                                                 gACGATGAATACGCTCAAGCGCGTCCAGAGGTACCTGTTAGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACAATCCGGCAATACCACAGGCGCACGCCAGATCAGAAAACAACGAAAATCACTCGACCTCCGACACTTcagcag                                                                                                  >  1:1984068/1‑145 (MQ=255)
                                                   cGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATACTGTCGCGGTCCGACACTTCCGGCata                                                                                                >  2:1548212/1‑150 (MQ=255)
                                                      tGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCg                                                                                             >  1:1414992/1‑150 (MQ=255)
                                                      tGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGAACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGCTCCGACACTTCCGGCATATCg                                                                                             >  1:937306/1‑150 (MQ=255)
                                                       gAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCgg                                                                                            <  2:202850/150‑1 (MQ=255)
                                                        aaTACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGc                                                                                           >  1:1949828/1‑150 (MQ=255)
                                                        aaTACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTCCCGGCATATCGGc                                                                                           >  2:485163/1‑150 (MQ=255)
                                                         aTACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCa                                                                                                  <  2:1821351/142‑1 (MQ=255)
                                                         aTACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCa                                                                                                  >  1:1821351/1‑142 (MQ=255)
                                                         aTACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCa                                                                                          >  1:1211071/1‑150 (MQ=255)
                                                         aTACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCa                                                                                          >  2:755398/1‑150 (MQ=255)
                                                             gCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGc                                                                                      <  2:149174/150‑1 (MQ=255)
                                                              cTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCt                                                                                     >  1:1470674/1‑150 (MQ=255)
                                                              cTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCt                                                                                     <  2:306087/150‑1 (MQ=255)
                                                              cTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCt                                                                                     >  1:373616/1‑150 (MQ=255)
                                                                cAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTa                                                                                   >  2:297754/1‑150 (MQ=255)
                                                                     gcgTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATc                                                                              <  1:1853155/150‑1 (MQ=255)
                                                                      cgTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCg                                                                             <  1:1609742/150‑1 (MQ=255)
                                                                       gTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCga                                                                            <  2:1707275/150‑1 (MQ=255)
                                                                           agagGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAAc                                                                        >  2:579027/1‑150 (MQ=255)
                                                                                    tgtTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGtt                                                               >  2:1498590/1‑150 (MQ=255)
                                                                                     gtTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGttt                                                              <  2:768079/150‑1 (MQ=255)
                                                                                        tGCCTTCCGCGTCTAAATCTAATGCCTGACCAGCGCACATTCCGGCAATAACACTGGCGCTCGCCAGTTCAGAATTCATCGAAATTCTGTGGCTGTGCGACAGTTCCGGCATATGGGCTTCGCTTATAATCGAGAACGGCATCGTTTGTa                                                           <  2:589387/150‑1 (MQ=255)
                                                                                         gCCTTCCGCGTCTAAATCTAATGCCTGACAACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTACGACACTTCCGGCATATCGGAATCGCTTAAAATCGAGAACGCCAGCGTTTGTaa                                                          >  1:597542/1‑150 (MQ=255)
                                                                                            ttCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGc                                                       >  1:1146015/1‑150 (MQ=255)
                                                                                             tCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCg                                                      >  2:1264973/1‑150 (MQ=255)
                                                                                             tCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCg                                                      >  2:719150/1‑150 (MQ=255)
                                                                                               cgcgTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTc                                                    >  2:312924/1‑150 (MQ=255)
                                                                                                 cgTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGc                                                  >  2:155370/1‑150 (MQ=255)
                                                                                                   tCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCa                                                >  2:1694654/1‑150 (MQ=255)
                                                                                                      aaaTCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCTTTTGTAAAGCGTCGCCAGCg                                             >  2:737036/1‑150 (MQ=255)
                                                                                                      aaaTCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCg                                             <  2:788257/150‑1 (MQ=255)
                                                                                                        aTCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCgag                                           >  2:183463/1‑150 (MQ=255)
                                                                                                           tAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAAt                                        <  1:978488/150‑1 (MQ=255)
                                                                                                               gCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCatat                                                                                               <  2:164483/91‑1 (MQ=255)
                                                                                                               gCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCatat                                                                                               >  1:164483/1‑91 (MQ=255)
                                                                                                                ccTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGTCACTTCCGGCATATCGGCTTCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGt                                   >  1:1297364/1‑150 (MQ=255)
                                                                                                                ccTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGt                                   >  1:1219832/1‑150 (MQ=255)
                                                                                                                  tGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGttt                                 <  1:454430/150‑1 (MQ=255)
                                                                                                                   gACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGAGTCGCCAGCGAGAATCGCGTTTg                                <  1:168728/150‑1 (MQ=255)
                                                                                                                     ccaccGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCt                              >  1:451028/1‑150 (MQ=255)
                                                                                                                       accGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTc                            <  2:1146015/150‑1 (MQ=255)
                                                                                                                        ccGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCg                           >  2:843414/1‑150 (MQ=255)
                                                                                                                        ccGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCg                           >  2:1561215/1‑150 (MQ=255)
                                                                                                                         cGCACATTCCGGCATTACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCGTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGc                          <  2:1840166/150‑1 (MQ=255)
                                                                                                                          gCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGcc                         >  2:802850/1‑150 (MQ=255)
                                                                                                                             caTTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCaaa                      >  2:585881/1‑150 (MQ=255)
                                                                                                                              aTTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAAc                     <  1:1828692/150‑1 (MQ=255)
                                                                                                                              aTTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGGCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAAc                     <  1:1611657/150‑1 (MQ=255)
                                                                                                                               ttCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACt                    <  1:1542687/150‑1 (MQ=255)
                                                                                                                                 ccGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTc                  <  1:312924/150‑1 (MQ=255)
                                                                                                                                 ccGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTc                  <  1:843414/150‑1 (MQ=255)
                                                                                                                                 ccGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTc                  >  2:1061033/1‑150 (MQ=255)
                                                                                                                                 ccGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTc                  <  2:1165939/150‑1 (MQ=255)
                                                                                                                                   ggCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGATCGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTcac                <  2:398095/150‑1 (MQ=255)
                                                                                                                                        ctgaaCTGGCGAGCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATAGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGc           >  2:1620559/5‑150 (MQ=255)
                                                                                                                                         aCCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCa          >  1:89898/1‑150 (MQ=255)
                                                                                                                                          gaaCTGGCGAGCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAg         >  1:1950560/3‑150 (MQ=255)
                                                                                                                                          ccACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAg         >  1:1546180/1‑150 (MQ=255)
                                                                                                                                           cACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAgg        <  1:755398/150‑1 (MQ=255)
                                                                                                                                               ggCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACt                    >  1:373021/1‑134 (MQ=255)
                                                                                                                                               ggCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACt                    <  2:373021/134‑1 (MQ=255)
                                                                                                                                                 cgcTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCa  <  1:585881/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTTAATGCACCAAGGCGAACGGCGGCGCGAATCAATGCGCCGGTTTTATGACGATGAATACGCTCAAGCGCGTCCAGAGGTACGTGTTTGCCTTCCGCGTCTAAATCTAATGCCTGACCACCGCACATTCCGGCAATACCACTGGCGCTCGCCAGTTCAGAAATCATCGAAATTCTGTCGCGGTCCGACACTTCCGGCATATCGGCATCGCTTAAAATCGAGAACGCCAGCGTTTGTAAAGCGTCGCCAGCGAGAATCGCGTTTGCTTCGCCAAACTTCACATGGCAGGTTGGCA  >  NC_000913/440515‑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: