Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 314,578 G→A 5.3% pseudogene (222/872 nt) eaeH → pseudogene, attaching and effacing protein homology;factor; Not classified

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913314,5780GA 5.3% 365.8 / 5.9 115pseudogene (222/872 nt)eaeHpseudogene, attaching and effacing protein homology;factor; Not classified
Reads supporting (aligned to +/- strand):  ref base G (56/52);  new base A (2/4);  total (59/56)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 116 total aligned reads displayed.

ATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAAT  >  NC_000913/314433‑314725
                                                                                                                                                 |                                                                                                                                                   
atatCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACGAAAGTAATGGCGGCACTTGCTAAGCAAGCGGTTCAGCCAAGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTg                                                                                                                                                 >  1:290248/1‑150 (MQ=255)
atatCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTg                                                                                                                                                 <  1:398530/150‑1 (MQ=255)
 tatCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGc                                                                                                                                                <  2:600070/150‑1 (MQ=255)
  atCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGTTTGAGCATGGGAAATACAACGGTAACTGCTGATAATAACGAGGAGAAAAAAGTCGCGTCCTTTccc                                                                                                                                               >  1:95455/1‑150 (MQ=255)
  atCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGcc                                                                                                                                               <  1:36785/150‑1 (MQ=255)
    ctctGTTCAGGTTCTTTTTCCACTCGCTGTCACGTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGTTAACTGCTGATAATAACGTGGTGAAAAATGTCCCGGCATTTGCGGc                                                                                                                                             <  1:779121/150‑1 (MQ=255)
    ctctGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGc                                                                                                                                             >  1:1691579/1‑150 (MQ=255)
       tGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCTTGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCaaa                                                                                                                                          >  2:1685611/1‑150 (MQ=255)
       tGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCaaa                                                                                                                                          <  1:1557206/150‑1 (MQ=255)
        gTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCCCGGTAGAGCATGGGAAATACTACGGTAACTGATGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAAt                                                                                                                                         >  1:888145/1‑150 (MQ=255)
        gTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAAt                                                                                                                                         <  1:1307370/150‑1 (MQ=255)
           cAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGGAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGcc                                                                                                                                      >  2:1137634/1‑150 (MQ=255)
             ggTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGTTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCgg                                                                                                                                    >  2:856675/1‑150 (MQ=255)
             ggTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCgg                                                                                                                                    <  2:36424/150‑1 (MQ=255)
              gTTCATTTTCCACTCGCTGTCTCCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGATCAGCCACGGTTGTTCATGGGATATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCggg                                                                                                                                   <  2:1451983/150‑1 (MQ=255)
                   ttttCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACAtt                                                                                                                              >  2:463649/1‑150 (MQ=255)
                     ttCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACAtttt                                                                                                                            >  1:759966/1‑150 (MQ=255)
                      tCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACAttttt                                                                                                                           >  2:653100/1‑150 (MQ=255)
                      tCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACAttttt                                                                                                                           >  2:685846/1‑150 (MQ=255)
                      tCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACAttttt                                                                                                                           <  2:1024896/150‑1 (MQ=255)
                       ccACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTa                                                                                                                          >  1:1383580/1‑150 (MQ=255)
                       ccACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTa                                                                                                                          >  2:393450/1‑150 (MQ=255)
                        cACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGATCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCCCGGCATTTGCGGCAAACGACGCGACAtttttct                                                                                                                         <  2:897874/150‑3 (MQ=38)
                         aCTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCCTGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGATTGCCGCAAATGCCGGGACATTTTTAAg                                                                                                                        >  2:1155608/1‑150 (MQ=255)
                         aCTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTTGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAg                                                                                                                        <  2:192334/150‑1 (MQ=255)
                          cTCGCTGTCACCTTTACCCCAGTAATGGCGGCACTTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCCCGGCATTTGCGGCAAACGACGCGACAtttttctcc                                                                                                                       <  1:1530856/150‑5 (MQ=37)
                          cTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCCCGGCATTTGCGGCAAACGACGCGACAtttttctcc                                                                                                                       <  2:498371/150‑5 (MQ=38)
                            cGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAg                                                                                                                     <  2:1205176/150‑1 (MQ=255)
                            cGCTCTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATCTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAg                                                                                                                     <  1:1024809/150‑1 (MQ=255)
                             gCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGt                                                                                                                    >  1:899062/1‑150 (MQ=255)
                              cTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTc                                                                                                                   >  1:1504915/1‑150 (MQ=255)
                                  cACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGcc                                                                                                               <  2:1666338/150‑1 (MQ=255)
                                          ccccAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCg                                                                                                       >  1:1697758/1‑150 (MQ=255)
                                          ccccAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCg                                                                                                       >  1:1527251/1‑150 (MQ=255)
                                           cccAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGa                                                                                                      <  1:193885/150‑1 (MQ=255)
                                            ccAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGAt                                                                                                     >  1:1192476/1‑150 (MQ=255)
                                                 aaTGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAACTGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGac                                                                                                >  1:1066559/1‑150 (MQ=255)
                                                 aaTGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGac                                                                                                >  1:865840/1‑150 (MQ=255)
                                                 aaTGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGac                                                                                                <  1:790823/150‑1 (MQ=255)
                                                 aaTGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGac                                                                                                <  1:890971/150‑1 (MQ=255)
                                                    ggcggcAGGTGCGCAGCATGCGGGGCAGCCACGGGTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGGCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGAGGCGACACg                                                                                             <  2:95847/150‑1 (MQ=255)
                                                         cACGTGCGCAGCAGGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACt                                                                                        <  2:918472/150‑1 (MQ=255)
                                                            gTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAGTACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTtta                                                                                     >  1:401185/1‑150 (MQ=255)
                                                              gcgcAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTttatt                                                                                   <  1:1717274/150‑1 (MQ=255)
                                                                 cagcaTGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGTTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGTACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTAcc                                                                                <  1:1070132/150‑1 (MQ=255)
                                                                    caTGCGGTTCGGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAACCAGTCCGCCAGATAGCGATGCGACACGTAACTTTATTACAGGa                                                                             >  2:1366445/1‑150 (MQ=255)
                                                                    caTGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGa                                                                             <  1:1238414/150‑1 (MQ=255)
                                                                    caTGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGa                                                                             >  2:169948/1‑150 (MQ=255)
                                                                    caTGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGa                                                                             >  2:237647/1‑150 (MQ=255)
                                                                      tGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAAt                                                                           >  2:1242267/1‑150 (MQ=255)
                                                                       gCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATg                                                                          <  1:622998/150‑1 (MQ=255)
                                                                         ggTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGc                                                                        <  1:591166/150‑1 (MQ=255)
                                                                          gTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGcc                                                                       <  2:1294323/150‑1 (MQ=255)
                                                                             cAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCcaca                                                                    <  1:793859/150‑1 (MQ=255)
                                                                              aGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACag                                                                   >  2:1005916/1‑150 (MQ=255)
                                                                               gCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGCTGCGACACGTAACTGTATTACCGGAATGGCCACagc                                                                  >  2:967709/1‑150 (MQ=255)
                                                                               gCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCTACACGTAACTTTATTACCGGAATGGCcacatc                                                                  >  2:974825/1‑148 (MQ=255)
                                                                               gCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACagc                                                                  >  1:1466668/1‑150 (MQ=255)
                                                                               gCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACagc                                                                  >  2:721284/1‑150 (MQ=255)
                                                                               gCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGAAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACAGGAATGGCCACagc                                                                  >  2:945901/1‑150 (MQ=255)
                                                                                    gcttGAGCATGGGAAATACTACGGTAACTGCAGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACagctaaag                                                             <  2:873196/148‑1 (MQ=255)
                                                                                     gTTGAGCATGGGATATCATACGGTAACTGCTGATAATAAAGTGGAGAAAAATGTCGCGGCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAAGTTTATTACCGGAATGGCCACagctaaagc                                                            <  1:281082/150‑1 (MQ=255)
                                                                                         aGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAAc                                                        <  2:1571561/150‑1 (MQ=255)
                                                                                             tGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATCCGCCCCGTAACATTATTACCGGAATGGCCACAGCTAAAGCTAACCAgg                                                    >  2:1204157/1‑150 (MQ=255)
                                                                                              gggAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCGACAGCTAATGCTTACCAGGa                                                   >  1:1149635/1‑150 (MQ=255)
                                                                                              gggAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGa                                                   >  1:377723/1‑150 (MQ=255)
                                                                                              gggAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGa                                                   >  2:1727175/1‑150 (MQ=255)
                                                                                                gAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGaaa                                                 >  1:1484316/1‑150 (MQ=255)
                                                                                                   aTACTACGGTAACTGCTGATAATCACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGCGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATAc                                              <  2:348101/150‑1 (MQ=255)
                                                                                                   aTACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAATCAGGAAATAc                                              <  1:653100/150‑1 (MQ=255)
                                                                                                   aTACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATAc                                              <  1:685846/150‑1 (MQ=255)
                                                                                                      ctacGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCAAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACAAGGAAATACAgg                                           >  2:729671/1‑150 (MQ=255)
                                                                                                      ctaaGGTAACTGCTGATACAAACGTGGAGAATAATGTCGCGTCGTTTGCCGGAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGACCACAGCTAAAGCTAACCAGGAAATACAgg                                           <  2:720819/150‑1 (MQ=255)
                                                                                                        acGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAg                                         <  1:1643751/150‑1 (MQ=255)
                                                                                                          ggTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAATTACAGGAGTg                                       >  2:801690/1‑150 (MQ=255)
                                                                                                          ggTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTg                                       >  2:1541479/1‑150 (MQ=255)
                                                                                                          ggTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTg                                       <  1:635622/150‑1 (MQ=255)
                                                                                                                    gATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGCTAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGaaa                             >  1:947514/1‑150 (MQ=255)
                                                                                                                     ataataACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAat                            >  2:491545/1‑150 (MQ=255)
                                                                                                                     ataataACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAat                            >  1:1680814/1‑150 (MQ=255)
                                                                                                                      taataaCGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAata                           >  1:1641279/1‑150 (MQ=255)
                                                                                                                       aataaCGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAatat                          <  2:1273248/150‑1 (MQ=255)
                                                                                                                        ataaCGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATg                         <  1:393450/150‑1 (MQ=255)
                                                                                                                         taaCGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTAGCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATgg                        <  1:1394558/150‑1 (MQ=255)
                                                                                                                         taaCGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATgg                        <  1:721284/150‑1 (MQ=255)
                                                                                                                         taaCGTGGAGAAAAATGTCGAGTCGTTTGCCGCAAATGCCGTGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACTCGTATCTTTATTTCCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATgg                        <  2:95455/150‑1 (MQ=255)
                                                                                                                            cGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTAc                     <  2:1383580/150‑1 (MQ=255)
                                                                                                                              tGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTg                   <  2:1527251/150‑1 (MQ=255)
                                                                                                                               ggAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTgc                  <  2:552693/150‑1 (MQ=255)
                                                                                                                               ggAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACAAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTgc                  >  2:967061/1‑150 (MQ=255)
                                                                                                                               gcttaaaaaTGTCCCGGCATTTGCGGCAAACGACGCGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTgc                  >  1:1241122/5‑150 (MQ=38)
                                                                                                                               gcttaaaaaTGTCCCGGCATTTGCGGCAAACGACGCGACAGTTTTAAGCAGTCAGCCAGAGAGCGATGCGACACGTAACTTTATTCCCGGAATGGCCACAGCTAAAGATAACCAGGACATACAGGAGTGGCTCGTGTAATATGGTCCTgc                  >  1:81046/5‑150 (MQ=25)
                                                                                                                                 agaAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACAAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTgcgc                >  2:1045000/1‑150 (MQ=255)
                                                                                                                                 agaAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCAAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATTGTACTgcgc                >  2:158980/1‑150 (MQ=255)
                                                                                                                                 agaAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTGAAGCAGTCAGCCAGATAGCGATGCGTCACGTAACTTTATTACCGGAATGGCCCCAGCTTAGGCTATCGAGGAAATACAGGAGTGGCTCGGGAAATAGGGTACTgcgc                >  1:318573/1‑150 (MQ=255)
                                                                                                                                     aaaTGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGTCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGt            <  2:888145/150‑1 (MQ=255)
                                                                                                                                         gTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAACCAGTCAGCCAGATAGCGATGCGACAAGTAACTTTATTACCGGAATGGCCACACCAAAAGCTAACCAGGAAATACAGGAGTCGGTCGGGAAATATGGTAATGCGAGCGACaaa        >  2:124225/1‑150 (MQ=255)
                                                                                                                                           cgcgTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACt      <  2:1477720/150‑1 (MQ=255)
                                                                                                                                            gcgTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAAccg     >  1:372720/1‑148 (MQ=255)
                                                                                                                                               tcgtTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAAt  <  1:686200/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ATATCTCTGTTCAGGTTCTTTTTCCACTCGCTGTCACCTTTACCCCAGTAATGGCGGCACGTGCGCAGCATGCGGTTCAGCCACGGTTGAGCATGGGAAATACTACGGTAACTGCTGATAATAACGTGGAGAAAAATGTCGCGTCGTTTGCCGCAAATGCCGGGACATTTTTAAGCAGTCAGCCAGATAGCGATGCGACACGTAACTTTATTACCGGAATGGCCACAGCTAAAGCTAACCAGGAAATACAGGAGTGGCTCGGGAAATATGGTACTGCGCGCGTCAAACTGAAT  >  NC_000913/314433‑314725

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: