Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 852,372 A→G 5.4% F76L (TTC→CTC)  opgE ← OPG biosynthetic transmembrane phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913852,3720AG 5.4% 362.3 / 8.0 110F76L (TTC→CTC) opgEOPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base A (38/65);  new base G (3/3);  total (41/69)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 111 total aligned reads displayed.

TTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGTTAATTAATAGCGAC  >  NC_000913/852225‑852520
                                                                                                                                                   |                                                                                                                                                    
ttACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAGAt                                                                                                                                                    >  1:2169186/1‑150 (MQ=255)
ttACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAGAt                                                                                                                                                    >  1:970342/1‑150 (MQ=255)
ttACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAGAt                                                                                                                                                    >  1:1455942/1‑150 (MQ=255)
  aCGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATaa                                                                                                                                                  <  2:914595/150‑1 (MQ=255)
      cGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAg                                                                                                                                              >  1:360779/1‑150 (MQ=255)
      cGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAg                                                                                                                                              <  2:2010030/150‑1 (MQ=255)
          cgcAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAgccgc                                                                                                                                          <  2:7204/150‑1 (MQ=255)
          cgcAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAgccgc                                                                                                                                          >  2:2063355/1‑150 (MQ=255)
             aaaaaTAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGccac                                                                                                                                       >  2:1262585/1‑150 (MQ=255)
             aaaaaTAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGccac                                                                                                                                       >  2:1882692/1‑150 (MQ=255)
             aaaaaTAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCAGAGCAAAAGGGAAATAACAAGACGccac                                                                                                                                       >  2:1969812/1‑150 (MQ=255)
               aaataaatAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGccacca                                                                                                                                     >  1:82836/1‑150 (MQ=255)
               aaataaatAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGccacca                                                                                                                                     >  1:82833/1‑150 (MQ=255)
                    aataaGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGaa                                                                                                                                >  2:179362/1‑150 (MQ=255)
                         gCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGac                                                                                                                           >  1:1068498/1‑150 (MQ=255)
                         gCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGac                                                                                                                           >  1:1213819/1‑150 (MQ=255)
                               ctggTTGCTGTACCACGTAAAAATCTTCTGGATTTGGTTCGACTCTTCCTTTTTGGTGGTGTGCAATGCCAGAATTTTTTTTAAATTTGGCGCGCCTTAGGCCTGAGCAAAAGGTAAATAACAAGCCGCCACCAGCGAACTGACACCGAc                                                                                                                     <  1:7162/148‑1 (MQ=255)
                               ctggTTGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATTTTTTGAATTTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGAc                                                                                                                     <  1:1235035/148‑1 (MQ=255)
                               cAGGTATCTGTACCCCGGAAAAATCGGCCGGATTTCGGTCGACTCTTCCATCTTGGTGGAGTGAAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTTGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGAc                                                                                                                     <  1:2352860/150‑1 (MQ=255)
                                  gTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAg                                                                                                                  >  1:2357750/1‑150 (MQ=255)
                                   gaGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGa                                                                                                                 <  2:2238455/149‑1 (MQ=255)
                                       tGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACtt                                                                                                             >  1:528312/1‑150 (MQ=255)
                                         tACCACGTAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACtttt                                                                                                           <  1:740159/150‑1 (MQ=255)
                                         tACCACGGAAAAATCGTCAGGATTTCGGTCTACTATTCCATATTGTTGGAGTGCAATGCCAGCAATGTATTGAATTTAGGCTCTCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACtttt                                                                                                           <  1:1247389/150‑1 (MQ=255)
                                         tACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACtttt                                                                                                           <  1:146300/150‑1 (MQ=255)
                                         tACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACtttt                                                                                                           <  1:1540392/150‑1 (MQ=255)
                                            cACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGt                                                                                                        >  1:572221/1‑150 (MQ=255)
                                                gAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACAGACCAGAACTTTTTGTAcgc                                                                                                    <  1:302893/150‑1 (MQ=255)
                                                     aTCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGc                                                                                               <  2:663042/150‑1 (MQ=255)
                                                        gTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAAt                                                                                            <  2:360779/150‑1 (MQ=255)
                                                         tCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATg                                                                                           <  1:1878212/150‑1 (MQ=255)
                                                         tCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATg                                                                                           <  1:1101059/150‑1 (MQ=255)
                                                         tCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATg                                                                                           <  1:2396106/150‑1 (MQ=255)
                                                              aTTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGc                                                                                      <  2:1582532/150‑1 (MQ=255)
                                                              aTTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGc                                                                                      <  2:1478310/150‑1 (MQ=255)
                                                              aTTTCGGTCGACTCTTCCATATTGGTGGAGGGCAATGCCAGCAATGTAGTGAAAGTAGGCGCGCCGTAGGCCTGAGCAACAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCCGAACTTTTTGTACGCGAGGCAATGTTCGc                                                                                      <  2:1249825/150‑1 (MQ=255)
                                                                ttCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCa                                                                                    <  1:1959764/150‑1 (MQ=255)
                                                                ttCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCa                                                                                    <  1:1959765/150‑1 (MQ=255)
                                                                  cGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAACCGTGACCGCCCCCGCCGCCACCAGCGAACTGGCACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATa                                                                                  <  2:2104849/150‑1 (MQ=255)
                                                                  cGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGGCACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATa                                                                                  <  2:440892/150‑1 (MQ=255)
                                                                  cGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGGCACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATa                                                                                  <  2:186523/150‑1 (MQ=255)
                                                                    gTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATagc                                                                                >  1:887275/1‑150 (MQ=255)
                                                                        aCTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATagcaaaa                                                                            <  1:2119733/150‑1 (MQ=255)
                                                                            ttCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATagcaaaagcaa                                                                        <  1:596223/150‑1 (MQ=255)
                                                                            ttCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATagcaaaagcaa                                                                        <  1:618371/150‑1 (MQ=255)
                                                                             tCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATagcaaaagcaaa                                                                       >  1:632441/1‑150 (MQ=255)
                                                                               cATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAAt                                                                     >  1:1733720/1‑150 (MQ=255)
                                                                                 tatTGGTGGAGTGCTATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAg                                                                   <  2:1480198/150‑1 (MQ=255)
                                                                                   ttGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGcc                                                                 <  1:99504/150‑1 (MQ=255)
                                                                                   ttGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGcc                                                                 <  2:1046568/150‑1 (MQ=255)
                                                                                   ttGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGcc                                                                 <  2:82833/150‑1 (MQ=255)
                                                                                   ttGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGcc                                                                 <  2:82836/150‑1 (MQ=255)
                                                                                       tggAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTaa                                                             <  1:604767/150‑1 (MQ=255)
                                                                                        ggAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAGATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTaaa                                                            <  2:970342/150‑1 (MQ=255)
                                                                                        ggAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAGATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTaaa                                                            <  2:2169186/150‑1 (MQ=255)
                                                                                        ggAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAGATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTaaa                                                            <  2:1455942/150‑1 (MQ=255)
                                                                                        ggAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTaaa                                                            >  1:1064528/1‑150 (MQ=255)
                                                                                        ggAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTaaa                                                            >  2:2130652/1‑150 (MQ=255)
                                                                                        ggAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTaaa                                                            >  2:2130656/1‑150 (MQ=255)
                                                                                        ggAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTaaa                                                            >  1:1127358/1‑150 (MQ=255)
                                                                                            tGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCa                                                        <  1:168577/150‑1 (MQ=255)
                                                                                            tGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCa                                                        >  2:923928/1‑150 (MQ=255)
                                                                                            tGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCa                                                        <  1:274096/150‑1 (MQ=255)
                                                                                            tGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCa                                                        <  1:1888959/150‑1 (MQ=255)
                                                                                            tGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAAAGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCa                                                        >  2:923797/1‑150 (MQ=255)
                                                                                              cAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGc                                                      <  1:2347600/150‑1 (MQ=255)
                                                                                                     aGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTag                                               <  1:1262999/150‑1 (MQ=255)
                                                                                                         aTGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCa                                           <  1:45973/150‑1 (MQ=255)
                                                                                                         aTGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCa                                           <  1:672365/150‑1 (MQ=255)
                                                                                                           gTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCaaa                                         <  1:1407072/150‑1 (MQ=255)
                                                                                                           gTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCaaa                                         <  1:1407075/150‑1 (MQ=255)
                                                                                                              ttGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTg                                      >  1:1523719/1‑150 (MQ=255)
                                                                                                               tGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGa                                     >  2:1698010/1‑150 (MQ=255)
                                                                                                                 aaaTTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAg                                   >  1:1987638/1‑150 (MQ=255)
                                                                                                                 aaaTTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAg                                   <  1:1647668/150‑1 (MQ=255)
                                                                                                                   aTTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAggg                                 <  1:56333/150‑1 (MQ=255)
                                                                                                                        gcgcgcCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAAc                            <  1:139594/150‑1 (MQ=255)
                                                                                                                        gcgcgcCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAAc                            >  2:1878342/1‑150 (MQ=255)
                                                                                                                        gcgcgcCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAAc                            >  2:2013127/1‑150 (MQ=255)
                                                                                                                        gcgcgcCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAAc                            <  1:1660334/150‑1 (MQ=255)
                                                                                                                            gcCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTaa                        >  1:2106514/1‑150 (MQ=255)
                                                                                                                            gcCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTaa                        >  1:1424331/1‑150 (MQ=255)
                                                                                                                            gcCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTaa                        >  1:956145/1‑150 (MQ=255)
                                                                                                                                 aGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGa                   <  1:1201907/150‑1 (MQ=255)
                                                                                                                                 aGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGa                   <  1:1697988/150‑1 (MQ=255)
                                                                                                                                 aGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGa                   <  1:1697980/150‑1 (MQ=255)
                                                                                                                                 aGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGa                   <  1:906050/150‑1 (MQ=255)
                                                                                                                                   gCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAgg                 <  1:1698010/150‑1 (MQ=255)
                                                                                                                                   gCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAgg                 >  2:1423317/1‑150 (MQ=255)
                                                                                                                                      tGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGtta              <  2:465090/150‑1 (MQ=255)
                                                                                                                                      tGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGtta              <  2:1380310/150‑1 (MQ=255)
                                                                                                                                      tGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGtta              >  1:976716/1‑150 (MQ=255)
                                                                                                                                        aGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGttaat            <  2:944924/150‑1 (MQ=255)
                                                                                                                                            aaaGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGTtaattaat        <  1:322379/150‑1 (MQ=255)
                                                                                                                                             aaGGGAAATAACAAGCCGCCACCAGCTAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGAAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGTTAATTAATa       >  1:968337/1‑150 (MQ=255)
                                                                                                                                             aaGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGTTAATTAATa       <  2:1538825/150‑1 (MQ=255)
                                                                                                                                             aaGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGTTAATTAATa       <  1:923728/150‑1 (MQ=255)
                                                                                                                                              aGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTATGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGTTAATTAATAg      <  2:1213819/150‑1 (MQ=255)
                                                                                                                                               gggAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGTTAATTAATAGc     <  1:1845804/150‑1 (MQ=255)
                                                                                                                                                  aaaTAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGTTAATTAATAGCGAc  >  1:1204386/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCGTAAACGCAGCGGTGTAGAGCAAACTGAAGGGGTAACCTAAGCCGAGGTTAATTAATAGCGAC  >  NC_000913/852225‑852520

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