Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,246,188 A→G E213E (GAA→GAG lamB → phage lambda receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,246,1880AG99.5% 1376.1 / ‑4.5 423E213E (GAA→GAGlamBphage lambda receptor protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (235/186);  total (235/188)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 436 total aligned reads displayed.

AACAATATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGA  >  NC_000913/4246096‑4246277
                                                                                            |                                                                                         
  caaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCt                                                                                         >  2:4547365/1‑93 (MQ=255)
  caaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCt                                                                                         >  2:7674312/1‑93 (MQ=255)
   aaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTg                                                                                        >  2:4207488/1‑93 (MQ=255)
   aaTATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTg                                                                                        >  1:4789190/1‑93 (MQ=255)
   aaTATTTATGACTACACCACCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCACCCCGGGCGGCACATTAGAGCTg                                                                                        <  2:573129/93‑1 (MQ=255)
    atatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTgg                                                                                       >  1:8859026/1‑93 (MQ=255)
    atatTTATGACTATACCAACGAAACCGCGAACCACGGTTTCCATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCttg                                                                                       >  2:8998255/1‑91 (MQ=255)
     tatTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTggg                                                                                      <  2:8387208/93‑1 (MQ=255)
     tatTTATGACTACACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTggg                                                                                      <  2:3771930/93‑1 (MQ=255)
      atTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgt                                                                                     <  1:1665383/93‑1 (MQ=255)
       tttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtg                                                                                    <  2:7044802/93‑1 (MQ=255)
       ggtATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGAGGGAAATCAACCCGGGCGGCACAACAGAGCTGGgtg                                                                                    >  1:7710922/3‑93 (MQ=255)
        ttATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGgtgt                                                                                   <  2:4753447/93‑1 (MQ=255)
         tATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTc                                                                                  <  2:535204/93‑1 (MQ=255)
            gACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGAc                                                                               <  1:5726401/93‑1 (MQ=255)
              cTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTa                                                                             <  2:4440586/93‑1 (MQ=255)
               tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCACACCGGGCGGCACATTAGAGCTGGGTGTCGACTAc                                                                            >  1:9681189/1‑93 (MQ=255)
               tataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGCGGACTAc                                                                            >  2:2672230/1‑93 (MQ=255)
                ataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACg                                                                           >  2:8635597/1‑93 (MQ=255)
                ataCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACg                                                                           >  1:8036982/1‑93 (MQ=255)
                 taCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACgg                                                                          >  1:3367111/1‑93 (MQ=255)
                 taCCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGGGTCGCCTACgg                                                                          >  1:7330102/1‑93 (MQ=255)
                   ccAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTc                                                                        >  1:3929072/1‑93 (MQ=255)
                    cAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCg                                                                       <  1:4199473/93‑1 (MQ=255)
                    cAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGCCGGCACATTAGAGCTGGGTGTCGACTAAGGACg                                                                       >  1:1972050/1‑93 (MQ=255)
                      aCGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTg                                                                     >  1:8078185/1‑93 (MQ=255)
                        gAATCCGCGATCTACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGcc                                                                   <  1:9462832/93‑1 (MQ=255)
                        gAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGcc                                                                   >  2:2117600/1‑93 (MQ=255)
                        gAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGcc                                                                   >  2:6141071/1‑93 (MQ=255)
                         aaaCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCa                                                                  >  1:9880326/1‑93 (MQ=255)
                         aaaCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCa                                                                  >  1:6945712/1‑93 (MQ=255)
                          aaCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCaa                                                                 >  2:1347884/1‑93 (MQ=255)
                           aCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAAc                                                                <  2:6919761/93‑1 (MQ=255)
                           aCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAAc                                                                >  2:6641748/1‑93 (MQ=255)
                           aCCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAAc                                                                <  2:283407/93‑1 (MQ=255)
                            ccGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACt                                                               <  1:8592446/93‑1 (MQ=255)
                            ccGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACt                                                               >  1:6763329/1‑93 (MQ=255)
                             cgcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACtt                                                              >  1:3611189/1‑93 (MQ=255)
                              gcgAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTg                                                             <  2:1588541/93‑1 (MQ=255)
                                 aaCGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGt                                                          >  1:10016565/1‑93 (MQ=255)
                                     acgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATa                                                      >  2:2896211/1‑93 (MQ=255)
                                     acgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATa                                                      >  2:3663965/1‑93 (MQ=255)
                                     acgTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATa                                                      >  2:1029056/1‑93 (MQ=255)
                                      cgTTTTCGATGTGCGTTGAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATaa                                                     <  2:6410837/93‑1 (MQ=255)
                                       gTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAAc                                                    <  1:5379644/93‑1 (MQ=255)
                                        ttttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACt                                                   >  2:3498871/1‑93 (MQ=255)
                                        ttttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACt                                                   >  1:490053/1‑93 (MQ=255)
                                         tttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTa                                                  <  1:7128662/93‑1 (MQ=255)
                                          ttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATCACTAt                                                 >  1:290469/1‑93 (MQ=255)
                                          ttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTAt                                                 >  1:5839614/1‑93 (MQ=255)
                                          ttCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTAt                                                 <  1:2573723/93‑1 (MQ=255)
                                           tCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATc                                                >  1:1277870/1‑93 (MQ=255)
                                            cGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCg                                               >  2:3924780/1‑93 (MQ=255)
                                              aTGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                             <  1:7790486/93‑1 (MQ=255)
                                              aTGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                             <  2:3170958/93‑1 (MQ=255)
                                              aTGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                             >  1:1426535/1‑93 (MQ=255)
                                              aTGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTc                                             <  2:9601591/93‑1 (MQ=255)
                                               tgtgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCt                                            >  2:4940342/1‑93 (MQ=255)
                                                gggCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTg                                           <  2:1893403/91‑1 (MQ=255)
                                                 tgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTgg                                          >  2:9696281/1‑93 (MQ=255)
                                                 tgCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTgg                                          >  1:5129507/1‑93 (MQ=255)
                                                  gCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGt                                         <  1:1503187/93‑1 (MQ=255)
                                                     tttAGCGCAGATGGAAATCAGCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGa                                      <  2:7302527/93‑1 (MQ=255)
                                                     tttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGa                                      <  1:1122203/93‑1 (MQ=255)
                                                     tttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGa                                      >  2:5145023/1‑93 (MQ=255)
                                                     tttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGa                                      <  2:2014451/93‑1 (MQ=255)
                                                     tttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGa                                      <  1:456829/93‑1 (MQ=255)
                                                     tttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGa                                      >  1:6276348/1‑93 (MQ=255)
                                                      ttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGAt                                     <  1:4612025/93‑1 (MQ=255)
                                                      ttAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGAt                                     <  2:5880333/93‑1 (MQ=255)
                                                       tAGCGCAGAGGGAAATCATGCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATg                                    <  2:118013/93‑1 (MQ=255)
                                                         gcgcAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGgc                                  <  2:9876186/93‑1 (MQ=255)
                                                         gcgcAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGgc                                  >  1:7886260/1‑93 (MQ=255)
                                                           gcAGATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGgcgc                                <  2:4294711/93‑1 (MQ=255)
                                                              gATGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATc                             >  1:8137502/1‑93 (MQ=255)
                                                               aTGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCg                            >  2:924246/1‑93 (MQ=255)
                                                               aTGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCg                            >  2:7574741/1‑93 (MQ=255)
                                                               aTGGAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCg                            <  1:4504688/93‑1 (MQ=255)
                                                               aTGCAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCg                            <  2:2533424/93‑1 (MQ=255)
                                                                 ggAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGaa                          >  2:8111746/1‑93 (MQ=255)
                                                                  gAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGaaa                         <  1:6550172/93‑1 (MQ=255)
                                                                  gAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGaaa                         <  2:6797678/93‑1 (MQ=255)
                                                                  gAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGaaa                         <  1:5884917/93‑1 (MQ=255)
                                                                  gAAATCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGaaa                         <  1:8346140/93‑1 (MQ=255)
                                                                    aaTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGa                       <  1:1769955/93‑1 (MQ=255)
                                                                    aaTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGa                       <  1:432049/93‑1 (MQ=255)
                                                                    aaTCAACCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGa                       <  1:2598486/93‑1 (MQ=255)
                                                                      tCAACCCGGGCGGCACATTAGAGCTGGGTGGCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACg                     >  1:9182746/1‑93 (MQ=255)
                                                                        aaCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTCGTTGATGGCGCATCGAAAGACGGc                   >  1:642736/1‑93 (MQ=255)
                                                                         aCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCt                  >  1:3027272/1‑93 (MQ=255)
                                                                         aCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCt                  >  1:2361122/1‑93 (MQ=255)
                                                                         aCCCGGGCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCt                  <  1:1164105/93‑1 (MQ=255)
                                                                             gggCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGtt              >  1:3999676/1‑93 (MQ=255)
                                                                             gggCGGCACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGtt              >  2:8816504/1‑93 (MQ=255)
                                                                              ggcggcACATTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTa             <  2:8287047/93‑1 (MQ=255)
                                                                                    acaTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCAct       >  2:9341798/1‑93 (MQ=255)
                                                                                     caTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActg      >  1:6091488/1‑93 (MQ=255)
                                                                                     caTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActg      >  2:6279976/1‑93 (MQ=255)
                                                                                      aTTAGAGCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCActgc     <  1:5492197/93‑1 (MQ=255)
                                                                                         agagCTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGa  >  2:7497244/1‑93 (MQ=255)
                                                                                            |                                                                                         
AACAATATTTATGACTATACCAACGAAACCGCGAACGACGTTTTCGATGTGCGTTTAGCGCAGATGGAAATCAACCCGGGCGGCACATTAGAACTGGGTGTCGACTACGGTCGTGCCAACTTGCGTGATAACTATCGTCTGGTTGATGGCGCATCGAAAGACGGCTGGTTATTCACTGCTGA  >  NC_000913/4246096‑4246277

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