Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,092,780 C→A 6.9% G15V (GGC→GTC)  frvA ← putative PTS enzyme IIA component FrvA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,092,7800CA 6.9% 1013.8 / 38.8 298G15V (GGC→GTC) frvAputative PTS enzyme IIA component FrvA
Reads supporting (aligned to +/- strand):  ref base C (158/119);  new base A (19/2);  total (177/121)
Fisher's exact test for biased strand distribution p-value = 2.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 306 total aligned reads displayed.

CGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAA  >  NC_000913/4092631‑4092928
                                                                                                                                                     |                                                                                                                                                    
   gCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCt                                                                                                                                                   <  1:171585/150‑1 (MQ=255)
    cGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTg                                                                                                                                                  >  1:304394/1‑150 (MQ=255)
    cGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTg                                                                                                                                                  <  2:653161/150‑1 (MQ=255)
    cGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTg                                                                                                                                                  <  1:1240347/150‑1 (MQ=255)
     ggCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGa                                                                                                                                                 <  2:1491311/150‑1 (MQ=255)
      gCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGaa                                                                                                                                                <  2:1960687/150‑1 (MQ=255)
       cACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAat                                                                                                                                               <  1:815147/150‑1 (MQ=255)
              acacCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAgg                                                                                                                                        <  1:487807/150‑1 (MQ=255)
                 ccGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCa                                                                                                                                     <  2:971982/150‑1 (MQ=255)
                   ggAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCa                                                                                                                                          >  2:1106609/1‑143 (MQ=255)
                   ggAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCa                                                                                                                                          <  1:1106609/143‑1 (MQ=255)
                   ggAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAAt                                                                                                                                   >  1:889189/1‑150 (MQ=255)
                   ggAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAAt                                                                                                                                   >  1:1753015/1‑150 (MQ=255)
                       ccAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGACCTGAATATTCAGGGCAAtggcg                                                                                                                               <  1:1161148/150‑5 (MQ=255)
                            tCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGc                                                                                                                          <  1:778195/150‑1 (MQ=255)
                              cAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAg                                                                                                                        <  2:1009864/150‑1 (MQ=255)
                                  agagTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAag                                                                                                                    <  2:466622/150‑1 (MQ=255)
                                            ttttCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCAtt                                                                                                          <  2:1224492/150‑1 (MQ=255)
                                                 gcgcAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCg                                                                                                     <  2:1921820/150‑1 (MQ=255)
                                                 gcgcAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCg                                                                                                     >  2:529971/1‑150 (MQ=255)
                                                   gcagcagcaGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGa                                                                                                   >  2:72351/1‑150 (MQ=255)
                                                    cagcagcagGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGAc                                                                                                  <  2:365649/150‑1 (MQ=255)
                                                         gcagGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGTCTCCAt                                                                                             <  1:690751/150‑1 (MQ=255)
                                                         gcagGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCAt                                                                                             >  1:1639329/1‑150 (MQ=255)
                                                              gTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAat                                                                                        >  2:211082/1‑150 (MQ=255)
                                                              gTTTGCAGGAACTGGTGTCAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAat                                                                                        >  2:800532/1‑150 (MQ=255)
                                                                 tGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAatata                                                                                     <  1:1414255/150‑1 (MQ=255)
                                                                  gCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAatatat                                                                                    >  2:620235/1‑150 (MQ=255)
                                                                     ggAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATAtttt                                                                                 >  1:1887770/1‑150 (MQ=255)
                                                                       aaCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTaa                                                                               <  2:949672/150‑1 (MQ=255)
                                                                        aCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAAt                                                                              <  1:1805209/150‑1 (MQ=255)
                                                                            gtgtGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAAt                                                                          >  1:25749/1‑150 (MQ=255)
                                                                            gtgtGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAAt                                                                          <  2:90558/150‑1 (MQ=255)
                                                                             tgtgAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAAtt                                                                         <  2:864165/150‑1 (MQ=255)
                                                                              gtgAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAgg                                                                                                                                        <  1:1196405/86‑1 (MQ=255)
                                                                              gtgAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAgg                                                                                                                                        >  2:1196405/1‑86 (MQ=255)
                                                                              gtgAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTa                                                                                                                      >  2:1598910/1‑104 (MQ=255)
                                                                              gtgAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTg                                                                        >  1:1951982/1‑150 (MQ=255)
                                                                                 aGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTaa                                                                     >  1:990093/1‑150 (MQ=255)
                                                                                    cGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtgg                                                                  <  2:304394/150‑1 (MQ=255)
                                                                                    cGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtgg                                                                  <  1:654725/150‑1 (MQ=255)
                                                                                     ggTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtggt                                                                 >  1:799906/1‑150 (MQ=255)
                                                                                      gTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtggtg                                                                <  2:424677/150‑1 (MQ=255)
                                                                                         aTAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATaaaa                                                                                         <  2:1652000/122‑1 (MQ=255)
                                                                                          tAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGc                                                            >  1:1134760/1‑150 (MQ=255)
                                                                                           aaaaCCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCg                                                           >  1:503964/1‑150 (MQ=255)
                                                                                           aaaaCCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCg                                                           <  1:529971/150‑1 (MQ=255)
                                                                                               ccGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTa                                                       <  2:1446721/150‑1 (MQ=255)
                                                                                                    ttGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtga                                                  <  2:742709/150‑1 (MQ=255)
                                                                                                       tAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtgatga                                               >  1:1607811/1‑150 (MQ=255)
                                                                                                        aaCGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATgag                                              <  2:508721/150‑1 (MQ=255)
                                                                                                        aaCGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATgag                                              >  2:1347599/1‑150 (MQ=255)
                                                                                                         aCGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATgaga                                             <  2:575759/150‑1 (MQ=255)
                                                                                                            cTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGc                                          <  2:1821976/150‑1 (MQ=255)
                                                                                                            cTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGc                                          <  2:1436656/150‑1 (MQ=255)
                                                                                                            cTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGc                                          >  1:428186/1‑150 (MQ=255)
                                                                                                               ttGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGat                                       <  1:620235/150‑1 (MQ=255)
                                                                                                                tGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGata                                      <  1:368223/150‑1 (MQ=255)
                                                                                                                 gTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGatat                                     >  2:1844462/1‑150 (MQ=255)
                                                                                                                   cGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGt                                   >  1:1583293/1‑150 (MQ=255)
                                                                                                                    gCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTc                                  <  2:1759977/150‑1 (MQ=255)
                                                                                                                    gCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTCGAGCGATTAAGTGATGAGAGGCGATATGTc                                  >  1:1346123/1‑150 (MQ=255)
                                                                                                                      cAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCgg                                >  2:1546797/1‑150 (MQ=255)
                                                                                                                           gCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAAc                           <  2:1068023/150‑1 (MQ=255)
                                                                                                                           gCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAAc                           >  1:517347/1‑150 (MQ=255)
                                                                                                                              tCAGAACGGAATAAGCGCCATTGCGCTGAATATTCAGGTCAATACAGCTTTCAGTAAGAGCTGCCATTAATCTGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAAt                        <  1:1289753/150‑1 (MQ=255)
                                                                                                                              tCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAAt                        >  1:1435815/1‑150 (MQ=255)
                                                                                                                              tCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAAt                        >  1:601474/1‑150 (MQ=255)
                                                                                                                               cAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATg                       >  1:1848599/1‑150 (MQ=255)
                                                                                                                               cAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATg                       >  1:196009/1‑150 (MQ=255)
                                                                                                                                aGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCg                                                                                                     >  1:1210732/1‑71 (MQ=255)
                                                                                                                                aGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGa                      >  1:1975347/1‑150 (MQ=255)
                                                                                                                                aGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGa                      <  2:1583293/150‑1 (MQ=255)
                                                                                                                                aGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGa                      >  2:1046443/1‑150 (MQ=255)
                                                                                                                                aGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGa                      >  1:4622/1‑150 (MQ=255)
                                                                                                                                 gAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGAc                     >  1:570286/1‑150 (MQ=255)
                                                                                                                                 gAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGAc                     >  1:292325/1‑150 (MQ=255)
                                                                                                                                 gAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGAc                     >  1:1517722/1‑150 (MQ=255)
                                                                                                                                 gAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGAc                     >  1:93600/1‑150 (MQ=255)
                                                                                                                                   aCGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACaa                   >  1:1082376/1‑150 (MQ=255)
                                                                                                                                   aCGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACaa                   >  1:167532/1‑150 (MQ=255)
                                                                                                                                      gAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAttt                >  2:1712983/1‑150 (MQ=255)
                                                                                                                                      gAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAttt                >  1:1128316/1‑150 (MQ=255)
                                                                                                                                       aaTAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTg               <  1:1014112/150‑1 (MQ=255)
                                                                                                                                         tAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTgag             <  2:130353/150‑1 (MQ=255)
                                                                                                                                          aaGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTgaga            >  1:1031728/1‑150 (MQ=255)
                                                                                                                                          aaGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTgaga            >  1:963109/1‑150 (MQ=255)
                                                                                                                                           aGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTgagaa           >  1:1319093/1‑149 (MQ=255)
                                                                                                                                            gcgcCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACa          >  1:551154/1‑150 (MQ=255)
                                                                                                                                             ctgtaTTGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACag         >  1:1565304/5‑150 (MQ=255)
                                                                                                                                             cgcCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACag         >  2:702820/1‑150 (MQ=255)
                                                                                                                                              tgtaTTGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACaga        >  1:1496593/4‑150 (MQ=255)
                                                                                                                                               gtaTTGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACagag       >  1:712917/3‑150 (MQ=255)
                                                                                                                                                cATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGt      <  1:322012/150‑1 (MQ=255)
                                                                                                                                                 aTTGACCTGAATATTTAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTg     >  1:1703263/1‑150 (MQ=255)
                                                                                                                                                 aTTGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTg     >  1:838748/1‑150 (MQ=255)
                                                                                                                                                  ttGACCTTAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGa    >  2:1264784/1‑150 (MQ=255)
                                                                                                                                                    gCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGaaa  >  1:1277871/1‑150 (MQ=255)
                                                                                                                                                    gCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGaaa  >  1:1392486/1‑150 (MQ=255)
                                                                                                                                                    gCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGaaa  >  1:1801669/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAA  >  NC_000913/4092631‑4092928

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: