Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 561,696:1 (G)6→7 100% coding (2000/2604 nt) sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913561,6901.G92.9% 273.4 / 13.7 85coding (1994/2604 nt)sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base G (34/44);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGC  >  NC_000913/561593‑561784
                                                                                                  |                                                                                               
ggCCACAACCTGAGTTATCGCGTAAGCCAGGGTGATACCTCGAATAATGGGTACAGGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCt                                                                                                  >  2:258337/1‑98 (MQ=255)
 gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggg                                                                                              >  1:855510/1‑100 (MQ=255)
 gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggg                                                                                              >  2:2309153/1‑100 (MQ=255)
  ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggg                                                                                             >  2:2238136/1‑100 (MQ=255)
  ccACAACCTGAGTGATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggg                                                                                             <  2:1603070/100‑1 (MQ=255)
   cacaACCTGAGTTATCACGTAAGCCAGGGTGGTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggggg                                                                                            <  2:1988008/100‑1 (MQ=255)
   cacaACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggggg                                                                                            <  2:2077845/100‑1 (MQ=255)
   cacaACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggggg                                                                                            >  1:1789098/1‑100 (MQ=255)
   cacaACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggggg                                                                                            <  2:2035224/100‑1 (MQ=255)
     caACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATGATCGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggggggg                                                                                         <  1:1122237/100‑2 (MQ=255)
       aCCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATTGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcg                                                                                       >  2:1203910/1‑97 (MQ=255)
       aCCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcg                                                                                       >  1:319032/1‑97 (MQ=255)
        ccTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcgg                                                                                      <  2:2092664/100‑5 (MQ=255)
         cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcggg                                                                                     <  1:730758/100‑6 (MQ=255)
         cTGAGTTATCACGTAAGCCAGGGGGATACGTCGAATAATGGGTACACGGGCAGCGCAACGGCAAGCTGGCAGGCCGCTTACGGTACGCTgggggggcgg                                                                                      >  1:1591131/1‑95 (MQ=255)
          tGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCCGGGGGGGCGGGt                                                                                     <  2:1701920/100‑1 (MQ=255)
           gAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGta                                                                                    <  2:1360706/100‑1 (MQ=255)
           gAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGta                                                                                    <  2:1511729/100‑1 (MQ=255)
            aGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtat                                                                                   <  2:798752/100‑1 (MQ=255)
             gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtata                                                                                  >  2:1569256/1‑100 (MQ=255)
              ttATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCCAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGGATaa                                                                                 >  1:1713475/1‑100 (MQ=255)
              ttATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATaa                                                                                 <  1:2428079/100‑1 (MQ=255)
              ttATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATaa                                                                                 <  1:179562/100‑1 (MQ=255)
               tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAAc                                                                                <  1:1834473/100‑1 (MQ=255)
                 tCACGTAAGCCAGGGTGATACCTCGAATAATGGGTAGACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTa                                                                              <  2:1816475/100‑1 (MQ=255)
                 tCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTa                                                                              >  2:923345/1‑100 (MQ=255)
                  cACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTAc                                                                             >  2:1077687/1‑100 (MQ=255)
                   aCGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACg                                                                            >  1:1776017/1‑100 (MQ=255)
                   aCGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACg                                                                            >  2:1691369/1‑100 (MQ=255)
                     gTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGAc                                                                          <  1:114399/100‑1 (MQ=255)
                      tAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGAcc                                                                         >  1:1866219/1‑100 (MQ=255)
                      tAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGAcc                                                                         >  1:1393347/1‑100 (MQ=255)
                        aGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGCCcg                                                                        >  1:2394801/1‑99 (MQ=255)
                         gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACcgcg                                                                      <  1:1694113/100‑1 (MQ=255)
                         gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACcgcg                                                                      <  2:1789107/100‑1 (MQ=255)
                           cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGAt                                                                    <  2:2397130/100‑1 (MQ=255)
                            aGGGTGATACCTCTAATAATGGGTACAGGGGGAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATc                                                                   <  1:480632/100‑1 (MQ=255)
                             gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCa                                                                  <  2:915576/100‑1 (MQ=255)
                             gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCa                                                                  <  1:1685705/100‑1 (MQ=255)
                              ggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCaa                                                                 <  1:1851657/100‑1 (MQ=255)
                                tGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACa                                                               <  2:2386732/100‑1 (MQ=255)
                                tGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACa                                                               <  1:1884715/100‑1 (MQ=255)
                                 gATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACAt                                                              <  2:1179330/100‑1 (MQ=255)
                                   tACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                           <  1:1114720/100‑1 (MQ=255)
                                    aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                           <  2:683348/100‑1 (MQ=255)
                                    aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                           <  1:2112644/100‑1 (MQ=255)
                                    aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACGGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                           <  1:1061471/100‑1 (MQ=255)
                                    aCCTCGAATAATGGGTACAAGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                           >  2:2377838/1‑100 (MQ=255)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGa                                                            >  2:1772940/1‑98 (MQ=255)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                          >  1:1895888/1‑100 (MQ=255)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                          >  2:788942/1‑100 (MQ=255)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                          <  2:1568746/100‑1 (MQ=255)
                                      cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                         >  2:1035353/1‑100 (MQ=255)
                                      cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                         >  2:2577827/1‑100 (MQ=255)
                                      cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                         >  2:2057525/1‑100 (MQ=255)
                                      cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                         >  1:1795296/1‑100 (MQ=255)
                                      cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                         >  1:2272455/1‑100 (MQ=255)
                                       tCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGtt                                                        >  2:2277392/1‑100 (MQ=255)
                                         gAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTaa                                                      >  1:1813572/1‑100 (MQ=255)
                                          aataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACt                                                    <  2:2059857/100‑1 (MQ=255)
                                           ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACt                                                    >  1:1152440/1‑100 (MQ=255)
                                           ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTg                                                   >  1:138556/1‑100 (MQ=255)
                                            taatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTg                                                   <  1:378538/100‑1 (MQ=255)
                                              atGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGc                                                 <  2:54717/100‑1 (MQ=255)
                                               tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                                <  2:572659/100‑1 (MQ=255)
                                                  gTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCt                                             <  1:2140586/100‑1 (MQ=255)
                                                  gTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCt                                             <  1:2139492/100‑1 (MQ=255)
                                                   tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                            >  1:538771/1‑100 (MQ=255)
                                                     cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTc                                          <  2:880980/100‑1 (MQ=255)
                                                      acGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGGACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCt                                         <  2:1020328/100‑1 (MQ=255)
                                                         ggCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTAGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGc                                      >  1:1414836/1‑100 (MQ=255)
                                                             gcgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtg                                  <  1:1099806/100‑1 (MQ=255)
                                                               gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                                >  2:1166762/1‑100 (MQ=255)
                                                                ccACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTgg                               >  2:2252146/1‑100 (MQ=255)
                                                                  aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACGGCGATCAACATGAGGTTAACTGGCTGctggctggggttgggg                               >  2:1677122/1‑90 (MQ=255)
                                                                  aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGGTGTCTGACGGTGTGGTc                             >  2:1121082/1‑100 (MQ=255)
                                                                  aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTgg                               >  1:1573046/1‑98 (MQ=255)
                                                                  aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTc                             <  2:1248734/100‑1 (MQ=255)
                                                                   cGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCg                            >  2:1638570/1‑100 (MQ=255)
                                                                   cGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCg                            >  2:1259264/1‑100 (MQ=255)
                                                                     gCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                          <  2:215969/100‑1 (MQ=255)
                                                                       aaaCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCa                        <  1:2094349/100‑1 (MQ=255)
                                                                         aCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATg                      >  2:1081590/1‑100 (MQ=255)
                                                                           tGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaa                    <  2:1859690/100‑1 (MQ=255)
                                                                            ggCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATAACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaaaa                  >  1:1654992/1‑100 (MQ=255)
                                                                                ggCCGCTAACTGTACGCTGGGGGAGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATg                <  1:1361554/99‑1 (MQ=255)
                                                                                 ggcGCTAACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg               <  1:658122/97‑1 (MQ=255)
                                                                                 gCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGCCCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATg                >  1:1355486/1‑98 (MQ=255)
                                                                                 gCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGc              <  2:2288334/100‑1 (MQ=255)
                                                                                  ccGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa             <  2:2032538/100‑1 (MQ=255)
                                                                                   cGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCAt            <  2:561770/100‑1 (MQ=255)
                                                                                      ttACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAAc         <  1:1166758/100‑1 (MQ=255)
                                                                                      ttACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAAc         <  1:2309139/100‑1 (MQ=255)
                                                                                      tctcGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAAc         <  1:1878375/97‑1 (MQ=255)
                                                                                         cGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAAcaga      >  2:1878385/1‑97 (MQ=255)
                                                                                         cGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCt      >  2:1296263/1‑100 (MQ=255)
                                                                                           gTACGCTGGGGGGCGGGTATAACTACGACCGCGATGAAGATGACGTTAACTGGTAGCTGTCTGGAGGTGTGGTCGGGCATGAAAATGGCATAACGCTGag   <  1:650317/100‑1 (MQ=255)
                                                                                            tACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGag   >  1:2534259/1‑100 (MQ=255)
                                                                                             aCGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagc  <  2:855513/100‑1 (MQ=255)
                                                                                             aCGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGagc  <  2:2107137/100‑1 (MQ=255)
                                                                                                  |                                                                                               
GGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGC  >  NC_000913/561593‑561784

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