Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,507,500 G→T 7.8% K240N (AAG→AAT ydcP → putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,507,5000GT 7.8% 364.3 / 16.7 102K240N (AAG→AATydcPputative peptidase
Reads supporting (aligned to +/- strand):  ref base G (47/47);  new base T (4/4);  total (51/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 103 total aligned reads displayed.

TTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGAC  >  NC_000913/1507353‑1507643
                                                                                                                                                   |                                                                                                                                               
ttGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAg                                                                                                                                                 >  2:872194/1‑148 (MQ=255)
     cGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                            <  3:1018558/148‑1 (MQ=255)
     cGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                            <  2:499189/148‑1 (MQ=255)
           cGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCg                                                                                                                                      <  1:872194/148‑1 (MQ=255)
                cGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGaa                                                                                                                                           >  1:1100057/1‑138 (MQ=255)
                cGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGaa                                                                                                                                           <  2:1100057/138‑1 (MQ=255)
                   ttGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAg                                                                                                                              <  1:781580/148‑1 (MQ=255)
                           cacaTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                            <  3:391055/126‑1 (MQ=255)
                           cacaTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                            >  4:391055/1‑126 (MQ=255)
                               ttGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAgg                                                                                                                                         >  2:904823/1‑125 (MQ=255)
                               ttGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAgg                                                                                                                                         <  1:904823/125‑1 (MQ=255)
                                  aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                   <  1:896859/127‑7 (MQ=255)
                                  aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                   >  2:896859/1‑121 (MQ=255)
                                  aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                   <  3:387115/127‑7 (MQ=255)
                                  aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                   >  4:387115/1‑121 (MQ=255)
                                   aaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                       <  1:183985/139‑1 (MQ=255)
                                   aaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                       >  2:183985/1‑139 (MQ=255)
                                              gggCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagc                                                                                                                                      >  4:788636/1‑109 (MQ=255)
                                              gggCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagc                                                                                                                                      <  3:788636/112‑4 (MQ=255)
                                                          tCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACaa                                                                                                                                <  3:545203/107‑1 (MQ=255)
                                                          tCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACaa                                                                                                                                >  4:545203/1‑107 (MQ=255)
                                                           ccTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGat                                                                                                                            >  2:979422/1‑110 (MQ=255)
                                                           ccTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGat                                                                                                                            <  1:979422/110‑1 (MQ=255)
                                                             tATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGa                                                                                    >  4:161892/1‑148 (MQ=255)
                                                              aTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAggg                                                                                                                                        >  3:788030/1‑95 (MQ=255)
                                                              aTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAggg                                                                                                                                        <  4:788030/95‑1 (MQ=255)
                                                              aTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaagaa                                                                                                           <  1:285331/124‑1 (MQ=255)
                                                              aTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaagaa                                                                                                           >  2:285331/1‑124 (MQ=255)
                                                               tGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                              <  1:165876/120‑1 (MQ=255)
                                                                 aaaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAATCTTGAAgg                                                                                                                                         >  4:87318/1‑91 (MQ=255)
                                                                 aaaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAATCTTGAAgg                                                                                                                                         <  3:87318/91‑1 (MQ=255)
                                                                  aaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGc                                                                                        <  3:573483/139‑1 (MQ=255)
                                                                  aaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGc                                                                                        >  4:573483/1‑139 (MQ=255)
                                                                    aaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGa                                                                             >  4:1012206/1‑148 (MQ=255)
                                                                    aaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGa                                                                             <  3:161892/148‑1 (MQ=255)
                                                                        tCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGa                                                                                                                        >  4:885824/1‑101 (MQ=255)
                                                                        tCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGa                                                                                                                        <  3:885824/101‑1 (MQ=255)
                                                                           gctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                     <  3:11041/85‑1 (MQ=255)
                                                                           gctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                     >  4:11041/1‑85 (MQ=255)
                                                                              gTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaaga                                                                                                            <  2:981050/107‑1 (MQ=255)
                                                                              gTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaaga                                                                                                            >  1:981050/1‑107 (MQ=255)
                                                                                  aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                            <  1:592837/71‑1 (MQ=255)
                                                                                  aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                            >  2:592837/1‑71 (MQ=255)
                                                                                  aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                     >  2:401218/1‑78 (MQ=255)
                                                                                  aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                     <  1:401218/78‑1 (MQ=255)
                                                                                    gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                       >  3:929447/1‑90 (MQ=255)
                                                                                    gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                       <  4:929447/90‑1 (MQ=255)
                                                                                    gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTga                                                                                                               >  3:577113/1‑98 (MQ=255)
                                                                                    gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTga                                                                                                               <  4:577113/98‑1 (MQ=255)
                                                                                    gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaagaa                                                             <  3:1012206/148‑1 (MQ=255)
                                                                                     aaaGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGc                                                                                                  <  3:6314/110‑1 (MQ=255)
                                                                                     aaaGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGc                                                                                                  >  4:6314/1‑110 (MQ=255)
                                                                                          tAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTgg                                                        <  1:473993/147‑1 (MQ=255)
                                                                                          tAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTgg                                                        >  2:473993/1‑147 (MQ=255)
                                                                                            aCGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                              >  4:877607/1‑91 (MQ=255)
                                                                                            aCGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                              <  3:877607/91‑1 (MQ=255)
                                                                                               aTCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTg                                                                              <  1:1083617/120‑1 (MQ=255)
                                                                                               aTCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTg                                                                              >  2:1083617/1‑120 (MQ=255)
                                                                                               aTCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                    <  1:1106294/130‑1 (MQ=255)
                                                                                               aTCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                    >  2:1106294/1‑130 (MQ=255)
                                                                                                 cAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGc                                                                                        <  4:39311/108‑1 (MQ=255)
                                                                                                 cAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGc                                                                                        >  3:39311/1‑108 (MQ=255)
                                                                                                 cAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                      <  4:830042/126‑1 (MQ=255)
                                                                                                 cAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAAGATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                      >  3:830042/1‑126 (MQ=255)
                                                                                                  aGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                    >  4:19799/1‑127 (MQ=255)
                                                                                                  aGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                    <  3:19799/127‑1 (MQ=255)
                                                                                                       gCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCt                                                                                >  1:553888/1‑110 (MQ=255)
                                                                                                       gCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCt                                                                                <  2:553888/110‑1 (MQ=255)
                                                                                                       gCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTTAAGAACGTg                                                         >  1:772256/1‑133 (MQ=255)
                                                                                                       gCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTg                                                         <  2:772256/133‑1 (MQ=255)
                                                                                                         cAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCt                                                                                <  2:825248/108‑1 (MQ=255)
                                                                                                         cAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCt                                                                                >  1:825248/1‑108 (MQ=255)
                                                                                                          aaCCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAg                                                                                                                                          >  1:104839/1‑49 (MQ=255)
                                                                                                          aaCCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAg                                                                                                                                          <  2:104839/49‑1 (MQ=255)
                                                                                                           aCCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTtc                                      >  4:281245/1‑148 (MQ=255)
                                                                                                              tCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTtcatc                                   <  3:281245/148‑1 (MQ=255)
                                                                                                                          ttgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTAt                                                                                           <  4:543310/80‑1 (MQ=255)
                                                                                                                               gCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTg                                                         <  1:390934/109‑1 (MQ=255)
                                                                                                                               gCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTg                                                         >  2:390934/1‑109 (MQ=255)
                                                                                                                                cTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTga                                                                 <  2:923970/100‑1 (MQ=255)
                                                                                                                                cTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTga                                                                 >  1:923970/1‑100 (MQ=255)
                                                                                                                                  ggTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTg                        >  2:747669/1‑139 (MQ=255)
                                                                                                                                  ggTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTg                        <  1:747669/139‑1 (MQ=255)
                                                                                                                                  ggTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCt               >  1:36836/1‑148 (MQ=255)
                                                                                                                                         gCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaagaa                                                             <  1:823452/95‑1 (MQ=255)
                                                                                                                                         gCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaagaa                                                             >  2:823452/1‑95 (MQ=255)
                                                                                                                                          cTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTg            >  1:670742/1‑143 (MQ=255)
                                                                                                                                          cTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTg            <  2:670742/143‑1 (MQ=255)
                                                                                                                                          cTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGtt          <  1:113838/145‑1 (MQ=255)
                                                                                                                                          cTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGtt          >  2:113838/1‑145 (MQ=255)
                                                                                                                                           tCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                    >  1:662064/1‑86 (MQ=255)
                                                                                                                                           tCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                    <  2:662064/86‑1 (MQ=255)
                                                                                                                                               tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTg                                                         <  3:932789/93‑1 (MQ=255)
                                                                                                                                               tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCt               >  2:816476/1‑135 (MQ=255)
                                                                                                                                               tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCt               <  1:816476/135‑1 (MQ=255)
                                                                                                                                               tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGAc  >  3:178242/1‑148 (MQ=255)
                                                                                                                                               tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGAc  >  2:906428/1‑148 (MQ=255)
                                                                                                                                               tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGAc  >  1:96661/1‑148 (MQ=255)
                                                                                                                                                 aaGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                      <  4:201028/78‑1 (MQ=255)
                                                                                                                                                 aaGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                      >  3:201028/1‑78 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGAC  >  NC_000913/1507353‑1507643

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