Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,408 A→T 15.4% intergenic (+35/+459) icd → / ← ymfD isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,4080AT15.4% 310.1 / 46.5 123intergenic (+35/+459)icd/ymfDisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (44/60);  new base T (10/9);  total (54/69)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

Only 100 of 130 total aligned reads displayed.

TGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGG  >  NC_000913/1196309‑1196506
                                                                                                   |                                                                                                  
tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTa                                                                                                    >  1:2748642/1‑100 (MQ=255)
tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTa                                                                                                    <  2:472591/100‑1 (MQ=255)
 gatggatgGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTaa                                                                                                   >  2:1490891/1‑100 (MQ=255)
 gatggatgGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTaa                                                                                                   >  1:920490/1‑100 (MQ=255)
 gatggatgGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTaa                                                                                                   >  2:2025961/1‑100 (MQ=255)
  atggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAAc                                                                                                  <  2:1809205/100‑1 (MQ=255)
   tggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACg                                                                                                 <  2:1400914/100‑1 (MQ=255)
     gatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCt                                                                                               >  1:482837/1‑100 (MQ=255)
       tggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGATTGTTAAATTTATTAACGGGAGCGTAACGCTcc                                                                                             <  1:1928758/100‑1 (MQ=255)
       tagcgcTAAACTGCTGAACTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTACCTTTATTAACGGTCGCGTAACGCt                                                                                               >  1:1671058/3‑98 (MQ=21)
         gcgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCg                                                                                           <  2:657813/100‑1 (MQ=255)
          cgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgt                                                                                          <  1:1273061/100‑1 (MQ=255)
           gcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgtt                                                                                         <  2:2260037/100‑1 (MQ=255)
            cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgtta                                                                                        <  1:499375/100‑2 (MQ=255)
            cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgttg                                                                                        <  2:1720717/100‑1 (MQ=255)
            cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgaga                                                                                        >  1:423976/1‑97 (MQ=255)
              taaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaat                                                                                      <  2:2544052/99‑4 (MQ=255)
              aaaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaat                                                                                      <  2:2072087/100‑4 (MQ=255)
              aaaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgttgtt                                                                                      <  2:2355304/100‑1 (MQ=255)
               aaCTGCTGAAATGTACAGAGCTTGCTGACGCGATCATCGAAAACGTGTAATGCCGTAGTTTGTTAAATTGATTAACGGGAGAGTAACGCTCCGgttg                                                                                        >  1:776477/1‑97 (MQ=25)
                atcgctgAAATGTTCAGAGGTTGTTGACGCCAGCATCGAAAACATGGAATGCCGTACTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaat                                                                                      <  1:192263/95‑4 (MQ=25)
                aCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAACTTTATTAACGGGAGCGTAACGCTCCCGTTGtttt                                                                                    >  1:778327/1‑100 (MQ=255)
                  tgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTACCGCTCCCGTTGtttttt                                                                                  <  1:61815/100‑1 (MQ=255)
                   gctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTg                                                                                 <  1:1680673/100‑1 (MQ=255)
                    ctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTTACGCTCCCgttaataaattt                                                                                <  2:314894/100‑10 (MQ=25)
                     tgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTtttttttt                                                                               >  1:2435609/1‑100 (MQ=255)
                     tgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGtt                                                                               <  2:2266376/100‑1 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTa                                                                              >  2:1551578/1‑100 (MQ=255)
                      gAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTa                                                                              >  2:374602/1‑100 (MQ=255)
                        aaTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAgg                                                                            <  2:329532/100‑1 (MQ=255)
                         aTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATTTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgc                                                                           <  2:2049374/100‑1 (MQ=255)
                         aTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgc                                                                           <  2:2309136/100‑1 (MQ=255)
                           gTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctg                                                                         >  2:1091284/1‑100 (MQ=255)
                           gTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctg                                                                         <  2:738098/100‑1 (MQ=255)
                            ttCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctgc                                                                        <  1:1236680/100‑1 (MQ=255)
                             tCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGgctgct                                                                       <  2:2313110/100‑1 (MQ=255)
                               agagTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTATGCTGCTaa                                                                     >  1:2063329/1‑100 (MQ=255)
                               agagTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTaa                                                                     >  1:2399503/1‑100 (MQ=255)
                                  gtgtGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACgg                                                                  <  1:1354474/97‑1 (MQ=255)
                                     tGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTa                                                               >  2:481517/1‑100 (MQ=255)
                                     tGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTa                                                               <  1:642815/100‑1 (MQ=255)
                                      ggTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTAt                                                              >  1:1472909/1‑100 (MQ=255)
                                      ggTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTAt                                                              >  1:330806/1‑100 (MQ=255)
                                       gTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATc                                                             >  2:889329/1‑100 (MQ=255)
                                       gTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATc                                                             <  1:1008638/100‑1 (MQ=255)
                                       gTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATc                                                             >  2:1399108/1‑100 (MQ=255)
                                         gACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCaa                                                           >  2:2265346/1‑100 (MQ=255)
                                             cgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAtt                                                       <  1:1096362/100‑1 (MQ=255)
                                              gATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAttt                                                      <  2:1248652/100‑1 (MQ=255)
                                              gATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAttt                                                      >  2:64069/1‑100 (MQ=255)
                                                tcatcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTa                                                    <  1:1499916/100‑1 (MQ=255)
                                                tcatcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTa                                                    <  2:1294047/100‑1 (MQ=255)
                                                  atcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATc                                                  >  1:453678/1‑100 (MQ=255)
                                                  atcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATc                                                  >  2:382462/1‑100 (MQ=255)
                                                   tcGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCa                                                 <  1:136516/100‑1 (MQ=255)
                                                    cGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaa                                                <  2:1806254/100‑1 (MQ=255)
                                                     gAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCaaa                                               >  2:1410789/1‑100 (MQ=255)
                                                     gAAAACATGTAATGCCGTAGTTTGTTAAACTTATTAACGGGAGCGTAACGTTCCCGTTGTTTTTTGTTACGCTGCTAACGGTTATCAAAATTTTATCaaa                                               >  1:2609751/1‑100 (MQ=255)
                                                            tGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGtt                                        <  2:96146/100‑1 (MQ=255)
                                                             gTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTa                                       <  1:585849/100‑1 (MQ=255)
                                                               aaTGACGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAAAAAAAAAAAAAAAAGTTATc                                     <  1:2349142/100‑1 (MQ=255)
                                                                 tGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCTAAATTTTATCAAAAAAAGTTATCaa                                   <  1:2845906/100‑1 (MQ=255)
                                                                  gCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCaaa                                  >  1:1330680/1‑100 (MQ=255)
                                                                   ccGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCaaaa                                 >  1:2102019/1‑100 (MQ=255)
                                                                    cGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAc                                >  1:2802368/1‑100 (MQ=255)
                                                                    cGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAAc                                <  2:1043855/100‑1 (MQ=255)
                                                                        gTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCt                            >  2:1284772/1‑100 (MQ=255)
                                                                         tttGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTc                           <  2:1543211/100‑1 (MQ=255)
                                                                         tttGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTc                           <  1:1133209/100‑1 (MQ=255)
                                                                           tGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCgg                         <  2:1520228/100‑1 (MQ=255)
                                                                             ttAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGTCTGCTACCGGTTATCAAACTTTTATCAAAAAAAGTTATCAAAACCCCTCGGTa                       >  1:2155753/1‑100 (MQ=255)
                                                                              tAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAg                      >  2:880017/1‑100 (MQ=255)
                                                                                  tttATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGtttt                  <  2:641589/100‑1 (MQ=255)
                                                                                  aaaaacaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGtttt                  >  2:498769/7‑100 (MQ=255)
                                                                                   ttattaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTg                 <  2:725937/100‑1 (MQ=255)
                                                                                    tattaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgg                <  1:2779172/100‑1 (MQ=255)
                                                                                    tattaACGGGAGCGTAAAGCTCGCGTAGATTTTTGTAAGGCTGCTAACGGTTAACAAAAAATAAAAAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgg                <  1:727137/100‑1 (MQ=255)
                                                                                    tattaACGAGAGCGTAACGCTCCAATTATTATATCTAAGACTTCTTATGGTTATTAAAAAATTAAAAAAAAAAGTTATCAAAACCCCTCGGt                        <  1:1987780/92‑1 (MQ=255)
                                                                                     attaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTggg               <  1:1262624/100‑1 (MQ=255)
                                                                                      ttaACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgggg              <  1:1677975/100‑1 (MQ=255)
                                                                                      acaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAGACCCCTCGGTAGTTTTg                 >  1:340264/3‑97 (MQ=255)
                                                                                      acaaCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTgggg              >  2:552297/3‑100 (MQ=255)
                                                                                         aCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGTGGTAg           >  1:1401003/1‑100 (MQ=255)
                                                                                         aCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAg           >  2:36046/1‑100 (MQ=255)
                                                                                         aCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAg           >  2:2387518/1‑100 (MQ=255)
                                                                                         aCGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAg           >  2:2247416/1‑100 (MQ=255)
                                                                                          cGGGAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTGATCAAAACCCCTCTGTAGTTTTGGGGTAgg          >  2:132910/1‑100 (MQ=255)
                                                                                           gggAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGc         <  1:2515845/100‑1 (MQ=255)
                                                                                           gggAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGGTATCAAACCCCCTCGGTAGTTTTGTGGTAg           >  1:2122994/1‑98 (MQ=255)
                                                                                            ggAGCGTTACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCt        >  2:2317379/1‑100 (MQ=255)
                                                                                            ggAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCt        >  1:223393/1‑100 (MQ=255)
                                                                                             gAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTg       <  2:334824/100‑1 (MQ=255)
                                                                                              aGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTgg      <  1:161763/100‑1 (MQ=255)
                                                                                               gCGTAACGCTCCCGTTGTTTTTTATTAGCTTCCTCAGGGGTACGAAAGTATTATCACCATATGTTGCCACAACCGCTCGATAGCGTTGGGGttgtctagc     >  1:464706/1‑92 (MQ=255)
                                                                                                cGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGcc    <  2:555752/100‑1 (MQ=255)
                                                                                                cGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGcc    <  2:319043/100‑1 (MQ=255)
                                                                                                cGTAACGCTCCCGTTGTTTTTTGATAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGc         >  1:1082472/1‑95 (MQ=255)
                                                                                                  tAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTCTCAAAATTTTATCAAAAAAAGTTCTCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCgg  <  1:2434641/100‑1 (MQ=255)
                                                                                                  tAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCgg  <  1:2190768/100‑1 (MQ=255)
                                                                                                   aaCGCGCCCGTTGTTTTTTGTTAGGTTGCTAACGGTTATCAATATTTTATCAAAAAAAGTTATCAAAACTCCTCGGTAGTTTTGGGGTAGGCTGGcc    <  1:2610180/97‑1 (MQ=255)
                                                                                                   |                                                                                                  
TGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTTTTTTGTTAGGCTGCTAACGGTTATCAAAATTTTATCAAAAAAAGTTATCAAAACCCCTCGGTAGTTTTGGGGTAGGCTGGCCGG  >  NC_000913/1196309‑1196506

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