Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 515,263 C→T P41S (CCT→TCT)  ybbL → putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913515,2630CT99.8% 1763.0 / ‑6.2 502P41S (CCT→TCT) ybbLputative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (209/291);  total (210/292)
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.11e-01

Only 100 of 515 total aligned reads displayed.

CAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGAT  >  NC_000913/515173‑515352
                                                                                          |                                                                                         
cAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCt                                                                                         >  2:902071/1‑93 (MQ=255)
 aaaaCGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCtt                                                                                        <  1:4596031/93‑1 (MQ=255)
 aaaaCGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCtt                                                                                        <  2:5306517/93‑1 (MQ=255)
  aaaCGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGttcttc                                                                                       <  1:2326465/93‑1 (MQ=255)
   aaCGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTtcttct                                                                                      <  1:1009830/93‑1 (MQ=255)
       tAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGtt                                                                                  >  1:19200/1‑93 (MQ=255)
         ggATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTtgt                                                                                >  1:9582835/1‑93 (MQ=255)
         ggATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTtgt                                                                                >  2:1076866/1‑93 (MQ=255)
          gATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTtgtg                                                                               <  1:6716034/93‑1 (MQ=255)
           atatCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTgg                                                                              <  1:6391491/93‑1 (MQ=255)
             atCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCa                                                                            <  2:9719978/93‑1 (MQ=255)
                tGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCaaaa                                                                         <  1:2226977/93‑1 (MQ=255)
                tGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCaaaa                                                                         >  2:1618221/1‑93 (MQ=255)
                 ggCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAg                                                                        <  2:9522254/93‑1 (MQ=255)
                   cGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTa                                                                      <  1:7183971/93‑1 (MQ=255)
                   cGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTa                                                                      >  2:9231469/1‑93 (MQ=255)
                      gTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACgc                                                                   >  2:2533219/1‑93 (MQ=255)
                      gTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTATGGTTGTGGCAAAAGTACgc                                                                   <  1:9452844/93‑1 (MQ=255)
                       tGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACgct                                                                  <  1:6950174/93‑1 (MQ=255)
                        gATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACgctg                                                                 <  2:8643142/93‑1 (MQ=255)
                          tGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACgctgct                                                               >  2:3563210/1‑93 (MQ=255)
                          tGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACgctgct                                                               <  1:2409129/93‑1 (MQ=255)
                           gCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTa                                                              <  2:3730289/93‑1 (MQ=255)
                                gATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAAt                                                         >  1:3411963/1‑93 (MQ=255)
                                 aTTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATa                                                        >  1:5600210/1‑93 (MQ=255)
                                  ttCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAg                                                       >  2:4483510/1‑93 (MQ=255)
                                  ttCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAg                                                       <  2:6644398/93‑1 (MQ=255)
                                  ttCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAg                                                       >  2:7728429/1‑93 (MQ=255)
                                   tCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGt                                                      >  1:9201410/1‑93 (MQ=255)
                                    cTTAATAACATCAATTTGTCGCTGCGTGCTAGCGAATTTAGGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGtt                                                     <  2:1443364/93‑1 (MQ=255)
                                     ttAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTg                                                    >  2:837143/1‑93 (MQ=255)
                                     ttAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTg                                                    <  2:7605243/93‑1 (MQ=255)
                                     ttAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTg                                                    <  2:3316626/93‑1 (MQ=255)
                                     ttAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTg                                                    >  2:6483626/1‑93 (MQ=255)
                                      taataaCATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGc                                                   <  1:4863139/93‑1 (MQ=255)
                                      taataaCATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGc                                                   >  1:726516/1‑93 (MQ=255)
                                       aataaCATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCt                                                  >  1:3324865/1‑93 (MQ=255)
                                       aataaCATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGCTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCt                                                  >  2:8502068/1‑93 (MQ=255)
                                        ataaCATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCtt                                                 <  1:1494690/93‑1 (MQ=255)
                                         taaCATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTc                                                <  1:1649193/93‑1 (MQ=255)
                                           aCATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCAt                                              <  1:5706714/93‑1 (MQ=255)
                                           aCATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCAt                                              <  2:8309349/93‑1 (MQ=255)
                                           aCATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCAt                                              <  1:212986/93‑1 (MQ=255)
                                            cATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATGACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCAtt                                             <  2:4960613/93‑1 (MQ=255)
                                             aTCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTg                                            <  2:1836452/93‑1 (MQ=255)
                                               cAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGAt                                          >  2:7495333/1‑93 (MQ=255)
                                               cAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGAt                                          >  1:14388/1‑93 (MQ=255)
                                                aaTTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATc                                         <  2:1745213/93‑1 (MQ=255)
                                                 aTTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCa                                        <  2:6793024/93‑1 (MQ=255)
                                                   ttttCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGc                                      <  1:5266309/93‑1 (MQ=255)
                                                    tttCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGcc                                     <  2:8097567/93‑1 (MQ=255)
                                                     ttCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGccc                                    <  1:8284259/93‑1 (MQ=255)
                                                        gCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAAc                                 <  2:9843989/93‑1 (MQ=255)
                                                        gCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAAc                                 <  2:6193233/93‑1 (MQ=255)
                                                        gCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAAc                                 <  2:2813673/93‑1 (MQ=255)
                                                        gCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAAc                                 <  2:2187055/93‑1 (MQ=255)
                                                        gCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAAc                                 >  1:6147713/1‑93 (MQ=255)
                                                         cTGCGTGCTGGCGAATTTAAGTTCATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAAcc                                <  2:5507882/93‑1 (MQ=255)
                                                         cTGCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAAcc                                <  2:5913090/93‑1 (MQ=255)
                                                           gCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAg                              >  1:4439291/1‑93 (MQ=255)
                                                           gCGTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAg                              <  1:1728822/93‑1 (MQ=255)
                                                             gTGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCg                            >  1:139243/1‑93 (MQ=255)
                                                              tGCTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCgg                           >  2:8794724/1‑93 (MQ=255)
                                                                cTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGaa                         <  1:7951407/93‑1 (MQ=255)
                                                                cTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGaa                         <  2:9428509/93‑1 (MQ=255)
                                                                cTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGaa                         <  1:2828299/93‑1 (MQ=255)
                                                                cTGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGaa                         <  1:300975/93‑1 (MQ=255)
                                                                 tGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAAc                        >  1:5386337/1‑93 (MQ=255)
                                                                 tGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAAc                        >  1:7732947/1‑93 (MQ=255)
                                                                 tGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAAc                        >  2:6249863/1‑93 (MQ=255)
                                                                 tGGCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAAc                        >  2:4633995/1‑93 (MQ=255)
                                                                  ggCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACg                       >  1:8068670/1‑93 (MQ=255)
                                                                  ggCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACg                       <  2:716256/93‑1 (MQ=255)
                                                                  ggCGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACg                       >  1:5072740/1‑93 (MQ=255)
                                                                    cGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGtt                     >  2:3378882/1‑93 (MQ=255)
                                                                    cGAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGtt                     >  1:7875838/1‑93 (MQ=255)
                                                                     gAATTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTa                    >  1:8492756/1‑93 (MQ=255)
                                                                      aaTTTAAGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGTCCAACCAGCGGAACGTTAc                   >  1:3068062/1‑93 (MQ=255)
                                                                            aGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTg             <  2:4564362/93‑1 (MQ=255)
                                                                            aGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTg             <  2:6873212/93‑1 (MQ=255)
                                                                            aGTTAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTg             <  1:2472232/93‑1 (MQ=255)
                                                                             gTTAATTACCGGATCTTATGGTTGTGGCAAAAGTACGCTGCTAAAAATAGATGCTTCATTGATCAGCACAACCCGCGGAACGTTACTGGTTGa            >  2:5389030/1‑93 (MQ=255)
                                                                              gtAATTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGaa           <  1:1143586/92‑1 (MQ=255)
                                                                                aaTTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAgg         >  1:8783518/1‑93 (MQ=255)
                                                                                 aTTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGt        <  2:5437972/93‑1 (MQ=255)
                                                                                 aTTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGt        <  1:3702502/93‑1 (MQ=255)
                                                                                 aTTACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGt        <  1:5772645/93‑1 (MQ=255)
                                                                                  ttACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTg       <  1:9689723/93‑1 (MQ=255)
                                                                                   tACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGa      >  2:3232630/1‑93 (MQ=255)
                                                                                   tACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGa      <  2:2293293/93‑1 (MQ=255)
                                                                                   tACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGa      <  1:6618641/93‑1 (MQ=255)
                                                                                   tACCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGa      <  1:7457297/93‑1 (MQ=255)
                                                                                    aCCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAg     >  1:3563998/1‑93 (MQ=255)
                                                                                    aCCGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTAATCAGCCCAACCACCGGAAGGTTACTGTTGGAAGGTGAg     >  1:4341618/1‑93 (MQ=255)
                                                                                     ccGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAgg    <  1:4102147/93‑1 (MQ=255)
                                                                                     ccGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAgg    <  2:4343230/93‑1 (MQ=255)
                                                                                      cGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGa   <  1:9935440/93‑1 (MQ=255)
                                                                                      cGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGa   <  1:3816332/93‑1 (MQ=255)
                                                                                      cGGTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGa   <  2:2008986/93‑1 (MQ=255)
                                                                                       ggTTCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGTTCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGAt  <  2:1244743/93‑1 (MQ=255)
                                                                                          |                                                                                         
CAAAACGTAGGATATCTGGCGGGTGATGCGAAGATTCTTAATAACATCAATTTTTCGCTGCGTGCTGGCGAATTTAAGTTAATTACCGGTCCTTCTGGTTGTGGCAAAAGTACGCTGCTAAAAATAGTTGCTTCATTGATCAGCCCAACCAGCGGAACGTTACTGTTTGAAGGTGAGGAT  >  NC_000913/515173‑515352

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑