Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,453,610 G→T 5.1% pseudogene (1392/2646 nt) yfcU ← pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,453,6100GT 5.1% 284.7 / 2.5 92pseudogene (1392/2646 nt)yfcUpseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (36/51);  new base T (3/2);  total (39/53)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCAGCGCCGCAGACTGATA  >  NC_000913/2453462‑2453759
                                                                                                                                                    |                                                                                                                                                     
aaGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACTGACCATATCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTGCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTgg                                                                                                                                                      <  1:341998/150‑1 (MQ=255)
aaGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCTGATACTCGCTCATGGGCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTgg                                                                                                                                                      <  1:251576/150‑1 (MQ=255)
aaGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTgg                                                                                                                                                      <  2:805425/150‑1 (MQ=255)
aaGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATACAGATACTCGATCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTgt                                                                                                                                                      >  1:735075/1‑149 (MQ=255)
 aggggCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGa                                                                                                                                                     <  1:1296366/147‑1 (MQ=255)
 aGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCATCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGa                                                                                                                                                     >  1:345090/1‑150 (MQ=255)
    gttggTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATa                                                                                                                                                  >  2:503715/4‑150 (MQ=255)
        gtgtACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAAc                                                                                                                                                >  2:1708139/1‑148 (MQ=255)
        gtgtACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAAc                                                                                                                                                <  1:1708139/148‑1 (MQ=255)
          gtACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGAGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTcac                                                                                                                                            >  2:913759/1‑150 (MQ=255)
           tACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTcaca                                                                                                                                           <  1:720553/150‑1 (MQ=255)
               tctcTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGa                                                                                                                                       >  2:801138/1‑150 (MQ=255)
                ctctTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGaa                                                                                                                                      >  2:1170293/1‑150 (MQ=255)
                 tctTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAAAAACTCACACGGAAg                                                                                                                                     >  1:244842/1‑150 (MQ=255)
                    ttGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAg                                                                                                                                  <  1:1485091/150‑1 (MQ=255)
                    ttGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAg                                                                                                                                  <  1:472450/150‑1 (MQ=255)
                     tGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGt                                                                                                                                 <  1:137614/150‑1 (MQ=255)
                        cGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTAc                                                                                                                              >  2:1685830/1‑150 (MQ=255)
                           tGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGt                                                                                                                           <  1:903804/150‑1 (MQ=255)
                            gCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGGCGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTc                                                                                                                          >  1:1543829/1‑150 (MQ=255)
                             cccGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTcc                                                                                                                         <  1:1387747/150‑1 (MQ=255)
                             cccGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTcc                                                                                                                         >  1:515023/1‑150 (MQ=255)
                                gTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAATTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGc                                                                                                                      >  2:1420560/1‑150 (MQ=255)
                                gTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGc                                                                                                                      <  1:801138/150‑1 (MQ=255)
                                   cGGACTATTTCGCTGTCACTGGCATCCAGATACTCGCTAATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGaa                                                                                                                   <  2:1148373/150‑1 (MQ=255)
                                     gACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAAcc                                                                                                                 >  2:534107/1‑150 (MQ=255)
                                      aCCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCt                                                                                                                >  1:126936/1‑150 (MQ=255)
                                       ccATTTCGCTGTCACTGGCATCCAGATACTCGCGCATTGTCATAACGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGGTACCGTCCAGCGAACCtt                                                                                                               <  1:289623/150‑1 (MQ=255)
                                         aTTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTATCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTTTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTg                                                                                                             <  2:1166274/150‑1 (MQ=255)
                                          tttCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGc                                                                                                            <  1:1593874/150‑1 (MQ=255)
                                           ttCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGcc                                                                                                           <  2:400511/150‑1 (MQ=255)
                                              gCTGTCACTGGCATCCAGATACTCGCTCTTGGTCATAAAGTTCTCTTACGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATa                                                                                                        <  2:132679/150‑1 (MQ=255)
                                              gCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATa                                                                                                        >  2:1176717/1‑150 (MQ=255)
                                              gCTGTCACTGGCATCCAGATACTCGCTCATGGTCAAAAAGTTATCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATa                                                                                                        <  2:1265439/150‑1 (MQ=255)
                                                  tCACTGGCATCCAGATACTCGCTCATTGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCg                                                                                                    <  1:1550302/150‑1 (MQ=255)
                                                  tCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCg                                                                                                    <  1:176504/150‑1 (MQ=255)
                                                  tCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCg                                                                                                    <  2:126936/150‑1 (MQ=255)
                                                    ccTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGGGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGgt                                                                                                  <  1:113601/149‑1 (MQ=255)
                                                     cTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGATGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGgtg                                                                                                 >  2:367712/1‑150 (MQ=255)
                                                     cTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGgtg                                                                                                 <  1:304811/150‑1 (MQ=255)
                                                       ggCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTc                                                                                               <  2:962715/150‑1 (MQ=255)
                                                         cATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCtt                                                                                             >  2:1312091/1‑150 (MQ=255)
                                                           tCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCTAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTa                                                                                           <  2:735075/150‑1 (MQ=255)
                                                           tCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTa                                                                                           <  1:969178/150‑1 (MQ=255)
                                                           cccAGATAGTCGCTCATGGTCATAAAGTGCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGTTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTa                                                                                           <  2:1067597/149‑1 (MQ=255)
                                                            ccAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTAt                                                                                          >  1:786491/1‑150 (MQ=255)
                                                             cAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATc                                                                                         >  1:1549937/1‑150 (MQ=255)
                                                              aGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATcc                                                                                        <  1:1021444/150‑1 (MQ=255)
                                                                    tCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTg                                                                                  <  1:1382225/150‑1 (MQ=255)
                                                                     cGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTgg                                                                                 <  1:1170293/150‑1 (MQ=255)
                                                                         cATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATg                                                                             <  1:1640994/150‑1 (MQ=255)
                                                                           ttgTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCg                                                                           <  1:1617137/148‑1 (MQ=255)
                                                                           tGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTGCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCg                                                                           <  1:963123/150‑1 (MQ=255)
                                                                             gTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTTGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGaa                                                                         >  1:400505/1‑150 (MQ=255)
                                                                             gTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTTGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGaa                                                                         >  1:1680147/1‑150 (MQ=255)
                                                                             gTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGaa                                                                         >  2:814631/1‑150 (MQ=255)
                                                                              tCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAAt                                                                        <  2:515023/150‑1 (MQ=255)
                                                                               cATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATg                                                                       <  1:1274498/150‑1 (MQ=255)
                                                                               cATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATg                                                                       <  2:1169880/150‑1 (MQ=255)
                                                                                 tAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAg                                                                                      >  2:859851/1‑133 (MQ=255)
                                                                                 tAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAg                                                                                      <  1:859851/133‑1 (MQ=255)
                                                                                 tAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAg                                                                     >  2:955293/1‑150 (MQ=255)
                                                                                   aaGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTa                                                                   <  1:1685830/150‑1 (MQ=255)
                                                                                   aaGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTa                                                                   >  1:423806/1‑150 (MQ=255)
                                                                                     tttCTCTTCCGAGAAGCGATATCCAGCGAAAGGGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAAc                                                                 <  1:661852/149‑1 (MQ=255)
                                                                                     gTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGTTTGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAAc                                                                 <  1:28364/150‑1 (MQ=255)
                                                                                       tctcTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACAt                                                               <  2:734125/150‑1 (MQ=255)
                                                                                          ctTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCaa                                                            >  2:1135560/1‑150 (MQ=255)
                                                                                             ccGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAAcg                                                         <  2:1208583/150‑1 (MQ=255)
                                                                                             ccGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAAcg                                                         <  1:1261157/150‑1 (MQ=255)
                                                                                              cGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAAcgc                                                        <  1:503715/150‑1 (MQ=255)
                                                                                                 gaAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCTTAAGAGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGAc                                                     <  2:1358654/150‑1 (MQ=255)
                                                                                                     cGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTTGAATAAATCACACGGAACGAGTTACCGTCCAGCGAAACTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTTACATCAAACGCGACcgcc                                                 >  1:1369155/1‑149 (MQ=255)
                                                                                                     cGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACcgcg                                                 >  1:1717708/1‑150 (MQ=255)
                                                                                                       atatCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGcc                                               <  1:534107/150‑1 (MQ=255)
                                                                                                           ccAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAg                                           >  1:804531/1‑150 (MQ=255)
                                                                                                            cAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAgg                                          >  2:1401693/1‑150 (MQ=255)
                                                                                                                 aaaGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGGCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGAc                                     <  2:843963/150‑1 (MQ=255)
                                                                                                                  aaGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACa                                    <  1:803358/150‑1 (MQ=255)
                                                                                                                   aGTGACGCGGCTGTTGCGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACACCAAACGCGACCGCGCCGAAGGTAGACaa                                   <  1:173495/150‑1 (MQ=255)
                                                                                                                   aGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACaa                                   >  1:616788/1‑150 (MQ=255)
                                                                                                                     tGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTATCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAAt                                 >  1:256829/1‑150 (MQ=255)
                                                                                                                             cTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGcgaccg                         >  1:111626/1‑150 (MQ=255)
                                                                                                                              tGTTGAGCTGGTCAAAGTCTTTTGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGcgaccga                        <  2:1680147/150‑1 (MQ=255)
                                                                                                                              tGTTGAGCTGGTCAAAGTCTTTTGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGcgaccga                        <  2:400505/150‑1 (MQ=255)
                                                                                                                              tGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGcgaccga                        <  2:721574/150‑1 (MQ=255)
                                                                                                                               gTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGcgaccgac                       >  1:645379/1‑150 (MQ=255)
                                                                                                                                    gCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCa                  >  1:1420339/1‑150 (MQ=255)
                                                                                                                                      tGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAATCGGTGTCTTTATCCAGTTTGGTATGAGAATGAGTAACTTCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCACAgc                >  1:956431/1‑150 (MQ=255)
                                                                                                                                      tGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCAgc                >  1:963018/1‑150 (MQ=255)
                                                                                                                                        gTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCAgcgc              >  2:128044/1‑150 (MQ=255)
                                                                                                                                           aatggcTTTGGAATAACTCACAGGGAAGTTGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGTATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCAGcgccgc           <  2:345090/145‑1 (MQ=255)
                                                                                                                                            aaGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCAGCGCCGCa          >  2:400434/1‑150 (MQ=255)
                                                                                                                                            aaGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCAGCGCCGCa          <  2:1543829/150‑1 (MQ=255)
                                                                                                                                             aGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCAGCGCCGCAg         >  2:663149/1‑150 (MQ=255)
                                                                                                                                               tCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCAGCGCCGCAGAc       <  2:1549937/150‑1 (MQ=255)
                                                                                                                                                    ggAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCAGCGCCGCAGACTGATa  >  2:900813/1‑150 (MQ=255)
                                                                                                                                                    ggAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCAGCGCCGCAGACTGATa  <  1:955293/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGTGGCGGTGTACATCTCTTTGTCGTTGCCCGTACGGACCATTTCGCTGTCACTGGCATCCAGATACTCGCTCATGGTCATAAAGTTCTCTTCCGAGAAGCGATATCCAGCGAAAGTGACGCGGCTGTTGAGCTGGTCAAAGTCTTTGGAATAACTCACACGGAAGGAGTTACCGTCCAGCGAACCTTTGCCATAAGCGGTGTCTTTATCCAGTTTGGTATGCGAATGAGTAACATCAAACGCGACCGCGCCGAAGGTAGACAAATCGCGACCGACGCCCAGCGCCGCAGACTGATA  >  NC_000913/2453462‑2453759

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