Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,450,334 G→A A64T (GCT→ACT)  ytfT → putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,450,3340GA99.6% 819.9 / NA 234A64T (GCT→ACT) ytfTputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (110/123);  total (110/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 235 total aligned reads displayed.

GTCGATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAGCTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCGCC  >  NC_000913/4450244‑4450423
                                                                                          |                                                                                         
gTCGATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACt                                                                                         <  1:3555573/93‑1 (MQ=255)
 tCGATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTc                                                                                        >  2:3585003/1‑93 (MQ=255)
  cGATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTcc                                                                                       <  2:4683632/93‑1 (MQ=255)
   gATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTccc                                                                                      <  2:5820048/93‑1 (MQ=255)
   gATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTccc                                                                                      <  1:5376914/93‑1 (MQ=255)
    aTAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCg                                                                                     <  1:5331902/93‑1 (MQ=255)
     tAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGt                                                                                    <  1:3318216/93‑1 (MQ=255)
     tAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGt                                                                                    <  2:587721/93‑1 (MQ=255)
      aGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGtt                                                                                   >  2:331623/1‑93 (MQ=255)
      aGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGtt                                                                                   <  2:6373781/93‑1 (MQ=255)
       gCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTg                                                                                  >  2:1048336/1‑93 (MQ=255)
       gCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTg                                                                                  <  1:73611/93‑1 (MQ=255)
          tGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGt                                                                               <  1:2391065/93‑1 (MQ=255)
          tGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGt                                                                               <  2:2905228/93‑1 (MQ=255)
           ggTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCAGTTGCGtt                                                                              <  1:6553020/93‑1 (MQ=255)
             tAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTAc                                                                            <  2:5090328/93‑1 (MQ=255)
              aGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACt                                                                           >  2:2600610/1‑93 (MQ=255)
               gCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTg                                                                          <  1:2629966/93‑1 (MQ=255)
               gCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTg                                                                          <  2:5887494/93‑1 (MQ=255)
                 cccGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGc                                                                        >  2:1901707/1‑93 (MQ=255)
                 cccGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGc                                                                        >  2:4400056/1‑93 (MQ=255)
                  ccGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCg                                                                       <  1:6048114/93‑1 (MQ=255)
                  ccGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCg                                                                       >  1:4057043/1‑93 (MQ=255)
                   cGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGa                                                                      <  1:6351538/93‑1 (MQ=255)
                     cATTTCTGGCAGGTGGTGCTTTAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGAtt                                                                    >  1:843511/1‑93 (MQ=255)
                     cATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGAtt                                                                    >  2:1425839/1‑93 (MQ=255)
                     cATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGAtt                                                                    >  1:4329450/1‑93 (MQ=255)
                     cATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGAtt                                                                    <  2:5535189/93‑1 (MQ=255)
                      aTTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTg                                                                   >  1:1572298/1‑93 (MQ=255)
                         tCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCa                                                                <  2:35076/93‑1 (MQ=255)
                         tCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCa                                                                <  1:390243/93‑1 (MQ=255)
                          cTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCAt                                                               <  1:5807930/93‑1 (MQ=255)
                            ggCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGa                                                             >  2:1768136/1‑93 (MQ=255)
                               aGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGc                                                          >  1:6166630/1‑93 (MQ=255)
                               aGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGc                                                          <  2:4292233/93‑1 (MQ=255)
                                ggtggtGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCt                                                         >  2:2950582/1‑93 (MQ=255)
                                ggtggtGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCt                                                         >  2:4483983/1‑93 (MQ=255)
                                  tggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTgg                                                       <  2:4566292/93‑1 (MQ=255)
                                  tggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTgg                                                       >  2:5979391/1‑93 (MQ=255)
                                   ggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGt                                                      >  1:4918868/1‑93 (MQ=255)
                                   ggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGt                                                      >  2:3467163/1‑93 (MQ=255)
                                   ggtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGt                                                      <  1:6253697/93‑1 (MQ=255)
                                    gtgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTg                                                     <  2:3969799/93‑1 (MQ=255)
                                     tgCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGa                                                    <  2:2802634/93‑1 (MQ=255)
                                      gCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGAt                                                   <  2:2174588/93‑1 (MQ=255)
                                      gCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGAt                                                   >  2:5471978/1‑93 (MQ=255)
                                       cTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATc                                                  >  1:2531764/1‑93 (MQ=255)
                                       cTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATc                                                  <  1:1111779/93‑1 (MQ=255)
                                         tCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGc                                                >  2:1546770/1‑93 (MQ=255)
                                         tCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGc                                                <  2:6664165/93‑1 (MQ=255)
                                           aGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCa                                              >  1:1427150/1‑93 (MQ=255)
                                           aGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCa                                              <  1:2810994/93‑1 (MQ=255)
                                           aGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCa                                              >  1:4964699/1‑93 (MQ=255)
                                               tGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCgg                                          <  1:3114553/93‑1 (MQ=255)
                                                gggCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGt                                         <  2:6706659/93‑1 (MQ=255)
                                                gggCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGt                                         <  1:2155313/93‑1 (MQ=255)
                                                   cGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTggg                                      <  2:3237878/93‑1 (MQ=255)
                                                    gTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGa                                     >  1:3666869/1‑93 (MQ=255)
                                                      ttGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGAtt                                   <  2:3359425/93‑1 (MQ=255)
                                                        gTTCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGa                                 <  1:3082533/93‑1 (MQ=255)
                                                        gTTCGGTAGCCCCATAGACATTCGTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGa                                 <  1:2825554/93‑1 (MQ=255)
                                                         ttCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAt                                >  2:575590/1‑93 (MQ=255)
                                                         ttCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAt                                <  2:264548/93‑1 (MQ=255)
                                                         ttCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAt                                >  1:3367003/1‑93 (MQ=255)
                                                         ttCGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTAAt                                <  2:4982649/93‑1 (MQ=255)
                                                           cGGTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtct                              >  1:1828312/1‑93 (MQ=255)
                                                            ggTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtctc                             <  2:4955422/93‑1 (MQ=255)
                                                             gTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtctct                            <  1:2141983/93‑1 (MQ=255)
                                                             gTAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGAGGGATTGAtctct                            >  1:4012388/1‑93 (MQ=255)
                                                              tAGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGAtctctc                           <  2:434458/93‑1 (MQ=255)
                                                               aGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTcc                          >  2:2512440/1‑93 (MQ=255)
                                                               aGCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATAGCCACCGGTGGGATTGATCTCTcc                          >  1:5526063/1‑93 (MQ=255)
                                                                gCCCCATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCg                         <  2:6480102/93‑1 (MQ=255)
                                                                 ccccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGt                        >  1:6413433/1‑93 (MQ=255)
                                                                   ccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAg                      <  1:4389797/93‑1 (MQ=255)
                                                                   ccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAg                      >  1:6101540/1‑93 (MQ=255)
                                                                   ccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAg                      <  2:571607/93‑1 (MQ=255)
                                                                   ccATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCAACGGTGGGATTGATCTCTCCGTAg                      >  1:5414639/1‑93 (MQ=255)
                                                                    cATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgg                     <  1:941903/93‑1 (MQ=255)
                                                                    cATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgg                     <  2:1480913/93‑1 (MQ=255)
                                                                    cATAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgg                     <  1:2117483/93‑1 (MQ=255)
                                                                     aTAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAggg                    >  2:6097996/1‑93 (MQ=255)
                                                                      tAGACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAgggg                   >  1:1543516/1‑93 (MQ=255)
                                                                        gACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCg                 <  1:3027451/93‑1 (MQ=255)
                                                                        gACATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCg                 <  1:5150512/93‑1 (MQ=255)
                                                                        gACATTCTTAACCGTGCAACTCCCGTTGCGTTACCGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCg                 >  1:6151830/1‑93 (MQ=255)
                                                                          cATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGt               <  2:6597159/93‑1 (MQ=255)
                                                                          cATTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGt               <  2:6034917/93‑1 (MQ=255)
                                                                           aTTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCgggg              >  1:1671806/1‑91 (MQ=255)
                                                                           aTTCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTg              <  1:4777399/93‑1 (MQ=255)
                                                                            ttCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGa             <  2:1282090/93‑1 (MQ=255)
                                                                             tCTTAACCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGAt            <  2:543501/93‑1 (MQ=255)
                                                                                  aCCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCtt       >  1:4221917/1‑92 (MQ=255)
                                                                                  aCCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTa       >  2:3897786/1‑93 (MQ=255)
                                                                                  aCCGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGGGATGGCTa       >  2:2013584/1‑93 (MQ=255)
                                                                                    cGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATc     <  1:5215385/93‑1 (MQ=255)
                                                                                    cGTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATc     <  1:270066/93‑1 (MQ=255)
                                                                                     gTGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCg    >  2:1973594/1‑93 (MQ=255)
                                                                                      tGCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCGc   >  2:5667797/1‑93 (MQ=255)
                                                                                       gCAACTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCGcc  <  1:3593700/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTCGATAGCCTGGTAGCCCCGCATTTCTGGCAGGTGGTGCTTCAGGATGGGCGTTTGTTCGGTAGCCCCATAGACATTCTTAACCGTGCAGCTCCCGTTGCGTTACTGGCGATTGGCATGACGCTGGTGATCGCCACCGGTGGGATTGATCTCTCCGTAGGGGCGGTGATGGCTATCGCC  >  NC_000913/4450244‑4450423

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