Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 477,192 A→G L219P (CTG→CCG)  ylaB ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913477,1920AG99.2% 864.3 / ‑2.9 265L219P (CTG→CCG) ylaBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (142/121);  total (142/123)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-01

Only 100 of 268 total aligned reads displayed.

AACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCGCT  >  NC_000913/477102‑477280
                                                                                          |                                                                                        
 aCCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGa                                                                                       <  2:375295/93‑1 (MQ=255)
 aCCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGa                                                                                       <  1:5198884/93‑1 (MQ=255)
  ccAAATAAAGTCTTGCCGATGCCAGCCTTTCTGCATCATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGat                                                                                      <  1:1290443/93‑1 (MQ=255)
  ccAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGat                                                                                      <  1:3223595/93‑1 (MQ=255)
  ccAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCCAAATATTCATCTCCGGTAAGGGCGGGat                                                                                      >  1:4082914/1‑93 (MQ=255)
   cAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGata                                                                                     <  1:1769712/93‑1 (MQ=255)
   cAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGata                                                                                     <  2:1395558/93‑1 (MQ=255)
   cAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGata                                                                                     >  2:2269072/1‑93 (MQ=255)
   cAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGata                                                                                     <  1:2526407/93‑1 (MQ=255)
     aaTAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATc                                                                                   <  1:5677701/93‑1 (MQ=255)
     aaTAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATc                                                                                   <  1:389366/93‑1 (MQ=255)
     aaTAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATc                                                                                   >  1:2217834/1‑93 (MQ=255)
     aaTAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATc                                                                                   <  2:5023774/93‑1 (MQ=255)
          aGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTaaa                                                                              >  1:840411/1‑93 (MQ=255)
          aGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTaaa                                                                              <  2:2935760/93‑1 (MQ=255)
             cTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGa                                                                           >  1:81408/1‑93 (MQ=255)
              ttGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGAt                                                                          <  1:2973888/93‑1 (MQ=255)
               tGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGGTATCGTAAATGAtt                                                                         >  1:6245487/1‑93 (MQ=255)
                 ccGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTcc                                                                       <  1:6391769/93‑1 (MQ=255)
                 ccGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTcc                                                                       <  2:529778/93‑1 (MQ=255)
                  cGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCa                                                                      >  1:4010352/1‑93 (MQ=255)
                  cGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCa                                                                      <  2:5400380/93‑1 (MQ=255)
                   gATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCat                                                                     <  2:345291/93‑1 (MQ=255)
                    aTGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCatt                                                                    >  2:3502796/1‑93 (MQ=255)
                     tGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCatta                                                                   <  2:1225141/93‑1 (MQ=255)
                     tGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCatta                                                                   <  1:2678914/93‑1 (MQ=255)
                     tGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCatta                                                                   <  1:470707/93‑1 (MQ=255)
                      gCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCattat                                                                  >  2:4026976/1‑93 (MQ=255)
                       ccAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCattatt                                                                 >  2:1883175/1‑93 (MQ=255)
                       ccAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCattatt                                                                 <  2:193157/93‑1 (MQ=255)
                       ccAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCattatt                                                                 <  2:474190/93‑1 (MQ=255)
                         aGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTAtttt                                                               <  2:1836256/93‑1 (MQ=255)
                         aGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTAtttt                                                               >  2:3715016/1‑93 (MQ=255)
                          gCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTc                                                              >  1:3633556/1‑93 (MQ=255)
                           ccTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCg                                                             >  1:3348994/1‑93 (MQ=255)
                            cTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGa                                                            <  1:5378676/93‑1 (MQ=255)
                            cTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGa                                                            >  2:5916532/1‑93 (MQ=255)
                             tttCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGat                                                           >  2:2224249/1‑93 (MQ=255)
                             tttCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGat                                                           >  2:1105891/1‑93 (MQ=255)
                              ttCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGata                                                          >  2:5891997/1‑93 (MQ=255)
                               tCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGatat                                                         <  1:5407202/93‑1 (MQ=255)
                               tCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGatat                                                         <  1:2569006/93‑1 (MQ=255)
                                 tGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGc                                                       >  2:5165350/1‑93 (MQ=255)
                                  gCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCt                                                      >  2:5988142/1‑93 (MQ=255)
                                  gCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCt                                                      >  1:6113704/1‑93 (MQ=255)
                                    aaCATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCa                                                    >  2:6280328/1‑93 (MQ=255)
                                    aaCATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCa                                                    >  2:3116154/1‑93 (MQ=255)
                                     aCATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCAc                                                   >  1:5518433/1‑93 (MQ=255)
                                        ttttCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCgg                                                >  1:1011299/1‑93 (MQ=255)
                                        ttttCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCgg                                                <  2:5645558/93‑1 (MQ=255)
                                          ttCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGgt                                              <  1:3519713/93‑1 (MQ=255)
                                           tCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGgtg                                             <  2:848899/93‑1 (MQ=255)
                                             gTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGtt                                           >  1:4507171/1‑93 (MQ=255)
                                             gTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGtt                                           >  1:6000584/1‑93 (MQ=255)
                                              ttGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGttt                                          >  1:1838174/1‑93 (MQ=255)
                                               tGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGtttt                                         <  1:3675666/93‑1 (MQ=255)
                                               tGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGtttt                                         <  2:1267692/93‑1 (MQ=255)
                                                gAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGttttt                                        >  1:5648697/1‑93 (MQ=255)
                                                gAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGttttt                                        >  1:4843591/1‑93 (MQ=255)
                                                gAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGttttt                                        >  2:2437693/1‑93 (MQ=255)
                                                  aGCCCATGTGATGAGCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTg                                      <  1:3807614/93‑1 (MQ=255)
                                                    cccATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTa                                    <  1:2817845/93‑1 (MQ=255)
                                                    cccATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTa                                    <  2:6197392/93‑1 (MQ=255)
                                                     ccATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAg                                   >  2:5518344/1‑93 (MQ=255)
                                                     ccATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAg                                   <  2:2106863/93‑1 (MQ=255)
                                                      cATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGc                                  >  1:5910710/1‑93 (MQ=255)
                                                      cATGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGc                                  >  1:3747660/1‑93 (MQ=255)
                                                       aTGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGcc                                 <  1:2787494/93‑1 (MQ=255)
                                                       aTGTGATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGcc                                 <  2:4947496/93‑1 (MQ=255)
                                                        tgtgATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCt                                >  2:330713/1‑93 (MQ=255)
                                                        tgtgATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCt                                >  2:4870971/1‑93 (MQ=255)
                                                         gtgATTATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTg                               <  1:4613231/93‑1 (MQ=255)
                                                         gtgATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTg                               >  2:2723341/1‑93 (MQ=255)
                                                         gtgATGATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTg                               <  1:1212733/93‑1 (MQ=255)
                                                         gtgATAATCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTg                               >  1:1043593/1‑93 (MQ=255)
                                                          tgatgaTCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTgg                              <  2:6644888/93‑1 (MQ=255)
                                                           gatgatCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGt                             >  1:209869/1‑93 (MQ=255)
                                                            atgatCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGta                            <  2:381741/93‑1 (MQ=255)
                                                             tgatCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGtaa                           >  2:1631018/1‑93 (MQ=255)
                                                             tgatCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGtaa                           <  1:6566604/93‑1 (MQ=255)
                                                              gatCGAAATATTCATCTCCGGTAAGGTCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGtaat                          <  1:1429827/93‑1 (MQ=255)
                                                              gatCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGtaat                          <  2:4418400/93‑1 (MQ=255)
                                                              gatCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGtaat                          <  1:4879245/93‑1 (MQ=255)
                                                               atCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGtaata                         <  1:5052307/93‑1 (MQ=255)
                                                               atCGAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGtaata                         >  1:4301518/1‑93 (MQ=255)
                                                                  gAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAAtt                      <  1:2071370/93‑1 (MQ=255)
                                                                  gAAATATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAAtt                      >  1:5852123/1‑93 (MQ=255)
                                                                   aaaTATTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTc                     >  2:3154072/1‑93 (MQ=255)
                                                                     atatTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTccc                   >  2:6464060/1‑93 (MQ=255)
                                                                      tatTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGCTTTTTGTAGCCTGGTAATAATTCCCt                  >  2:2584362/1‑93 (MQ=255)
                                                                       atTCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTg                 <  1:1570897/93‑1 (MQ=255)
                                                                        ttCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGa                >  2:1458344/1‑93 (MQ=255)
                                                                         tCATCTCCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAg               >  2:5688984/1‑93 (MQ=255)
                                                                              tcCGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAAtt          >  2:4632465/1‑93 (MQ=255)
                                                                                cGGTAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTc        <  2:90273/93‑1 (MQ=255)
                                                                                  gTAAGGGCGGGATATCGTAAATGATTCCCTTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTTAt      >  1:2341816/1‑93 (MQ=255)
                                                                                   tAAGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATc     >  2:6037140/1‑93 (MQ=255)
                                                                                     aGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCgc   >  1:4157281/1‑93 (MQ=255)
                                                                                     aGGGCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCgc   >  2:781917/1‑93 (MQ=255)
                                                                                      gggCGGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCgct  >  1:783692/1‑93 (MQ=255)
                                                                                          |                                                                                        
AACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGTAAGGGCAGGATATCGTAAATGATTCCATTATTTTCGATATGCTCACCGGGTGTTTTTTGTAGCCTGGTAATAATTCCCTGAGCAATTTCATCGCT  >  NC_000913/477102‑477280

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