Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,507,501 A→C 7.8% I241L (ATT→CTT)  ydcP → putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,507,5010AC 7.8% 362.0 / 20.4 102I241L (ATT→CTT) ydcPputative peptidase
Reads supporting (aligned to +/- strand):  ref base A (47/47);  new base C (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 = 6.52e-01

Only 100 of 106 total aligned reads displayed.

CGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACG  >  NC_000913/1507358‑1507644
                                                                                                                                               |                                                                                                                                               
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                                                                                                                                             >  4:391055/1‑126 (MQ=255)
                      cacaTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                             <  3:391055/126‑1 (MQ=255)
                          ttGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAgg                                                                                                                                          <  1:904823/125‑1 (MQ=255)
                          ttGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAgg                                                                                                                                          >  2:904823/1‑125 (MQ=255)
                             aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                    <  3:387115/127‑7 (MQ=255)
                             aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                    <  1:896859/127‑7 (MQ=255)
                             aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                    >  2:896859/1‑121 (MQ=255)
                             aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCCCTTCAATCTTGAAggagcgta                                                                                                                                    >  4:387115/1‑121 (MQ=255)
                              aaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                        <  1:183985/139‑1 (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                                                                                                                             <  1:979422/110‑1 (MQ=255)
                                                      ccTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGat                                                                                                                             >  2:979422/1‑110 (MQ=255)
                                                        tATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGa                                                                                     >  4:161892/1‑148 (MQ=255)
                                                         aTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAggg                                                                                                                                         <  4:788030/95‑1 (MQ=255)
                                                         aTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAggg                                                                                                                                         >  3:788030/1‑95 (MQ=255)
                                                         aTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaagaa                                                                                                            <  1:285331/124‑1 (MQ=255)
                                                         aTGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaagaa                                                                                                            >  2:285331/1‑124 (MQ=255)
                                                          tGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                               <  1:165876/120‑1 (MQ=255)
                                                          tGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                               >  2:165876/1‑120 (MQ=255)
                                                            aaaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAATCTTGAAgg                                                                                                                                          <  3:87318/91‑1 (MQ=255)
                                                            aaaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAATCTTGAAgg                                                                                                                                          >  4:87318/1‑91 (MQ=255)
                                                             aaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGc                                                                                         >  4:573483/1‑139 (MQ=255)
                                                               aaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGa                                                                              <  3:161892/148‑1 (MQ=255)
                                                               aaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGa                                                                              >  4:1012206/1‑148 (MQ=255)
                                                                   tCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGa                                                                                                                         >  4:885824/1‑101 (MQ=255)
                                                                   tCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGa                                                                                                                         <  3:885824/101‑1 (MQ=255)
                                                                      gctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                      >  4:11041/1‑85 (MQ=255)
                                                                      gctgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                      <  3:11041/85‑1 (MQ=255)
                                                                         gTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaaga                                                                                                             <  2:981050/107‑1 (MQ=255)
                                                                         gTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaaga                                                                                                             >  1:981050/1‑107 (MQ=255)
                                                                             aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                             >  2:592837/1‑71 (MQ=255)
                                                                             aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGa                                                                                                                                             <  1:592837/71‑1 (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                                                         >  2:473993/1‑147 (MQ=255)
                                                                                     tAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTgg                                                         <  1:473993/147‑1 (MQ=255)
                                                                                       aCGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                               <  3:877607/91‑1 (MQ=255)
                                                                                       aCGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTgaa                                                                                                               >  4:877607/1‑91 (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                                                                     <  3:19799/127‑1 (MQ=255)
                                                                                             aGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                     >  4:19799/1‑127 (MQ=255)
                                                                                                  gCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCt                                                                                 >  1:553888/1‑110 (MQ=255)
                                                                                                  gCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCt                                                                                 <  2:553888/110‑1 (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)
                                                                                                         tCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTtcatc                                    <  3:281245/148‑1 (MQ=255)
                                                                                                                     ttgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTAt                                                                                            <  4:543310/80‑1 (MQ=255)
                                                                                                                     ttgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTAt                                                                                            >  3:543310/1‑80 (MQ=255)
                                                                                                                          gCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTg                                                          >  2:390934/1‑109 (MQ=255)
                                                                                                                          gCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTg                                                          <  1:390934/109‑1 (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                                                              >  2:823452/1‑95 (MQ=255)
                                                                                                                                    gCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaagaa                                                              <  1:823452/95‑1 (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                                                                     <  2:662064/86‑1 (MQ=255)
                                                                                                                                      tCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattat                                                                     >  1:662064/1‑86 (MQ=255)
                                                                                                                                          tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTg                                                          <  3:932789/93‑1 (MQ=255)
                                                                                                                                          tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCt                <  1:816476/135‑1 (MQ=255)
                                                                                                                                          tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCt                >  2:816476/1‑135 (MQ=255)
                                                                                                                                          tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGAc   >  1:96661/1‑148 (MQ=255)
                                                                                                                                          tCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGAc   >  3:178242/1‑148 (MQ=255)
                                                                                                                                            aaGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                       <  4:201028/78‑1 (MQ=255)
                                                                                                                                            aaGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCatt                                                                       >  3:201028/1‑78 (MQ=255)
                                                                                                                                               aTTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTAt                                                                                            <  4:272258/54‑1 (MQ=255)
                                                                                                                                               aTTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTAt                                                                                            >  3:272258/1‑54 (MQ=255)
                                                                                                                                               aTTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACg  <  3:903265/144‑1 (MQ=255)
                                                                                                                                               aTTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACg  >  4:903265/1‑144 (MQ=255)
                                                                                                                                               |                                                                                                                                               
CGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTGAAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACG  >  NC_000913/1507358‑1507644

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