Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,817,948 G→T 16.4% E63D (GAG→GAT dinD → DNA damage‑inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,817,9480GT16.4% 490.3 / 75.4 191E63D (GAG→GATdinDDNA damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base G (77/82);  new base T (15/17);  total (92/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

Only 100 of 214 total aligned reads displayed.

CTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGAC  >  NC_000913/3817850‑3818047
                                                                                                  |                                                                                                   
ctGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGgaga                                                                                                    >  2:1259962/1‑100 (MQ=255)
ctGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGgaga                                                                                                    >  1:3816996/1‑100 (MQ=255)
 tGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGAt                                                                                                   >  2:2401688/1‑100 (MQ=255)
   ggTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATcc                                                                                                 >  2:4018445/1‑100 (MQ=255)
   ggTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATcc                                                                                                 <  1:903820/100‑1 (MQ=255)
      tATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTTTGACCATATCGAGGATCTCCTCGa                                                                                             <  2:2570492/100‑1 (MQ=37)
         tCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGata                                                                                           >  2:1117381/1‑100 (MQ=255)
          cAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGATCTCCTCGaaat                                                                                          <  2:1493871/100‑1 (MQ=255)
          cAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGatat                                                                                          >  2:78895/1‑100 (MQ=255)
          cAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGatat                                                                                          <  2:1799261/100‑1 (MQ=255)
          cAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGatat                                                                                          >  2:3154015/1‑100 (MQ=255)
          cAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGatat                                                                                          >  1:4178266/1‑100 (MQ=255)
           agggTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATGGATGACCATTTCGAGGAGATCCTCGATATg                                                                                         >  1:200579/4‑100 (MQ=255)
           agagTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTTTGACCATATCGAGGATCTCCTCGAAATgg                                                                                        <  1:1215182/100‑1 (MQ=25)
             agTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGt                                                                                       <  1:2699705/100‑1 (MQ=255)
              gTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGTc                                                                                      <  2:1063739/100‑1 (MQ=38)
                 tCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGTATAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCaaa                                                                                   >  1:3161283/1‑100 (MQ=255)
                  cGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGTCatca                                                                                  <  2:4034566/100‑4 (MQ=37)
                   gTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAAt                                                                                 <  1:1601746/100‑1 (MQ=255)
                   gTCACTTTACACCTGTATTAACCCGCGCGTGAGAAGCTTGTGAGAACAGTGGCCACAGCATTGATGACCATATCGATGATCTCCTCGAAATGGTCat                                                                                    <  1:1804194/97‑2 (MQ=16)
                    tCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAAtt                                                                                <  1:1439419/100‑1 (MQ=255)
                     cACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTg                                                                               <  1:2731866/100‑1 (MQ=255)
                     cACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCatcaattg                                                                              <  2:3479220/100‑8 (MQ=18)
                       cTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGc                                                                             >  2:2878896/1‑100 (MQ=255)
                       cTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGc                                                                             <  1:2040858/100‑1 (MQ=255)
                       cTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGc                                                                             <  2:1838489/100‑1 (MQ=255)
                       cTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTTTGACCATATCGAGGATCTCCTCGAAATGGTCatcaattgtg                                                                            <  2:1021338/100‑10 (MQ=17)
                        tttATACCTGTATTAACGGGCGCCAAAGAAGCGTGTGAGAAGAGTCGTCACACAATTGATGTCGATTTCGAGGAGATGCTCGAAATGGTCAAAATTGGCt                                                                            <  1:248151/100‑1 (MQ=255)
                        tttATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCt                                                                            <  2:2176241/100‑1 (MQ=255)
                         ttATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTc                                                                           <  1:1818160/100‑1 (MQ=255)
                         ttATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTc                                                                           <  2:1933273/100‑1 (MQ=255)
                          tataCCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCa                                                                          >  1:2468676/1‑100 (MQ=255)
                            taCCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCaaa                                                                        >  2:2430512/1‑100 (MQ=255)
                             aCCTGTATTAACGCGCGCCAAAGGAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGACTCGCGATAGGGTCAAAATTGGCTCcaa                                                                        >  1:3599327/1‑99 (MQ=255)
                              ccTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGTTGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATg                                                                      <  2:1920660/100‑1 (MQ=255)
                                 gTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCa                                                                   <  2:3273566/100‑1 (MQ=255)
                                 gTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCa                                                                   <  1:1346523/100‑1 (MQ=255)
                                  tATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCaa                                                                  >  1:2087533/1‑100 (MQ=255)
                                    ttAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCaaaa                                                                >  1:216263/1‑100 (MQ=255)
                                    ttAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCaaaa                                                                >  2:3794768/1‑100 (MQ=255)
                                      aaCGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAaga                                                              >  1:1684876/1‑100 (MQ=255)
                                         gcgcgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCa                                                           <  1:4224844/100‑1 (MQ=255)
                                           gcgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGATATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCAtt                                                         <  1:103473/100‑1 (MQ=255)
                                           gcgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCAtt                                                         <  1:756726/100‑1 (MQ=255)
                                            cgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTa                                                        <  1:1019638/100‑1 (MQ=255)
                                                 aagaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGa                                                   <  2:2167909/100‑1 (MQ=255)
                                                  agaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGAc                                                  >  2:755511/1‑100 (MQ=255)
                                                           tgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTActc                                         >  1:1303263/1‑100 (MQ=255)
                                                           tgAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATAATCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTAct                                          <  1:10763/99‑1 (MQ=255)
                                                              gaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTcc                                      >  2:511298/1‑100 (MQ=255)
                                                                 cAGTAGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGc                                   <  1:454113/100‑1 (MQ=255)
                                                                   gTGGTCACACAATTGATGCCCATTTCGACGAGATCCTCGATATGGTCAAAATTGACTGAAATGCCAAAAGAGCATTAAACGCCAAACAACTCTCCCGCTa                                 >  1:2774125/1‑100 (MQ=255)
                                                                    tGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTAt                                <  2:2074183/100‑1 (MQ=255)
                                                                     gttCACACAATTGATGACCATTTCGAGGATATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATg                               <  2:3972478/98‑1 (MQ=255)
                                                                      gTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGc                              >  1:4240425/1‑100 (MQ=255)
                                                                      gTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGc                              >  1:3949360/1‑100 (MQ=255)
                                                                      gTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGc                              >  1:3036101/1‑100 (MQ=255)
                                                                       tCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                             <  1:957801/100‑1 (MQ=255)
                                                                       tCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                             <  1:4018418/100‑1 (MQ=255)
                                                                       tCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGcc                             <  1:353632/100‑1 (MQ=255)
                                                                       tCACACAATTGATGACCATTTCGAGGAGATCCGCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGAATTAAAAGACATCGTACTCTCCCGCTATGcc                             <  1:3480329/100‑1 (MQ=255)
                                                                         acacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTg                           <  2:596941/100‑1 (MQ=255)
                                                                          cacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGt                          <  2:3848240/100‑1 (MQ=255)
                                                                          cacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGt                          >  2:3004377/1‑100 (MQ=255)
                                                                          cacaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGt                          >  1:468847/1‑100 (MQ=255)
                                                                           acaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGtt                         >  2:1664884/1‑100 (MQ=255)
                                                                              gatcttgaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTAcc                      <  2:1045472/95‑1 (MQ=255)
                                                                              aTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTAcc                      <  2:3519796/100‑1 (MQ=255)
                                                                              aTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTAc                       >  1:2310907/1‑99 (MQ=255)
                                                                               ttGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCt                     <  2:1478996/100‑1 (MQ=255)
                                                                                tgatgaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTg                    <  1:4212781/100‑1 (MQ=255)
                                                                                 gatgaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                   >  2:3093737/1‑100 (MQ=255)
                                                                                 gatgaCCATTTCGAGGAGACCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAACAGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                   >  1:339991/1‑100 (MQ=255)
                                                                                  ttgaCCATAACGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATCAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgt                  >  1:251733/2‑100 (MQ=37)
                                                                                  atgaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgt                  <  1:1050874/100‑1 (MQ=255)
                                                                                   tgaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCATATGCCAAAACAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTCCCTGgta                 >  1:1628093/1‑100 (MQ=255)
                                                                                    gaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtag                >  2:3390803/1‑100 (MQ=255)
                                                                                    gaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtag                >  2:1175024/1‑100 (MQ=255)
                                                                                    gaCCATATCGAGGATCTCCTCGAAATGGTCAAAATTGGTTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtag                >  1:2707652/1‑100 (MQ=37)
                                                                                    gaCCATATCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCtagatc                >  2:1148665/1‑95 (MQ=37)
                                                                                    gaCCATATCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtag                >  1:4124485/1‑100 (MQ=38)
                                                                                     aCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagt               >  2:4185063/1‑100 (MQ=255)
                                                                                      ccATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagta              >  1:2247974/1‑100 (MQ=255)
                                                                                      ccATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagta              >  1:4338589/1‑100 (MQ=255)
                                                                                      ccATATCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagta              >  1:1152397/1‑100 (MQ=38)
                                                                                       cATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTAc             <  2:4266788/100‑1 (MQ=255)
                                                                                        aTTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACa            >  1:689371/1‑100 (MQ=255)
                                                                                        aTTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACa            <  1:4220042/100‑1 (MQ=255)
                                                                                         tatCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACaa           >  2:1788482/3‑100 (MQ=39)
                                                                                          atCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACaa           >  1:1507263/2‑99 (MQ=39)
                                                                                            cGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAAc        >  1:262250/1‑100 (MQ=255)
                                                                                              aggagATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACgg      >  2:3395675/1‑100 (MQ=255)
                                                                                              aggagATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACgg      <  2:2953020/100‑1 (MQ=255)
                                                                                               tctgaTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGc     <  2:4060541/97‑1 (MQ=255)
                                                                                               ggagATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGc     <  2:2817067/100‑1 (MQ=255)
                                                                                               ggagATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGc     >  1:2298216/1‑100 (MQ=255)
                                                                                                 agaTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGa   <  2:3334343/100‑1 (MQ=255)
                                                                                                  gaTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGAc  >  2:3687489/1‑100 (MQ=255)
                                                                                                  gaTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGAc  >  2:3557095/1‑100 (MQ=255)
                                                                                                  gaTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGAc  >  2:867746/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGAC  >  NC_000913/3817850‑3818047

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