Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 559,101 A→G N61S (AAT→AGT)  sfmD → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913559,1010AG85.3% 643.2 / 99.8 272N61S (AAT→AGT) sfmDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (15/25);  new base G (126/106);  total (141/131)
Fisher's exact test for biased strand distribution p-value = 5.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-02

Only 100 of 284 total aligned reads displayed.

CTTTTTAAGCTATGCGGAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAATCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGAATCGACAACAGAA  >  NC_000913/559009‑559192
                                                                                            |                                                                                           
ctTTTTAAGCTATGCGGAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTaa                                                                                             <  2:5948892/93‑1 (MQ=255)
 tttttAAGCTATGCGGAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAat                                                                                            >  2:955769/1‑93 (MQ=255)
   tttAAGCTATGCGGAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAgtca                                                                                          <  2:5083075/93‑1 (MQ=255)
    ttAAGCTATGCGGAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAatcat                                                                                         >  2:2392792/1‑93 (MQ=255)
     tAAGCTATGCGGAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAgtcatc                                                                                        >  1:6168700/1‑93 (MQ=255)
      aaGCTATGCGGAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGtcatca                                                                                       >  1:6994414/1‑93 (MQ=255)
         cTATGCGGAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAAcc                                                                                    <  2:4317818/93‑1 (MQ=255)
         cTATGCGGAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAAcc                                                                                    <  1:5818297/93‑1 (MQ=255)
           aTGCGGAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAATCATCAACCTg                                                                                  <  2:5626178/93‑1 (MQ=255)
               ggAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCggg                                                                              >  1:2774301/1‑93 (MQ=255)
                gAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGc                                                                             >  2:1537148/1‑93 (MQ=255)
                gAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGc                                                                             >  2:2442210/1‑93 (MQ=255)
                gAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGc                                                                             >  2:1681222/1‑93 (MQ=255)
                gAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGc                                                                             >  1:1113090/1‑93 (MQ=255)
                gAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGc                                                                             >  1:644994/1‑93 (MQ=255)
                gAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGc                                                                             <  1:2545191/93‑1 (MQ=255)
                gAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGc                                                                             >  1:3638475/1‑93 (MQ=255)
                 aaaGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCg                                                                            >  1:3073164/1‑93 (MQ=255)
                 aaaGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCg                                                                            >  2:4345529/1‑93 (MQ=255)
                  aaGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCgt                                                                           >  2:7422556/1‑93 (MQ=255)
                   aGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCgtg                                                                          >  2:807662/1‑93 (MQ=255)
                   aGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCggg                                                                          >  2:5429993/1‑91 (MQ=255)
                    gTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCgtgt                                                                         <  2:2847276/93‑1 (MQ=255)
                    gTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAATCATCAACCTGCGGGCgtgt                                                                         >  1:3717246/1‑93 (MQ=255)
                      tATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTCCTGATTTATCGCGCTTTGAGAGAGGTAATCATCAACCTGCGGGCGTGTAt                                                                       <  1:7371506/93‑1 (MQ=255)
                       aTTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATc                                                                      <  2:3028692/93‑1 (MQ=255)
                       aTTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATc                                                                      <  2:4922211/93‑1 (MQ=255)
                        tttCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCg                                                                     <  1:7493602/93‑1 (MQ=255)
                        tttCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCg                                                                     >  2:3238740/1‑93 (MQ=255)
                        tttCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCg                                                                     <  1:1427519/93‑1 (MQ=255)
                         ttCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCgg                                                                    <  2:6845398/93‑1 (MQ=255)
                         ttCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCgg                                                                    <  1:484014/93‑1 (MQ=255)
                           cAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGt                                                                  <  1:5426098/93‑1 (MQ=255)
                           cAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGt                                                                  <  1:4224687/93‑1 (MQ=255)
                            aaCCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTg                                                                 >  1:4174094/1‑93 (MQ=255)
                            aaCCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTg                                                                 <  2:7471282/93‑1 (MQ=255)
                            aaCCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTg                                                                 <  2:5325568/93‑1 (MQ=255)
                            aaCCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTg                                                                 <  1:1509701/93‑1 (MQ=255)
                            aaCCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAATCATCAACCTGCGGGCGTGTATCGGGTg                                                                 >  2:4810763/1‑93 (MQ=255)
                             aCCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTgg                                                                >  2:4060228/1‑93 (MQ=255)
                             aCCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTgg                                                                >  2:3971954/1‑93 (MQ=255)
                             aCCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAATCATCAACCTGCGGGCGTGTATCGGGTgg                                                                >  1:1682278/1‑93 (MQ=255)
                              cccGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGa                                                               >  1:631094/1‑93 (MQ=255)
                              cccGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGa                                                               >  1:2881105/1‑93 (MQ=255)
                               ccGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGAt                                                              >  1:3359969/1‑93 (MQ=255)
                                cGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGAtc                                                             >  2:3926192/1‑93 (MQ=255)
                                 ggCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGAtct                                                            >  2:1619241/1‑93 (MQ=255)
                                 ggCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATAAACCTGCGGGCGTGTATCGGGTGGAtct                                                            >  2:715670/1‑93 (MQ=255)
                                   cATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATctct                                                          >  1:687235/1‑93 (MQ=255)
                                    aTTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTg                                                         >  1:6599890/1‑93 (MQ=255)
                                    aTTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTg                                                         >  1:5037790/1‑93 (MQ=255)
                                    aTTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAATCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTg                                                         <  2:7075120/93‑1 (MQ=255)
                                       tttATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCg                                                      <  2:5683971/93‑1 (MQ=255)
                                           ttaGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAat                                                  >  2:2304964/1‑93 (MQ=255)
                                           ttaGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTATATCGGGTGGATCTCTGGCGTAat                                                  <  1:474698/93‑1 (MQ=255)
                                             agagAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAATCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAatga                                                <  1:4241086/93‑1 (MQ=255)
                                              gagaATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAatgat                                               <  1:6779519/93‑1 (MQ=255)
                                                gaATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAAtgatga                                             >  2:3451970/1‑93 (MQ=255)
                                                gaATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAAtgatga                                             <  1:2324723/93‑1 (MQ=255)
                                                 aaTGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAATCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAg                                            <  2:2074679/93‑1 (MQ=255)
                                                  aTGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAATCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGt                                           <  2:4715294/93‑1 (MQ=255)
                                                   tGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGtt                                          >  2:1912558/1‑93 (MQ=255)
                                                   tGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGtt                                          <  1:546703/93‑1 (MQ=255)
                                                     gCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCa                                        <  2:1010659/93‑1 (MQ=255)
                                                      cacaTCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCAt                                       >  1:3853335/1‑93 (MQ=255)
                                                      cacaTCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCAt                                       <  1:2078693/93‑1 (MQ=255)
                                                        caTCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTg                                     >  1:1265570/1‑93 (MQ=255)
                                                        caTCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTg                                     <  2:2618602/93‑1 (MQ=255)
                                                           ccGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGtt                                  >  1:1866905/1‑93 (MQ=255)
                                                            cGTTGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTc                                 >  2:121993/1‑93 (MQ=255)
                                                              ttGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGc                               >  1:7178068/1‑93 (MQ=255)
                                                              ttGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGc                               >  2:7294973/1‑93 (MQ=255)
                                                              ttGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTATTGGTTCGCa                              >  1:4732273/1‑93 (MQ=255)
                                                               tGCTGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCa                              <  1:4023311/93‑1 (MQ=255)
                                                                  tGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGa                           >  2:2503489/1‑93 (MQ=255)
                                                                  tGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGa                           <  1:5741676/93‑1 (MQ=255)
                                                                  tGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGa                           <  2:6759694/93‑1 (MQ=255)
                                                                  tGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGa                           <  1:2266842/93‑1 (MQ=255)
                                                                  tGATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGa                           <  2:107122/93‑1 (MQ=255)
                                                                   gATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGat                          <  2:3690986/93‑1 (MQ=255)
                                                                   gATTTATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGat                          >  2:6326024/1‑93 (MQ=255)
                                                                    aTTTATCGCGCTTTGAGAGAGGTAATCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGata                         <  1:393679/93‑1 (MQ=255)
                                                                     tttATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGatat                        >  2:4631955/1‑93 (MQ=255)
                                                                       tATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCg                      <  2:3882652/93‑1 (MQ=255)
                                                                       tATCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCg                      >  2:2335932/1‑93 (MQ=255)
                                                                         tCGCGCTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTa                    >  2:5876071/1‑93 (MQ=255)
                                                                            cgcTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTAttt                 <  1:5777201/93‑1 (MQ=255)
                                                                             gcTTTGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTg                >  2:6626407/1‑93 (MQ=255)
                                                                               tttGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGaa              <  1:6051999/93‑1 (MQ=255)
                                                                               tttGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGaa              >  1:803642/1‑93 (MQ=255)
                                                                               tttGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGaa              <  1:2025530/93‑1 (MQ=255)
                                                                               tttGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGaa              <  2:759938/93‑1 (MQ=255)
                                                                                ttGAGAGAGGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGAAt             <  1:5532408/93‑1 (MQ=255)
                                                                                    gagagGTAGTCATCAACCTGCGGGCGTGTATCGGGTTGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGAATCGac         <  2:6252088/93‑1 (MQ=255)
                                                                                     agagGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGAATCGaca        <  2:4411858/93‑1 (MQ=255)
                                                                                      gagGTAATCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGAATCGacaa       <  2:2230973/93‑1 (MQ=255)
                                                                                       agGTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGAATCGacaac      <  2:4512329/93‑1 (MQ=255)
                                                                                        ggTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGAATCGacaaca     >  1:437562/1‑93 (MQ=255)
                                                                                        ggTAGTCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGAATCGacaaca     <  1:3745680/93‑1 (MQ=255)
                                                                                         gTAATCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGAATCGACAACAg    >  2:3852874/1‑93 (MQ=255)
                                                                                            |                                                                                           
CTTTTTAAGCTATGCGGAAAGTTATTTCAACCCGGCATTTTTATTAGAGAATGGCACATCCGTTGCTGATTTATCGCGCTTTGAGAGAGGTAATCATCAACCTGCGGGCGTGTATCGGGTGGATCTCTGGCGTAATGATGAGTTCATTGGTTCGCAGGATATCGTATTTGAATCGACAACAGAA  >  NC_000913/559009‑559192

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