Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,246,452 T→C G301G (GGT→GGC)  lamB → phage lambda receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,246,4520TC99.4% 1469.4 / ‑3.6 460G301G (GGT→GGC) lamBphage lambda receptor protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (214/242);  total (217/243)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 482 total aligned reads displayed.

GTCGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGTGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAGTGGTG  >  NC_000913/4246362‑4246541
                                                                                          |                                                                                         
gTCGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGgcgc                                                                                         >  2:3088883/1‑93 (MQ=255)
gTCGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGgcgc                                                                                         >  2:4473569/1‑93 (MQ=255)
  cGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGa                                                                                       >  1:8085716/1‑93 (MQ=255)
  cGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGa                                                                                       >  1:2811491/1‑93 (MQ=255)
   gCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGAt                                                                                      >  1:4377882/1‑93 (MQ=255)
   gCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGAt                                                                                      <  2:7816616/93‑1 (MQ=255)
    cAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGAtc                                                                                     >  2:3831054/1‑93 (MQ=255)
       ggTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTcc                                                                                  <  1:4169743/93‑1 (MQ=255)
       ggTTCTGGCGTTGCATTTGATAAAGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTcc                                                                                  <  2:5600579/93‑1 (MQ=255)
        gTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCa                                                                                 <  1:6877637/93‑1 (MQ=255)
         ttCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCAt                                                                                >  1:780692/1‑93 (MQ=255)
         ttCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCAt                                                                                >  1:9944289/1‑93 (MQ=255)
            tGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCGTCGACCACGGTGCGATCTCCATggg                                                                             >  2:3231654/1‑93 (MQ=255)
             ggCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGc                                                                            >  1:6235091/1‑93 (MQ=255)
                gTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGaa                                                                         >  2:753771/1‑92 (MQ=255)
                gTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGAc                                                                         >  1:5158403/1‑93 (MQ=255)
                gTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGAc                                                                         >  2:1404083/1‑93 (MQ=255)
                 ttGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACa                                                                        >  1:2437711/1‑93 (MQ=255)
                 ttGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACa                                                                        <  1:2195485/93‑1 (MQ=255)
                  tGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTAGCCTCGACCACGGTGCGATCTCCATGGGCGACaa                                                                       >  2:5306554/1‑93 (MQ=255)
                     aTTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTg                                                                    >  2:1267757/1‑93 (MQ=255)
                       ttGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTggg                                                                  <  1:5906660/93‑1 (MQ=255)
                        tGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGa                                                                 <  1:4660141/93‑1 (MQ=255)
                         gATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGATGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGAc                                                                <  2:901450/93‑1 (MQ=255)
                          aTAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACa                                                               <  1:6599030/93‑1 (MQ=255)
                             aCGACAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACctga                                                            <  2:8932888/93‑1 (MQ=255)
                             aCGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACatga                                                            >  1:7362650/1‑93 (MQ=255)
                             aCGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACatga                                                            >  2:1109995/1‑93 (MQ=255)
                              cGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACatgat                                                           <  2:1649744/93‑1 (MQ=255)
                                 aaaaTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGtc                                                        >  2:5980815/1‑92 (MQ=255)
                                 aaaaTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTa                                                        <  1:924246/93‑1 (MQ=255)
                                  aaaTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTAc                                                       <  2:8725299/93‑1 (MQ=255)
                                  aaaTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTAc                                                       <  1:1759594/93‑1 (MQ=255)
                                   aaTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACg                                                      >  1:5716665/1‑93 (MQ=255)
                                    aTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGt                                                     <  2:6224441/93‑1 (MQ=255)
                                     tttGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTg                                                    >  2:8043361/1‑93 (MQ=255)
                                     tttGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTg                                                    >  1:6103169/1‑93 (MQ=255)
                                     tttGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTg                                                    >  1:3775503/1‑93 (MQ=255)
                                     tttGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTg                                                    <  2:4339633/93‑1 (MQ=255)
                                     tttGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTg                                                    <  2:7765536/93‑1 (MQ=255)
                                     tttGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTg                                                    >  1:191301/1‑93 (MQ=255)
                                        gCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTggga                                                 >  2:3975587/1‑92 (MQ=255)
                                         ccTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTa                                                <  2:2064081/93‑1 (MQ=255)
                                         ccTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTa                                                <  2:6844345/93‑1 (MQ=255)
                                          cTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTAt                                               >  2:8117329/1‑93 (MQ=255)
                                          cTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTAt                                               >  1:1144660/1‑93 (MQ=255)
                                           tACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATg                                              >  2:2473989/1‑93 (MQ=255)
                                           tACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATg                                              >  2:7326746/1‑93 (MQ=255)
                                           tACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATg                                              >  1:552388/1‑93 (MQ=255)
                                           tACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATg                                              <  1:673728/93‑1 (MQ=255)
                                           caCAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATg                                              <  2:9147685/92‑1 (MQ=255)
                                            aCAATGTCAACAACAACGGTCACATGCTGCGTATCCTAGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGt                                             <  1:3450049/93‑1 (MQ=255)
                                            aCAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGt                                             >  1:7793048/1‑93 (MQ=255)
                                             cAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTa                                            <  1:10008337/93‑1 (MQ=255)
                                              aaTGTCAACAATCACGGTCACATGCTTGGTGTCCTCGCCCACGGCGCGATCACCATGGGAGACAACTGGGACATGATGTACGTGGGTATATAc                                           <  2:479342/93‑1 (MQ=255)
                                              aaTATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTAc                                           <  1:7279391/93‑1 (MQ=255)
                                              aaTATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTAc                                           <  1:4921939/93‑1 (MQ=255)
                                               atatCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTAcc                                          <  1:7115119/93‑1 (MQ=255)
                                               atatCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTAcc                                          >  1:3041254/1‑93 (MQ=255)
                                               atatCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTAcc                                          <  2:2692548/93‑1 (MQ=255)
                                                tatCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCa                                         >  2:2601095/1‑93 (MQ=255)
                                                tatCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCa                                         >  2:3390812/1‑93 (MQ=255)
                                                  tCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAgg                                       <  2:9485153/93‑1 (MQ=255)
                                                     acaacaacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGata                                    <  2:1513487/93‑1 (MQ=255)
                                                     acaacaacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGata                                    <  1:10041423/93‑1 (MQ=255)
                                                      caacaacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGatat                                   >  2:3651084/1‑93 (MQ=255)
                                                         caacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCaa                                <  2:6348661/93‑1 (MQ=255)
                                                          aacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAAc                               >  2:1350758/1‑93 (MQ=255)
                                                          aacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAAc                               >  1:3344216/1‑93 (MQ=255)
                                                          aacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAAc                               >  1:1555566/1‑93 (MQ=255)
                                                          aacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAAc                               >  1:7715056/1‑93 (MQ=255)
                                                          aacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAAc                               <  1:1660677/93‑1 (MQ=255)
                                                             ggTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTgg                            >  2:5867591/1‑93 (MQ=255)
                                                              gTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTggg                           >  1:3914710/1‑93 (MQ=255)
                                                              gTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTggg                           <  1:7497245/93‑1 (MQ=255)
                                                               tCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGa                          >  2:30161/1‑93 (MQ=255)
                                                                 acaTGCTGCGGATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATa                        <  1:3272861/93‑1 (MQ=255)
                                                                  caTGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATaa                       >  1:141666/1‑93 (MQ=255)
                                                                  caTGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACGGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATaa                       <  1:7910757/93‑1 (MQ=255)
                                                                    tgctgcGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACg                     >  2:9328884/1‑93 (MQ=255)
                                                                    tgctgcGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACg                     >  1:5930493/1‑93 (MQ=255)
                                                                    tgctgcGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACg                     <  1:8517509/93‑1 (MQ=255)
                                                                       tgcGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACa                  >  2:4974073/1‑93 (MQ=255)
                                                                        gcGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACaa                 <  2:5789631/93‑1 (MQ=255)
                                                                        gcGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACaa                 <  1:9375615/93‑1 (MQ=255)
                                                                          gTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGAATGTACCAGGATATCAACTGGGATAACGACAACg               >  1:8348832/1‑93 (MQ=255)
                                                                          gTATCCTCAACCACGGAGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATACCAACTGGGATAACGACAACg               <  2:9681188/93‑1 (MQ=255)
                                                                           tATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACgg              >  2:5538306/1‑93 (MQ=255)
                                                                           tATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACgg              >  1:8049872/1‑93 (MQ=255)
                                                                           tATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACgg              >  1:2020337/1‑93 (MQ=255)
                                                                             tCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCa            >  1:8893392/1‑93 (MQ=255)
                                                                               cTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCAcc          >  2:2287266/1‑93 (MQ=255)
                                                                                tCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCa         >  2:7639144/1‑93 (MQ=255)
                                                                                 cGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCaa        <  1:5361394/93‑1 (MQ=255)
                                                                                  gACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACACCGGCACCAAg       <  2:6548818/93‑1 (MQ=255)
                                                                                   aCCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAgt      >  1:6462101/1‑93 (MQ=255)
                                                                                    ccACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAgtg     <  1:2675330/93‑1 (MQ=255)
                                                                                     cACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAgtgg    <  1:9067050/93‑1 (MQ=255)
                                                                                     cACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAgtgg    >  1:4482172/1‑93 (MQ=255)
                                                                                       cGGCGCGGTCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAgtggtg  <  1:3658059/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTCGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGTGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAGTGGTG  >  NC_000913/4246362‑4246541

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