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,5640AG99.8% 1463.3 / ‑5.1 453Y117C (TAT→TGT) yagUhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (274/178);  total (274/179)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 485 total aligned reads displayed.

CCAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTT  >  NC_000913/302474‑302653
                                                                                          |                                                                                         
ccAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGgttgtt                                                                                         >  2:7727632/1‑93 (MQ=255)
   aTGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTtgtg                                                                                      >  1:1148625/1‑93 (MQ=255)
   aTGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGGGTTTGCTGTCGGTTGTtgtg                                                                                      >  1:6268379/1‑93 (MQ=255)
     gCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTtgtgtg                                                                                    <  1:6351509/93‑1 (MQ=255)
       ggCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGt                                                                                  >  1:7140144/1‑93 (MQ=255)
       ggCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGt                                                                                  <  1:3153335/93‑1 (MQ=255)
        gCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTc                                                                                 >  1:8842048/1‑93 (MQ=255)
          tGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGc                                                                               >  2:1523303/1‑93 (MQ=255)
           gTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCt                                                                              <  1:4123928/93‑1 (MQ=255)
           gTTTATACCTTTGCCGGGCATGTCTTTAACTGGGGTGGTGTTACGCCCATTATCTTTTCGATAGTGTTTGCTCTCGGTTGTTGTGTGGTCGCt                                                                              >  1:1232015/1‑93 (MQ=255)
            tttATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTg                                                                             >  2:4556431/1‑93 (MQ=255)
             ttATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGa                                                                            <  1:7587099/93‑1 (MQ=255)
             ttATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGa                                                                            >  1:5360234/1‑93 (MQ=255)
              tataCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGaa                                                                           <  2:9859/93‑1 (MQ=255)
              tataCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGaa                                                                           <  1:9804351/93‑1 (MQ=255)
              tataCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGaa                                                                           >  1:2541161/1‑93 (MQ=255)
               ataCCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAg                                                                          >  1:1541060/1‑93 (MQ=255)
                    tttGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAttt                                                                     <  2:9526286/93‑1 (MQ=255)
                    tttGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTAttt                                                                     >  1:4524018/1‑93 (MQ=255)
                    tttGCCGGGCATGTCTTTAACCGGGTTGGTTTTGCGCCCATTATCTTTTCGATAGTGTTTGGTGTCTGTTGTTGTGTGGTCGCTGAAGTAttt                                                                     <  2:8220091/93‑1 (MQ=255)
                      tGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTcc                                                                   <  1:1015675/93‑1 (MQ=255)
                      tGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTcc                                                                   >  2:9229287/1‑93 (MQ=255)
                       gCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCa                                                                  <  1:4354236/93‑1 (MQ=255)
                         cGGGCATGTCTTTAACTGGGGTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaaa                                                                >  2:1842174/1‑93 (MQ=255)
                          gggCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTGGCTGAAGTATTTCCaaaa                                                               >  1:559658/1‑93 (MQ=255)
                           ggCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCaaaaa                                                              <  2:2550634/93‑1 (MQ=255)
                            gCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTTCAACAAt                                                             >  1:9290553/1‑93 (MQ=255)
                            gCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAAt                                                             <  2:1277671/93‑1 (MQ=255)
                             cATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAAtt                                                            >  2:2411261/1‑93 (MQ=255)
                              aTGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTa                                                           >  1:6029985/1‑93 (MQ=255)
                              aTGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTa                                                           <  2:7429543/93‑1 (MQ=255)
                               tGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTGGCTGAAGTATTTCCAAAAAttta                                                          >  2:8319656/1‑91 (MQ=255)
                               tGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTaa                                                          >  1:2437241/1‑93 (MQ=255)
                               tGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTaa                                                          >  2:2196536/1‑93 (MQ=255)
                               tGTCTTTAACTGGGTTGGTGATACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTaa                                                          >  1:6947848/1‑93 (MQ=255)
                                gTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTaaa                                                         <  1:8243917/93‑1 (MQ=255)
                                 tCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAAc                                                        <  2:66571/93‑1 (MQ=255)
                                 tCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAAc                                                        >  2:2033386/1‑93 (MQ=255)
                                 tCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAAc                                                        >  1:2790890/1‑93 (MQ=255)
                                 tCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAAc                                                        >  2:2959329/1‑93 (MQ=255)
                                 tCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAAc                                                        >  1:9958077/1‑93 (MQ=255)
                                   tttAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAActc                                                      <  1:211995/93‑1 (MQ=255)
                                    ttAACTCGGTTGGTGTTACGCGCATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAActct                                                     <  1:4692446/93‑1 (MQ=255)
                                       aCTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGc                                                  <  1:8470756/93‑1 (MQ=255)
                                         tGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAg                                                >  1:9672275/1‑93 (MQ=255)
                                          gggTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAgg                                               <  2:3319086/93‑1 (MQ=255)
                                             ttGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCt                                            >  2:5737671/1‑93 (MQ=255)
                                             ttGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTCGCAGGGCt                                            >  2:4500696/1‑93 (MQ=255)
                                              tGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAGCTCTGGCAGGTCtt                                           <  2:5104862/93‑1 (MQ=255)
                                                gtgtTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTAc                                         <  2:5592322/93‑1 (MQ=255)
                                                 tgtTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGTCTTACt                                        >  1:4952580/1‑93 (MQ=255)
                                                 tgtTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACt                                        >  1:1986344/1‑93 (MQ=255)
                                                 tgtTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACt                                        >  2:8517598/1‑93 (MQ=255)
                                                 tgtTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACt                                        >  1:9086227/1‑93 (MQ=255)
                                                     aCGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCa                                    >  1:878860/1‑93 (MQ=255)
                                                     aCGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCa                                    <  1:6574890/93‑1 (MQ=255)
                                                      cGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAg                                   >  1:2037883/1‑93 (MQ=255)
                                                      cGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAg                                   >  2:2106628/1‑93 (MQ=255)
                                                        cacaTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAggg                                 >  2:1695725/1‑92 (MQ=255)
                                                         acaTTATCTTTTCGATAGTGTTTGCTGTCGGTTGATGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTg                                >  2:6377946/1‑93 (MQ=255)
                                                          caTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGc                               >  1:517579/1‑93 (MQ=255)
                                                          caTTAGCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGc                               <  2:619954/93‑1 (MQ=255)
                                                           aTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCt                              >  2:4028793/1‑93 (MQ=255)
                                                           aTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCt                              >  1:6166152/1‑93 (MQ=255)
                                                           aTTATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCt                              >  2:1093466/1‑93 (MQ=255)
                                                            ttATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTCTTTCCAAAAATTAACCTCTGGCAGGGCTTACTGGCAGGTGcct                             >  1:5233666/1‑91 (MQ=255)
                                                            ttATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCtt                             >  1:6699380/1‑93 (MQ=255)
                                                            ttATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCtt                             >  1:573449/1‑93 (MQ=255)
                                                            ttATCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCtt                             >  2:3031651/1‑93 (MQ=255)
                                                              aTCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGTCAGGTGCTTTa                           >  2:5209371/1‑93 (MQ=255)
                                                              aTCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGGGCTTTa                           >  1:7430569/1‑93 (MQ=255)
                                                               tCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCtttta                          >  2:7574111/1‑91 (MQ=255)
                                                               tCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg                          >  1:1797605/1‑93 (MQ=255)
                                                               tCTTTTCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAg                          >  2:8055293/1‑93 (MQ=255)
                                                                   ttCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTTGCCCa                      >  1:2698395/1‑93 (MQ=255)
                                                                   ttCGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCa                      <  1:5524032/93‑1 (MQ=255)
                                                                     cGATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGCGCTTTAGCCCaat                    >  2:3507224/1‑92 (MQ=255)
                                                                      gATAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACt                   <  1:8679881/93‑1 (MQ=255)
                                                                        tAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACttt                 <  2:7195142/93‑1 (MQ=255)
                                                                        tAGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACttt                 >  1:7982459/1‑93 (MQ=255)
                                                                         aGTGTTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACtttt                <  2:6912918/93‑1 (MQ=255)
                                                                          gtgtTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACttttt               <  2:9096522/93‑1 (MQ=255)
                                                                          gtgtTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACttttt               <  2:3132120/93‑1 (MQ=255)
                                                                            gtTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGt             <  2:9657431/93‑1 (MQ=255)
                                                                            gtTTGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGt             <  2:4711309/93‑1 (MQ=255)
                                                                             tttGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAAATTTTTGtt            >  1:3757083/1‑93 (MQ=255)
                                                                              tcgCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTc           >  1:9404416/3‑93 (MQ=255)
                                                                               tGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCCACTTTTTGTTCa          >  2:5451780/1‑93 (MQ=255)
                                                                               tGCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCa          <  1:3928861/93‑1 (MQ=255)
                                                                                gCTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCat         >  2:8723748/1‑93 (MQ=255)
                                                                                 cTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTTTTtata        >  1:2224472/1‑93 (MQ=255)
                                                                                 cTGTCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCata        <  1:192451/93‑1 (MQ=255)
                                                                                    tCGTTTGTTGTGTGGTCGCTGAAGTATTGCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGa     <  1:3223464/93‑1 (MQ=255)
                                                                                    tCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATgt     >  1:7774038/1‑92 (MQ=255)
                                                                                    tCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATgt     >  1:9579558/1‑92 (MQ=255)
                                                                                    tCGGTTGTTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGa     <  2:5871732/93‑1 (MQ=255)
                                                                                    tCGGTTGTTGTGTAGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGa     <  2:2805106/93‑1 (MQ=255)
                                                                                       gttgttGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAttt  >  1:587434/1‑93 (MQ=255)
                                                                                       gttgttGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAttt  >  1:2976945/1‑93 (MQ=255)
                                                                                       gttgttGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGAttt  >  1:3508256/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCAATGCGGCTGTTTATACCTTTGCCGGGCATGTCTTTAACTGGGTTGGTGTTACGCACATTATCTTTTCGATAGTGTTTGCTGTCGGTTATTGTGTGGTCGCTGAAGTATTTCCAAAAATTAAACTCTGGCAGGGCTTACTGGCAGGTGCTTTAGCCCAACTTTTTGTTCATATGATTT  >  NC_000913/302474‑302653

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