Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,288,033 T→G 5.6% intergenic (+23/‑57) yraH → / → yraI putative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,288,0330TG 5.6% 366.7 / 7.5 124intergenic (+23/‑57)yraH/yraIputative fimbrial‑like adhesin protein/putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (45/72);  new base G (3/4);  total (48/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

Only 100 of 124 total aligned reads displayed.

CGGCGGAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAAC  >  NC_000913/3287886‑3288181
                                                                                                                                                   |                                                                                                                                                    
cggcggAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGCCCTGt                                                                                                                                                    >  1:1762306/1‑150 (MQ=255)
cggcggAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGt                                                                                                                                                    >  2:515987/1‑150 (MQ=255)
       aTCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATaa                                                                                                                                             >  1:399078/1‑150 (MQ=255)
          cAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTTTCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGATGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGc                                                                                                                                          <  2:208396/150‑1 (MQ=255)
          cAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGCCCTGTCAGATAACGc                                                                                                                                          >  2:239376/1‑150 (MQ=255)
           aaTGGTATTGTCAATGTTGACCTCATCGCGTCCTACGTATCAACCGTAGTTCCATCAGAAATCTCCCCAGGGGAAGCAGACGCGATCGTTAACGTAACGCTGGATTACCGATACTACGTTAGGGCGTTATCTGACAGGGCAGATAACGcc                                                                                                                                         <  2:1004348/150‑1 (MQ=255)
            aTGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGccc                                                                                                                                        <  2:1656469/150‑1 (MQ=255)
            aTGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGccc                                                                                                                                        <  1:170715/150‑1 (MQ=255)
            aTGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACAGGTCAGATAACGccc                                                                                                                                        <  2:1537/150‑1 (MQ=255)
             tGGTATTGTCAATGTTGACCTCATCGCTTACTACGTTACAACCGAAGTTCCATCACAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCt                                                                                                                                       <  2:1682272/150‑1 (MQ=255)
              ggTATTGTCAATGTTGACCTCATCGCTTACTACGTATCATCCGTAGTTCCATCAGAAATCACCCCATGGGAAGCAGACGCTATCGTTAACGTACCGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGGCAGATAAGGCCCtt                                                                                                                                      >  2:203177/1‑150 (MQ=255)
               gtcttGGCAATGCTCACCTCATCGCTTACTACGTATCAACCGACTTACCATCACAAATAACCCCCGGCGCCCCAGACGCTAGCGTTAACGTCACGCTCGAATACCGATAATACGTTAGGGCGTAATCAGACCTCTCAGATAACGCCCtta                                                                                                                                     <  2:263663/147‑2 (MQ=255)
                  ttGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttcc                                                                                                                                  <  1:449408/150‑1 (MQ=255)
                     tCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCGGATAACTCCCTTttccttc                                                                                                                               >  1:1304145/1‑150 (MQ=255)
                      cAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTttccttcc                                                                                                                              >  1:193573/1‑150 (MQ=255)
                        aTGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctc                                                                                                                            <  1:1751481/150‑1 (MQ=255)
                         tGTTGACCTCATCGCTTACTCCGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCctct                                                                                                                           <  1:679086/150‑1 (MQ=255)
                          gTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCtt                                                                                                                          <  1:231207/150‑1 (MQ=255)
                            tGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTtc                                                                                                                        >  2:96497/1‑150 (MQ=255)
                             gcccTCATCGCATACTACGTATCAACCGTAGTTCCATCAGAAAGCACCCCAGCGGAAGCAGACGCTATCATTAACGTAACGCTGGATTACCGTTAAAACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTtct                                                                                                                       <  2:346183/148‑1 (MQ=255)
                                 tCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGtt                                                                                                                   <  2:701185/150‑1 (MQ=255)
                                   aTCGCTTACTTCGTATCAACCGGAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATTCGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCGTTTCCGTCCGGTTTCTCGTTGt                                                                                                                 >  2:525894/1‑150 (MQ=255)
                                    tCGCTTACTACGTATCATCCGTAGTTCCATCTGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTa                                                                                                                <  1:398612/150‑1 (MQ=255)
                                    tCGCTTACTACGTATCAACCGTAGTTCCATCAGAACTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTGTCTCGTTGtt                                                                                                                >  1:748883/1‑149 (MQ=255)
                                    tCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTa                                                                                                                >  1:1485774/1‑150 (MQ=255)
                                       cTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTTTCa                                                                                                             >  2:807234/1‑150 (MQ=255)
                                         tactacGTATCTACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTAAATACGTTTGGGCGTTATCTGACCTGTCATATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAgg                                                                                                           <  1:1539754/150‑1 (MQ=255)
                                          actacGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTCTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGt                                                                                                          <  2:1763916/150‑1 (MQ=255)
                                           ctacGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCGTTCCTCTTTCTCGTTGTATCAGGtt                                                                                                         >  1:312062/1‑150 (MQ=255)
                                           atacGTATCAACCGTAGTTCCATCAGAAATCACCCCTGGGGAAGCAGACGCTATCGTAAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGtt                                                                                                         <  2:873388/149‑1 (MQ=255)
                                            tacGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTg                                                                                                        <  2:729182/150‑1 (MQ=255)
                                            tacGTATCAACCGTAGTTCCATCACAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCACATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTg                                                                                                        <  2:624826/150‑1 (MQ=255)
                                              cGTATCAACCGTAGTTCCATCAGAAATCACCCCTGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGaa                                                                                                      <  2:1774099/150‑1 (MQ=255)
                                               gTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGaaa                                                                                                     <  1:1611061/150‑1 (MQ=255)
                                                 aTCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATATGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGaaaaa                                                                                                   >  2:1550607/1‑150 (MQ=255)
                                                   cAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATa                                                                                                 <  1:670183/150‑1 (MQ=255)
                                                     aCCTTAGTTCCATCAGAAATCACCCCTGGGGAAGCTGACGCTATCGTTAACGTAACGGTGTATTACCGTAAATACGTTAGGGCGTTATCTGACCTGTCATATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAAc                                                                                               <  2:1415844/150‑1 (MQ=255)
                                                     aCCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAAc                                                                                               <  2:237534/150‑1 (MQ=255)
                                                     aCCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTGCCTCATTCTCGATGTATCAGCTTGAAAAATAAc                                                                                               >  1:293073/1‑150 (MQ=255)
                                                       cGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCCTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTa                                                                                             <  1:170035/150‑1 (MQ=255)
                                                       cGTAGTTCCATCAGAAACCACCCCAGGCGAAGCAGACGCTATCCTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTa                                                                                             <  2:146352/150‑1 (MQ=255)
                                                         tAGTACCATCAGAAATCACCCCAGCGGAAGCAGACGCTTTCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCACTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATg                                                                                           <  1:1071284/150‑1 (MQ=255)
                                                           tttCCATCAGAAACCACCCCAGGGGAAGCTGACGCTATCGTTAACGTAACGCTTGATTACCTTTAATACGTTAGGGCGTTATCTGACCTGTCATATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTc                                                                                         <  1:1421216/149‑1 (MQ=255)
                                                            ttCCATCAGAAATCACCCCAGCGGAAGCAGACGCTATCGTTAACGTAACGCTGGATGACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATggcc                                                                                        <  1:514749/150‑4 (MQ=255)
                                                              ccATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTGTCTCGTTGTATCAGGTTGAAAAATAACTATGGCaaa                                                                                      >  1:107476/1‑150 (MQ=255)
                                                               cATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCaaaa                                                                                     >  1:1207910/1‑150 (MQ=255)
                                                                  cAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGa                                                                                  >  1:317512/1‑150 (MQ=255)
                                                                  cAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGCCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTTTCAGGTTGAAAAATAACTATGTCAAAACGa                                                                                  >  1:998140/1‑150 (MQ=255)
                                                                    taaaTCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAAc                                                                                <  1:151179/149‑1 (MQ=255)
                                                                        tCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCGTATTTCTCGTTGTATCAGGTTGAAAAATAACAATGTCAAAACGAACAtta                                                                            >  2:1524873/1‑149 (MQ=255)
                                                                         cACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTcg                                                                           <  2:935613/150‑1 (MQ=255)
                                                                          aCCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTcgc                                                                          <  1:498600/150‑1 (MQ=255)
                                                                           ccccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTcgcg                                                                         >  1:8412/1‑150 (MQ=255)
                                                                             ccAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGt                                                                       <  1:515987/150‑1 (MQ=255)
                                                                               aGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGa                                                                     <  2:1459057/150‑1 (MQ=255)
                                                                               aGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCGCTTTCTCGGTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTgt                                                                     >  2:561411/1‑149 (MQ=255)
                                                                                   gAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATAtt                                                                 >  2:1058649/1‑150 (MQ=255)
                                                                                    aaGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTa                                                                <  1:1453092/150‑1 (MQ=255)
                                                                                       cAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAcc                                                             >  2:1291904/1‑150 (MQ=255)
                                                                                       cAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAAcc                                                             <  1:510264/150‑1 (MQ=255)
                                                                                          aCGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttg                                                          <  1:409912/150‑1 (MQ=255)
                                                                                           cGCTATCGTTAACGTAACGCTGGATTACCGTTAATTCGTTAGGTCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgt                                                         <  1:1255364/150‑1 (MQ=255)
                                                                                             cTATCGTTAACGTAACGCTGGATTAGCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgttg                                                       <  2:1810014/150‑1 (MQ=255)
                                                                                             cTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCttgttg                                                       <  1:1830584/150‑1 (MQ=255)
                                                                                              tATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTACTCTTTATAGTTGTATCAGGTTGAAAAATAACTATGTAAAAACGAACATTCGCGGTGATATTTACCTtgttgt                                                      >  1:1615429/1‑150 (MQ=255)
                                                                                               aTCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTtgtg                                                     <  2:877569/150‑1 (MQ=255)
                                                                                               aTCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTtgtg                                                     <  2:193573/150‑1 (MQ=255)
                                                                                                tCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTgtgt                                                    <  1:670190/150‑1 (MQ=255)
                                                                                                tCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTa                                                   <  2:214481/150‑1 (MQ=255)
                                                                                                 cGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTa                                                   <  2:1110384/150‑1 (MQ=255)
                                                                                                   ttAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGc                                                 <  2:317512/150‑1 (MQ=255)
                                                                                                    taacgtaacgCTGGATTACCTTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCt                                                <  1:1354786/150‑1 (MQ=255)
                                                                                                      acgtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCATATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTc                                              <  1:1847170/150‑1 (MQ=255)
                                                                                                      acgtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTc                                              <  2:1625140/150‑1 (MQ=255)
                                                                                                       cgtaacgCTGGAGTACCGATAATACGTTAGGGCGTTAGCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTGGAAAAATAACTATGTCAAAACGTACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCt                                             <  2:27972/150‑1 (MQ=255)
                                                                                                        gtaacgCTGGATTACCGTTAATACGTTAGGGCGTTATATGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTg                                            <  2:1163806/150‑1 (MQ=255)
                                                                                                         taacgCTTGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTGATATATCGGCCTTTTCCTTCCTCTTTCTCGTTTTATCAGGTTTAAAAATAACTATTTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGt                                           <  2:1211936/150‑1 (MQ=255)
                                                                                                         taacgCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGt                                           <  1:839174/150‑1 (MQ=255)
                                                                                                            atcTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCATAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTAtt                                        <  1:1277820/148‑1 (MQ=255)
                                                                                                             gCTGTATAACCGTTAATACGTTAGGTCGTTATCTGACCTGTCAGATAACGCACTTTTCCTTCCTCTTTCTCGTTGTATCAGGCTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTg                                       <  2:155648/150‑1 (MQ=255)
                                                                                                             gCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTa                                                   >  1:1172733/1‑138 (MQ=255)
                                                                                                             gCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTg                                       <  1:1291904/150‑1 (MQ=255)
                                                                                                                   ttACCGTTTATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAgg                                 <  1:1795614/150‑1 (MQ=255)
                                                                                                                     aCCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGgcg                               >  2:208658/1‑150 (MQ=255)
                                                                                                                     aCCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTAGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGgcg                               >  2:343216/1‑150 (MQ=255)
                                                                                                                      ccGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTTTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATTTTAACATTGTTGTGTAGCTTCTGTATTGGCCAGgcgg                              >  1:1577470/1‑149 (MQ=255)
                                                                                                                             tACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCagg                       <  2:1485774/150‑1 (MQ=255)
                                                                                                                             tACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCagg                       <  2:1762306/150‑1 (MQ=255)
                                                                                                                               cGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTACTTCCTCTTTCTCGTTGTATCAAGTTGAAAAATAACTATGTCAAAACGAACATTCGCCGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCaggag                     >  1:572423/1‑150 (MQ=255)
                                                                                                                               cGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTGTTCCTACCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTCTAGCTTCTGTATTGGCCAGGCGCTTGCaggag                     >  1:303454/1‑150 (MQ=255)
                                                                                                                                      gCGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCTGGAGGAATCGt              >  1:175721/1‑150 (MQ=255)
                                                                                                                                       cGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGtt             <  1:344621/150‑1 (MQ=255)
                                                                                                                                       cGTTATCTGACAGGTCAGATAACGCCCTTTTCCTTCCTCTTTCACGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAAAATTCGCTGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGtt             >  2:1632495/1‑150 (MQ=255)
                                                                                                                                        gTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGttt            <  2:1758686/150‑1 (MQ=255)
                                                                                                                                        gTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTGGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGttt            <  2:1149129/150‑1 (MQ=255)
                                                                                                                                          tATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTa          >  1:1416360/1‑150 (MQ=255)
                                                                                                                                             cTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTGAAAACGATCATTCGCGGTGATATTATCCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAg       >  2:1231031/1‑150 (MQ=255)
                                                                                                                                              tGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGc      <  1:1300221/150‑1 (MQ=255)
                                                                                                                                               gACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACcgcg     >  2:1327734/1‑150 (MQ=255)
                                                                                                                                                  cTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAAc  >  2:559100/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGCGGAATCCAATGGTATTGTCAATGTTGACCTCATCGCTTACTACGTATCAACCGTAGTTCCATCAGAAATCACCCCAGGGGAAGCAGACGCTATCGTTAACGTAACGCTGGATTACCGTTAATACGTTAGGGCGTTATCTGACCTGTCAGATAACGCCCTTTTCCTTCCTCTTTCTCGTTGTATCAGGTTGAAAAATAACTATGTCAAAACGAACATTCGCGGTGATATTAACCTTGTTGTGTAGCTTCTGTATTGGCCAGGCGCTTGCAGGAGGAATCGTTTTACAGCGAAC  >  NC_000913/3287886‑3288181

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