Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 302,564 A→G Y117C (TAT→TGT)  yagU → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913302,5640AG98.8% 518.2 / ‑4.0 162Y117C (TAT→TGT) yagUhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (97/63);  total (97/65)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 170 total aligned reads displayed.

CCAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTT  >  NC_000913/302474‑302653
                                                                                          |                                                                                         
ccAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGgttgtt                                                                                         <  2:3153008/93‑1 (MQ=255)
  aaTGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTtgttgt                                                                                       <  1:4740542/93‑1 (MQ=255)
   aTGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTtgtg                                                                                      >  1:3588597/1‑93 (MQ=255)
    tGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTtgtgt                                                                                     <  2:147132/93‑1 (MQ=255)
    tGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTtgtgt                                                                                     <  1:4195826/93‑1 (MQ=255)
     gCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTtgtgtg                                                                                    <  2:2937857/93‑1 (MQ=255)
     gCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTtgtgtg                                                                                    >  1:2069181/1‑93 (MQ=255)
       ggCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGt                                                                                  >  1:1009936/1‑93 (MQ=255)
       ggCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGt                                                                                  >  1:2482417/1‑93 (MQ=255)
       ggCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGt                                                                                  <  2:4144124/93‑1 (MQ=255)
        gCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTc                                                                                 <  1:1595455/93‑1 (MQ=255)
         cTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCg                                                                                >  1:2897789/1‑93 (MQ=255)
           gTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCt                                                                              >  1:3923909/1‑93 (MQ=255)
            tttATCCCTTTGCCGGGCATGTCTTTAGCTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTATTGTGTGGTCGCTg                                                                             <  1:3109432/93‑1 (MQ=255)
            gttATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTGGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTg                                                                             <  1:1319151/92‑1 (MQ=255)
             ttATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGa                                                                            >  1:3807626/1‑93 (MQ=255)
              tataCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGaa                                                                           <  1:1421025/93‑1 (MQ=255)
              tataCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGaa                                                                           >  1:2390519/1‑93 (MQ=255)
                 aCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTa                                                                        <  2:1494613/93‑1 (MQ=255)
                  ccTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAt                                                                       >  1:775584/1‑93 (MQ=255)
                  ccTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAt                                                                       <  1:991927/93‑1 (MQ=255)
                   cTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAtt                                                                      >  1:2249205/1‑93 (MQ=255)
                    tttGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAttt                                                                     >  1:2711872/1‑93 (MQ=255)
                    tttGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAttt                                                                     >  1:4558896/1‑93 (MQ=255)
                    tttGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAttt                                                                     >  1:4034327/1‑93 (MQ=255)
                    tttGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAttt                                                                     >  2:2463427/1‑93 (MQ=255)
                     ttGCCGTGCATGTCTTTAACTGTGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTc                                                                    <  1:2246468/93‑1 (MQ=255)
                      tGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTcc                                                                   <  2:58778/93‑1 (MQ=255)
                        ccGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaa                                                                 >  1:1158268/1‑93 (MQ=255)
                         cGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaaa                                                                <  2:1396037/93‑1 (MQ=255)
                         cGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaaa                                                                >  2:3242092/1‑93 (MQ=255)
                          gggCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaaaa                                                               <  1:4332382/93‑1 (MQ=255)
                          gggCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaaaa                                                               >  1:3387860/1‑93 (MQ=255)
                          gggCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaaaa                                                               >  2:1967974/1‑93 (MQ=255)
                           ggCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaaaaa                                                              >  2:1451881/1‑93 (MQ=255)
                           ggCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaaaaa                                                              >  2:3602458/1‑93 (MQ=255)
                            gCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAAt                                                             >  1:812092/1‑93 (MQ=255)
                            gCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAAt                                                             <  1:3062234/93‑1 (MQ=255)
                               tGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTaa                                                          >  1:3778810/1‑93 (MQ=255)
                               tGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTaa                                                          >  2:480754/1‑93 (MQ=255)
                               tGTCTTGAACTGGGTTGGTGTTACGTACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTaa                                                          <  2:828330/93‑1 (MQ=255)
                                gTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTaaa                                                         >  1:1721988/1‑93 (MQ=255)
                                 tCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAAc                                                        >  2:3847148/1‑93 (MQ=255)
                                  cTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAAct                                                       <  2:2796314/93‑1 (MQ=255)
                                   tttAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAActc                                                      >  2:2434360/1‑93 (MQ=255)
                                    ttAGCTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAActct                                                     <  1:3531246/93‑1 (MQ=255)
                                    ttAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAActct                                                     <  2:3084745/93‑1 (MQ=255)
                                     tAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTg                                                    <  2:1849050/93‑1 (MQ=255)
                                      aaCTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTgg                                                   >  2:3427240/1‑93 (MQ=255)
                                          gggTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAgg                                               >  1:3408322/1‑93 (MQ=255)
                                          gggTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAgg                                               >  2:12481/1‑93 (MQ=255)
                                           ggTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAggg                                              <  1:4235535/93‑1 (MQ=255)
                                           ggTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAggg                                              <  1:3967347/93‑1 (MQ=255)
                                             ttGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCt                                            >  2:1249418/1‑93 (MQ=255)
                                             ttGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCt                                            <  1:1812074/93‑1 (MQ=255)
                                               ggTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTa                                          <  1:1911180/93‑1 (MQ=255)
                                                gtgtTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTAc                                         <  2:1713341/93‑1 (MQ=255)
                                                gtgtTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTAc                                         <  2:164058/93‑1 (MQ=255)
                                                 tgtTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACt                                        <  2:674241/93‑1 (MQ=255)
                                                     aCGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCa                                    >  1:1070689/1‑93 (MQ=255)
                                                      cGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAg                                   <  1:4503461/93‑1 (MQ=255)
                                                      cGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAg                                   <  1:711340/93‑1 (MQ=255)
                                                       gCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAgg                                  >  1:3591793/1‑93 (MQ=255)
                                                        cacaTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAggg                                 >  1:2733164/1‑92 (MQ=255)
                                                        cacaTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAggg                                 >  2:239123/1‑92 (MQ=255)
                                                         acaTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTg                                >  2:1151326/1‑93 (MQ=255)
                                                          caTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTTCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGc                               >  1:4248138/1‑93 (MQ=255)
                                                          caTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGc                               <  1:547272/93‑1 (MQ=255)
                                                           aTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCt                              >  2:3779859/1‑93 (MQ=255)
                                                           aTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAATTCCAAAAATTAAACTCTGTCAGCGCTTACTGGCAGGTgtt                              >  2:1513203/1‑91 (MQ=255)
                                                            ttATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCtt                             <  1:2757740/93‑1 (MQ=255)
                                                              aTCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTa                           >  2:1818083/1‑93 (MQ=255)
                                                              aTCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTa                           >  2:2063893/1‑93 (MQ=255)
                                                              aTCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTa                           >  2:345640/1‑93 (MQ=255)
                                                              aTCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTa                           >  2:3406543/1‑93 (MQ=255)
                                                               tCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTac                          >  2:4181137/1‑92 (MQ=255)
                                                               tCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg                          >  1:1316757/1‑93 (MQ=255)
                                                                cTTTTCGATAGTGTTTGCTGTCGGTTGGTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGTCAGCGGCTTTTGc                         >  1:2427838/1‑93 (MQ=255)
                                                                 ttttCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGcc                        >  1:3240629/1‑93 (MQ=255)
                                                                    tCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGcccca                     >  2:2313949/1‑91 (MQ=255)
                                                                    tCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCaa                     >  2:3068999/1‑93 (MQ=255)
                                                                      gATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACt                   <  2:3917709/93‑1 (MQ=255)
                                                                      gATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACt                   >  2:3601875/1‑93 (MQ=255)
                                                                       aTAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACtt                  <  2:540712/93‑1 (MQ=255)
                                                                        tAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACttt                 >  1:2400342/1‑93 (MQ=255)
                                                                         aGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCCTTAGCCCAACtttt                >  2:4339245/1‑93 (MQ=255)
                                                                         aGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACtttt                >  1:3661790/1‑93 (MQ=255)
                                                                           tgtTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACtttttt              >  1:3877262/1‑92 (MQ=255)
                                                                           tgtTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTg              <  2:4309712/93‑1 (MQ=255)
                                                                              ttGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTc           >  1:4405557/1‑93 (MQ=255)
                                                                                 cTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCata        <  1:1324813/93‑1 (MQ=255)
                                                                                  tGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCatat       <  2:4083022/93‑1 (MQ=255)
                                                                                   gTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTAGCCCAACTTTTTGTTTATATg      >  2:1030645/1‑93 (MQ=255)
                                                                                    tCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTGTTCATATGa     >  1:2268008/1‑93 (MQ=255)
                                                                                    tCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTTTTCATATGa     >  2:1785433/1‑93 (MQ=255)
                                                                                     cGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAt    <  1:1456142/93‑1 (MQ=255)
                                                                                      ggTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTGTTCATATGAtt   >  1:4716925/1‑93 (MQ=255)
                                                                                      ggTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtt   >  2:4341581/1‑93 (MQ=255)
                                                                                      ggTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAtt   <  2:3726084/93‑1 (MQ=255)
                                                                                       gttgttGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAttt  >  1:640400/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTT  >  NC_000913/302474‑302653

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: