Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,765,829 (C)5→4 intergenic (‑53/‑89) yfjM ← / → rnlA CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, RNase LS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,765,8250C.98.6% 284.3 / ‑2.8 71intergenic (‑49/‑93)yfjM/rnlACP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, RNase LS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (34/36);  minor base G (0/1);  total (34/37)
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

TGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCA  >  NC_000913/2765689‑2765959
                                                                                                                                        |                                                                                                                                      
tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCGGCCCTATAGTTGGAGTGccc                                                                                                                                      >  2:1111063/1‑139 (MQ=255)
        ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAg                                                                                                                             <  2:935115/139‑1 (MQ=255)
              taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCt                                                                                                                       <  2:1526730/139‑1 (MQ=255)
                         tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCaa                                                                                                            <  2:283996/139‑1 (MQ=255)
                         tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCaa                                                                                                            >  1:672706/1‑139 (MQ=255)
                           tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                          >  1:79181/1‑139 (MQ=255)
                           tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                          >  2:39274/1‑139 (MQ=255)
                              aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAGATCaa                                                                                                       >  1:224585/1‑139 (MQ=255)
                              aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                                                                                                       <  1:895907/139‑1 (MQ=255)
                               gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAAt                                                                                                      >  1:804786/1‑139 (MQ=255)
                                ttGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                     <  1:1111063/139‑1 (MQ=255)
                                ttGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                     <  2:79181/139‑1 (MQ=255)
                                  gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                                                                   <  2:385903/139‑1 (MQ=255)
                                   ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAGATCAATacac                                                                                                  <  2:224585/139‑1 (MQ=255)
                                   ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                                                                                                  >  1:1272985/1‑139 (MQ=255)
                                   ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                     <  1:753231/136‑1 (MQ=255)
                                   ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                     >  2:753231/1‑136 (MQ=255)
                                        tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                                             <  2:577025/139‑1 (MQ=255)
                                          gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                           <  2:672706/139‑1 (MQ=255)
                                           tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                                                                          <  1:48460/139‑1 (MQ=255)
                                               aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                          <  1:788762/119‑1 (MQ=255)
                                               aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                          >  2:788762/1‑119 (MQ=255)
                                                    gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt                                                                                 >  2:1033998/1‑139 (MQ=255)
                                                    gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt                                                                                 >  1:1077503/1‑139 (MQ=255)
                                                         ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                     <  1:676303/114‑1 (MQ=255)
                                                         ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                     >  2:676303/1‑114 (MQ=255)
                                                         ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTaga                                                                            >  1:32053/1‑139 (MQ=255)
                                                         ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTaga                                                                            >  1:1601591/1‑139 (MQ=255)
                                                          gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                                                                                                      <  1:519643/81‑1 (MQ=255)
                                                          gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                                                                                                      >  2:519643/1‑81 (MQ=255)
                                                          gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTagag                                                                           >  1:367478/1‑139 (MQ=255)
                                                          gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTagag                                                                           >  1:668799/1‑139 (MQ=255)
                                                          gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTagag                                                                           >  1:248935/1‑139 (MQ=255)
                                                          gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTagag                                                                           >  2:352596/1‑139 (MQ=255)
                                                           ttACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                                                          <  2:248935/139‑1 (MQ=255)
                                                               cGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCa                                                                                                                              >  2:240646/1‑83 (MQ=255)
                                                               cGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCa                                                                                                                              <  1:240646/83‑1 (MQ=255)
                                                               cGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTCTcc                                                                      <  2:367478/139‑1 (MQ=255)
                                                                  aaCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCa                                                                   <  1:1033998/139‑1 (MQ=255)
                                                                            atcatcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCCGCCCTATAGTTGGAG‑GGCCCCTATTCAGTTGGCTAAATTTAGCAATTCAATACACTTCAGGGGGGTATTAATGTAGAGTTTCCCCATATGTTTCTa                                                         <  2:1621265/139‑1 (MQ=255)
                                                                              catcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg                                                       <  1:39274/139‑1 (MQ=255)
                                                                                            tgGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGCGGGGTATTAATGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATg                                         <  2:1019749/139‑1 (MQ=255)
                                                                                                cTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                     <  2:1319433/75‑1 (MQ=255)
                                                                                                cTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                     >  1:1319433/1‑75 (MQ=255)
                                                                                                  gCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                                                                                                      <  1:149959/41‑1 (MQ=255)
                                                                                                  gCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                                                                                                      >  2:149959/1‑41 (MQ=255)
                                                                                                  gCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                                                                   >  2:1366803/1‑75 (MQ=255)
                                                                                                  gCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                                                                   <  1:1366803/75‑1 (MQ=255)
                                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAg                                                                                                                             >  2:346046/1‑48 (MQ=255)
                                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAg                                                                                                                             <  1:346046/48‑1 (MQ=255)
                                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg                                                                                                               >  1:1625173/1‑62 (MQ=255)
                                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg                                                                                                               <  2:1625173/62‑1 (MQ=255)
                                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCaa                                                                                                            <  1:717323/65‑1 (MQ=255)
                                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCaa                                                                                                            >  2:717323/1‑65 (MQ=255)
                                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCa                                  >  2:864038/1‑139 (MQ=255)
                                                                                                                    ccTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                        <  1:1242706/68‑1 (MQ=255)
                                                                                                                    ccTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                        >  2:1242706/1‑68 (MQ=255)
                                                                                                                    ccTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGt                                                              >  1:1175554/1‑94 (MQ=255)
                                                                                                                    ccTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGt                                                              <  2:1175554/94‑1 (MQ=255)
                                                                                                                     cTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATg                                         >  1:1025909/1‑114 (MQ=255)
                                                                                                                     cTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATg                                         <  2:1025909/114‑1 (MQ=255)
                                                                                                                      tGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                     >  1:1657001/1‑53 (MQ=255)
                                                                                                                      tGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                     <  2:1657001/53‑1 (MQ=255)
                                                                                                                             taGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                                      <  2:308244/61‑1 (MQ=255)
                                                                                                                             taGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                                      >  1:308244/1‑61 (MQ=255)
                                                                                                                             taGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTc        >  1:1278262/1‑139 (MQ=255)
                                                                                                                             taGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTc        <  2:1272985/139‑1 (MQ=255)
                                                                                                                               gTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAg      <  2:1601591/139‑1 (MQ=255)
                                                                                                                                  ggAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                           >  1:1059130/1‑51 (MQ=255)
                                                                                                                                  ggAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                           <  2:1059130/51‑1 (MQ=255)
                                                                                                                                  ggAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                                             >  1:136986/1‑49 (MQ=255)
                                                                                                                                  ggAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                                             <  2:136986/49‑1 (MQ=255)
                                                                                                                                  ggAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAAt                                    <  1:128919/106‑1 (MQ=255)
                                                                                                                                  ggAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAAt                                    >  2:128919/1‑106 (MQ=255)
                                                                                                                                  ggAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGc   >  1:459924/1‑139 (MQ=255)
                                                                                                                                   gAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCa  <  2:459924/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGCA  >  NC_000913/2765689‑2765959

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