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.7% 947.5 / NA 294G301G (GGT→GGC) lamBphage lambda receptor protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (146/147);  minor base A (0/1);  total (146/148)
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 308 total aligned reads displayed.

GTCGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGTGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAGTGGTG  >  NC_000913/4246362‑4246541
                                                                                          |                                                                                         
gTCGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGgcgc                                                                                         <  1:1041498/93‑1 (MQ=255)
gTCGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGgcgc                                                                                         <  1:458705/93‑1 (MQ=255)
gTCGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGgcgc                                                                                         <  1:5483761/93‑1 (MQ=255)
 tCGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGcgcg                                                                                        >  2:738569/1‑93 (MQ=255)
   gCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGAt                                                                                      >  2:5595967/1‑93 (MQ=255)
    cAGGTTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGAtc                                                                                     <  1:2336546/93‑1 (MQ=255)
     aGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGAtct                                                                                    >  1:3854913/1‑93 (MQ=255)
      gggTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGAtctc                                                                                   <  1:2224608/93‑1 (MQ=255)
      gggTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGAtctc                                                                                   >  2:2782450/1‑93 (MQ=255)
      gggTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGAtctc                                                                                   >  2:2954351/1‑93 (MQ=255)
        gTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCa                                                                                 >  2:5712605/1‑93 (MQ=255)
          tCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATg                                                                               <  1:962145/93‑1 (MQ=255)
            tGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATggg                                                                             <  1:6093448/93‑1 (MQ=255)
             ggCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGc                                                                            <  1:3529509/93‑1 (MQ=255)
               cGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGa                                                                          >  1:3443976/1‑93 (MQ=255)
                gTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGAc                                                                         >  1:1539252/1‑93 (MQ=255)
                gTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGAc                                                                         >  1:2978904/1‑93 (MQ=255)
                 ttGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACa                                                                        >  2:4117073/1‑93 (MQ=255)
                  tGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACaa                                                                       >  1:1931266/1‑93 (MQ=255)
                  tGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACaa                                                                       <  2:5328581/93‑1 (MQ=255)
                   gCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAAc                                                                      <  2:3589656/93‑1 (MQ=255)
                   gCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAAc                                                                      <  1:1236238/93‑1 (MQ=255)
                   gCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAAc                                                                      >  1:2816442/1‑93 (MQ=255)
                   gCATTTGATAACGAAAAATGTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAAc                                                                      <  1:262920/93‑1 (MQ=255)
                    cATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACt                                                                     >  2:6272094/1‑93 (MQ=255)
                    cATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACt                                                                     >  1:748678/1‑93 (MQ=255)
                     aTTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTg                                                                    >  1:132663/1‑93 (MQ=255)
                     aTTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTg                                                                    >  1:3607853/1‑93 (MQ=255)
                     aTTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTg                                                                    <  2:5295012/93‑1 (MQ=255)
                      tttGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTgg                                                                   >  1:624890/1‑93 (MQ=255)
                      tttGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTgg                                                                   >  1:2945589/1‑93 (MQ=255)
                        tGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGa                                                                 <  2:1206350/93‑1 (MQ=255)
                        tGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGa                                                                 <  1:5531843/93‑1 (MQ=255)
                         gATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGAc                                                                >  1:5116459/1‑93 (MQ=255)
                         gATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGAc                                                                <  1:5994806/93‑1 (MQ=255)
                          aTAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACa                                                               <  2:2271396/93‑1 (MQ=255)
                           tAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACat                                                              >  1:4094494/1‑93 (MQ=255)
                            aaCGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACatg                                                             >  2:3974185/1‑93 (MQ=255)
                             aCGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACatga                                                            <  2:1718771/93‑1 (MQ=255)
                               gAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGCGCGACAACTGGGACatgatg                                                          >  2:3273448/1‑93 (MQ=255)
                                aaaaaTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGt                                                         >  2:4544615/1‑93 (MQ=255)
                                  aaaTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGGCAACTGGGACATGATGTAc                                                       <  2:6148410/93‑1 (MQ=255)
                                  aaaTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTAc                                                       <  2:1349078/93‑1 (MQ=255)
                                    aTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACgg                                                     >  1:4927647/1‑92 (MQ=255)
                                    aTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGt                                                     >  2:37475/1‑93 (MQ=255)
                                    aTTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGt                                                     <  1:1367829/93‑1 (MQ=255)
                                     tttGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTg                                                    <  1:2160640/93‑1 (MQ=255)
                                      ttGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTgg                                                   >  1:4398749/1‑93 (MQ=255)
                                       tGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTggg                                                  <  2:2591273/93‑1 (MQ=255)
                                         ccTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTa                                                >  2:456297/1‑93 (MQ=255)
                                          cTACATTATCGACAACAACGGTCACATGCTGCGTATCATCGAACACCACACGATCTCCATGGGCGATAACTGGAGTATGATGTACGTGGGTAt                                               >  2:2873502/1‑93 (MQ=255)
                                            aCAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGt                                             >  2:865884/1‑93 (MQ=255)
                                              aaTATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTAc                                           >  1:3147089/1‑93 (MQ=255)
                                               atatCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTAcc                                          <  1:210148/93‑1 (MQ=255)
                                               atatCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTAcc                                          >  2:2176897/1‑93 (MQ=255)
                                                tatCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCa                                         >  1:2548039/1‑93 (MQ=255)
                                                 atCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAg                                        >  1:4243523/1‑93 (MQ=255)
                                                 atCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAg                                        <  2:1508633/93‑1 (MQ=255)
                                                  tCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAgg                                       >  2:4775311/1‑93 (MQ=255)
                                                   caacaacaaCGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGa                                      <  1:3677582/93‑1 (MQ=255)
                                                    aacaacaacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGat                                     <  1:1834133/93‑1 (MQ=255)
                                                       aacaacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATc                                  <  1:5742236/93‑1 (MQ=255)
                                                       aacaacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATc                                  <  2:1813524/93‑1 (MQ=255)
                                                       aacaacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATc                                  >  1:5607328/1‑93 (MQ=255)
                                                        acaacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCa                                 >  2:4010392/1‑93 (MQ=255)
                                                        acaacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCa                                 >  1:5406170/1‑93 (MQ=255)
                                                         caacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCaa                                >  1:3354913/1‑93 (MQ=255)
                                                          aacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAAc                               >  2:3201242/1‑93 (MQ=255)
                                                          aacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAAc                               <  2:4703614/93‑1 (MQ=255)
                                                          aacGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAAc                               >  2:2062745/1‑93 (MQ=255)
                                                           acGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACt                              >  2:2397909/1‑93 (MQ=255)
                                                           acGGTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACt                              <  2:4441261/93‑1 (MQ=255)
                                                              gTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTggg                           >  2:4993217/1‑93 (MQ=255)
                                                              gTCACATGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTggg                           <  1:4726910/93‑1 (MQ=255)
                                                                cacaTGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGAt                         <  1:4051562/93‑1 (MQ=255)
                                                                 acaTGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATa                        >  1:844268/1‑93 (MQ=255)
                                                                 acaTGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATa                        >  2:4467221/1‑93 (MQ=255)
                                                                  caTGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATaa                       <  2:4875731/93‑1 (MQ=255)
                                                                  caTGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATaa                       <  1:5029293/93‑1 (MQ=255)
                                                                  caTGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATaa                       <  2:5434050/93‑1 (MQ=255)
                                                                  caTGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATaa                       <  2:1137076/93‑1 (MQ=255)
                                                                   aTGCTGCGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAAc                      <  1:5656938/93‑1 (MQ=255)
                                                                    tgctgcGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACg                     >  1:5844309/1‑93 (MQ=255)
                                                                      ctgcGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGAc                   <  1:2691628/93‑1 (MQ=255)
                                                                       tgcGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACa                  <  2:1022948/93‑1 (MQ=255)
                                                                       tgcGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACa                  <  1:481791/93‑1 (MQ=255)
                                                                        gcGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACaa                 >  1:2040658/1‑93 (MQ=255)
                                                                         cGTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAAc                >  2:368461/1‑93 (MQ=255)
                                                                          gTATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACg               <  2:4937857/93‑1 (MQ=255)
                                                                           tATCCTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACgg              <  2:4330753/93‑1 (MQ=255)
                                                                              ccTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCAc           <  2:5097851/93‑1 (MQ=255)
                                                                              ccTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCAc           >  1:1619238/1‑93 (MQ=255)
                                                                              ccTCGACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCAc           >  2:6015257/1‑93 (MQ=255)
                                                                                  gACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAg       >  2:3492921/1‑93 (MQ=255)
                                                                                  gACCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAg       >  1:3303933/1‑93 (MQ=255)
                                                                                   aCCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAgt      >  2:3926687/1‑93 (MQ=255)
                                                                                   aCCACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAgt      <  2:5857287/93‑1 (MQ=255)
                                                                                    ccACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAgtg     <  2:1932314/93‑1 (MQ=255)
                                                                                     cACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAgtgg    >  2:5168774/1‑93 (MQ=255)
                                                                                     cACGGCGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAgtgg    >  1:5584037/1‑93 (MQ=255)
                                                                                          |                                                                                         
GTCGCAGGGTTCTGGCGTTGCATTTGATAACGAAAAATTTGCCTACAATATCAACAACAACGGTCACATGCTGCGTATCCTCGACCACGGTGCGATCTCCATGGGCGACAACTGGGACATGATGTACGTGGGTATGTACCAGGATATCAACTGGGATAACGACAACGGCACCAAGTGGTG  >  NC_000913/4246362‑4246541

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