Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,874,663 (G)5→6 intergenic (‑83/+32) yidA ← / ← b3698 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,874,6581.G98.7% 939.7 / 3.2 305intergenic (‑78/+37)yidA/b3698hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base G (153/148);  total (157/148)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 366 total aligned reads displayed.

GAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGT  >  NC_000913/3874566‑3874747
                                                                                               |                                                                                          
gAGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGaaaa                                                                                             <  2:5555033/93‑1 (MQ=255)
 aGTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAg                                                                                           <  1:1311133/93‑1 (MQ=255)
  gTTTAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgg                                                                                          <  1:2950750/93‑1 (MQ=255)
    ttAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgggg                                                                                       >  2:6060430/1‑93 (MQ=255)
     tAATAGCCGTGAGTTCCCCGGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACACCACCCGAGGGAAAAggggg                                                                                      <  1:7008412/93‑1 (MQ=255)
     tAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggggg                                                                                      >  1:7549026/1‑93 (MQ=255)
     tAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggggg                                                                                      <  1:5548125/93‑1 (MQ=255)
     tAATAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggggg                                                                                      <  2:5321777/93‑1 (MQ=255)
      aaTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgggggg                                                                                     <  2:2305376/93‑2 (MQ=255)
       aTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggggggg                                                                                    >  2:5488629/1‑91 (MQ=255)
       aTAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAgggggga                                                                                    <  2:365747/93‑3 (MQ=255)
        tAGCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAggggggag                                                                                   >  1:3997037/1‑90 (MQ=255)
          gCCGTGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGTGGGATGc                                                                                  >  1:7175344/1‑93 (MQ=255)
           ccgtGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCa                                                                                 <  1:7651704/90‑1 (MQ=255)
            cgtGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCaa                                                                                >  2:6883381/3‑93 (MQ=255)
             gtGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAAc                                                                               >  1:2145983/2‑93 (MQ=255)
             gtGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAAc                                                                               <  1:1979538/92‑1 (MQ=255)
              tGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAAcg                                                                              >  2:3436721/1‑93 (MQ=255)
               gAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAAcgc                                                                             <  1:4598247/93‑1 (MQ=255)
                aGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAAcgcg                                                                            >  2:5010852/1‑93 (MQ=255)
                 gTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAAcgcgc                                                                           >  1:4756379/1‑93 (MQ=255)
                  ttCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAAcgcgcg                                                                          >  1:3426463/1‑93 (MQ=255)
                   tCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGt                                                                         <  2:6985300/93‑1 (MQ=255)
                     cccAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGGAAAGGGGGGATGCAACGCGCGTTg                                                                       >  1:1517406/1‑93 (MQ=255)
                      ccAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGTAAAGGGGGGGTGCAACGCGCGTTGc                                                                      >  1:6318511/1‑93 (MQ=255)
                      ccAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGc                                                                      >  2:616599/1‑93 (MQ=255)
                       cAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCg                                                                     >  2:7530806/1‑93 (MQ=255)
                       cAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCg                                                                     >  2:1040396/1‑93 (MQ=255)
                        gggTAGAGTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGt                                                                    >  2:4797468/2‑93 (MQ=255)
                        aGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGt                                                                    <  1:5055092/93‑1 (MQ=255)
                         ggTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTg                                                                   >  1:5313301/1‑93 (MQ=255)
                          gTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGc                                                                  >  1:174746/1‑93 (MQ=255)
                          gTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGc                                                                  >  2:1349546/1‑93 (MQ=255)
                           tAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCt                                                                 <  2:3284335/93‑1 (MQ=255)
                           tAGATTGGTTACTGACTAATGCTAACGCGATTCACCTCAAAAATCAGTAAAACACCCGAGGGAAAAGGGGGGATGCAACGAGCGCTGCGTgat                                                                 >  2:5214313/1‑91 (MQ=255)
                             gATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTcc                                                               >  1:4980117/1‑93 (MQ=255)
                              aTTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTccc                                                              <  1:6591263/93‑1 (MQ=255)
                              aTTGGTTCCTGACCCATGCTAACGCGATTCCGCTAAAAAATCAGTACAACCCCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTccc                                                              <  2:3036515/93‑1 (MQ=255)
                               ttGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCt                                                             <  2:5370618/93‑1 (MQ=255)
                                 ggTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCttt                                                           <  2:152454/93‑1 (MQ=255)
                                 ggTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCttt                                                           <  1:2826540/93‑1 (MQ=255)
                                  gTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACCCCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCtttt                                                          <  2:3513441/93‑1 (MQ=255)
                                  gTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCtttt                                                          <  1:5211414/93‑1 (MQ=255)
                                   ttCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCttttt                                                         <  1:6845381/93‑1 (MQ=255)
                                     ccTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGc                                                       <  2:6729650/93‑1 (MQ=255)
                                     ccTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGc                                                       <  1:3833243/93‑1 (MQ=255)
                                       tGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCtt                                                     <  1:1353500/93‑1 (MQ=255)
                                        gccccATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTa                                                    <  2:562392/91‑1 (MQ=255)
                                          cccATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGc                                                  <  1:5646589/93‑1 (MQ=255)
                                           ccATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCg                                                 >  1:2179715/1‑93 (MQ=255)
                                            cATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCgg                                                >  1:4208226/1‑93 (MQ=255)
                                            cATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCgg                                                >  2:3067008/1‑93 (MQ=255)
                                              tGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGaa                                              >  1:4424485/1‑93 (MQ=255)
                                              tGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGaa                                              <  1:2538320/93‑1 (MQ=255)
                                                cTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAaga                                            <  2:2596429/93‑1 (MQ=255)
                                                   aCGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGtt                                         >  2:2532794/1‑93 (MQ=255)
                                                    cgcgATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGttt                                        <  1:4059001/93‑1 (MQ=255)
                                                    cgcgATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGttt                                        <  2:3825230/93‑1 (MQ=255)
                                                     gcgATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTc                                       >  1:579353/1‑93 (MQ=255)
                                                     gcgATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTc                                       <  1:2427954/93‑1 (MQ=255)
                                                       gATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTccc                                     >  1:5430262/1‑93 (MQ=255)
                                                       gATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTccc                                     <  2:4604581/93‑1 (MQ=255)
                                                        aTTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCt                                    <  2:4731221/93‑1 (MQ=255)
                                                        aTTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCt                                    <  1:7936495/93‑1 (MQ=255)
                                                        aTTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCt                                    >  1:7702434/1‑93 (MQ=255)
                                                        aTTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCt                                    >  2:7049525/1‑93 (MQ=255)
                                                         ttCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCtt                                   >  1:274561/1‑93 (MQ=255)
                                                           ccGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTc                                 >  2:1698562/1‑93 (MQ=255)
                                                             gCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTcag                               >  2:2911895/1‑93 (MQ=255)
                                                             gCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTcag                               >  2:1129906/1‑93 (MQ=255)
                                                               tCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTcagca                             >  2:7778165/1‑93 (MQ=255)
                                                               tCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTcagca                             <  1:11761/93‑1 (MQ=255)
                                                                cAAAAATCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTcagcag                            <  1:6405813/93‑1 (MQ=255)
                                                                 aaaaaTCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGt                           >  2:6583809/1‑93 (MQ=255)
                                                                 aaaaaTCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGt                           >  2:81511/1‑93 (MQ=255)
                                                                  aaaaTCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGtt                          <  2:2049871/93‑1 (MQ=255)
                                                                  aaaaTCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGtt                          >  1:3303849/1‑93 (MQ=255)
                                                                   aaaTCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTc                         >  2:4033237/1‑93 (MQ=255)
                                                                    aaTCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGGGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTcc                        <  1:1821599/93‑1 (MQ=255)
                                                                     aTCAGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCa                       >  1:1007875/1‑93 (MQ=255)
                                                                        aGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGc                    <  1:958372/93‑1 (MQ=255)
                                                                        aGTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGc                    >  1:4691199/1‑93 (MQ=255)
                                                                         gTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGcc                   <  2:2504374/93‑1 (MQ=255)
                                                                         gTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGcc                   >  1:1858970/1‑93 (MQ=255)
                                                                         gTACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGcc                   >  2:7484607/1‑93 (MQ=255)
                                                                          tACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCt                  <  2:491658/93‑1 (MQ=255)
                                                                          tACAACACCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCt                  <  1:36596/93‑1 (MQ=255)
                                                                             aacaCCCGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCg               >  1:4878244/1‑93 (MQ=255)
                                                                                 cccGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGaaa           >  2:6554775/1‑93 (MQ=255)
                                                                                  ccGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAg          >  1:621482/1‑93 (MQ=255)
                                                                                   cGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGc         <  2:5919513/93‑1 (MQ=255)
                                                                                   cGAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGc         >  2:7420962/1‑93 (MQ=255)
                                                                                    gAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCa        >  1:6694147/1‑93 (MQ=255)
                                                                                    gAGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCa        >  2:218764/1‑93 (MQ=255)
                                                                                     aGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAg       >  1:4362171/1‑93 (MQ=255)
                                                                                     aGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAg       <  2:1458157/93‑1 (MQ=255)
                                                                                     aGGGAAAAGGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAg       >  2:7733832/1‑93 (MQ=255)
                                                                                        ggaaaagGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGAt     <  1:3539655/90‑1 (MQ=255)
                                                                                        ggaaaagGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGAt     <  2:4265922/90‑1 (MQ=255)
                                                                                         gaaaagGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATc    <  2:6404144/91‑1 (MQ=255)
                                                                                               |                                                                                          
GAGTTTAATAGCCATGAGTTCCCCAGGTAGATTGGTTCCTGACCCATGCTAACGCGATTCCGCTCAAAAATCAGTACAACACCCGAGGGAAAAGGGGGATGCAACGCGCGTTGCGTGCTCCCTTTTTGCTTAGCGGAAGAGTTTCCCTTTCAGCAGTTCCATGCCTGCGGAAAGCAGATCGT  >  NC_000913/3874566‑3874747

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