Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,525,135 (A)6→5 100% coding (9/765 nt) yfeN → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,525,1300A.69.7% 155.5 / 136.0 99coding (4/765 nt)yfeNputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (13/17);  new base . (36/33);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

GGCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTG  >  NC_000913/2525031‑2525225
                                                                                                   |                                                                                               
ggCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAATCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGa                                                                                                 <  2:288153/100‑1 (MQ=255)
ggCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGa                                                                                                 <  1:2504952/100‑1 (MQ=255)
ggCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGa                                                                                                 >  2:2178757/1‑100 (MQ=255)
 gCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaa                                                                                                >  2:1863114/1‑100 (MQ=255)
  cGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaa                                                                                               >  2:110825/1‑100 (MQ=255)
   gATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaaa                                                                                              <  2:285007/100‑1 (MQ=255)
     tatTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaaaac                                                                                            >  2:1734690/1‑99 (MQ=255)
     tatTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaaaac                                                                                            <  1:1511211/100‑2 (MQ=255)
     tatTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaaaaa                                                                                            >  2:791523/1‑100 (MQ=255)
      atTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGaaaaaca                                                                                           >  2:2770506/1‑98 (MQ=255)
       ttAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACa                                                                                          <  1:2440129/100‑1 (MQ=255)
         aaTTCCTTTTCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATc                                                                                        <  2:517868/100‑1 (MQ=255)
         aaTTCCTTTGCATTAACAATATGTCGTTTTTAACAAACCACTTTATTAAACAATCTACCAATAACAATCGCCATAGGGAATTATGATATACAAAACATCt                                                                                       >  1:3201399/1‑100 (MQ=255)
           ttCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCtt                                                                                      >  2:616653/1‑100 (MQ=255)
           ttCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCttt                                                                                     <  2:657963/100‑1 (MQ=255)
              cTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTa                                                                                   >  2:1219078/1‑100 (MQ=255)
               tttGCAATAACAATATGTCGTTTTTAGGAAACCACTTTATTAAGCAACCTACGAATAACGATCGGCCTAGGGAATTATGCTATGAAAAACATCTTTTAAc                                                                                 <  1:2240037/100‑1 (MQ=255)
                 tGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAActc                                                                               >  2:579679/1‑100 (MQ=255)
                   cATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAActct                                                                              >  1:544508/1‑100 (MQ=255)
                    aTTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTg                                                                             <  1:819267/100‑1 (MQ=255)
                     ttAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGac                                                                           <  1:877046/100‑1 (MQ=255)
                     ttAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGac                                                                           >  1:1807138/1‑100 (MQ=255)
                       aaCAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGaca                                                                          <  1:2789999/100‑1 (MQ=255)
                       aaCAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGacac                                                                         >  1:1580424/1‑100 (MQ=255)
                          aaTATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACtt                                                                       >  2:3273049/1‑100 (MQ=255)
                          aaTATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACttt                                                                      >  1:2594220/1‑100 (MQ=255)
                           atatGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTc                                                                     >  1:2523990/1‑100 (MQ=255)
                            tatGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTcc                                                                    <  2:2942306/100‑1 (MQ=255)
                              tGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCctc                                                                  <  1:1690987/100‑1 (MQ=255)
                                  gTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGGAAAACATCTTTTAACTCTGATACTTTCCTCtata                                                              <  2:2299101/100‑1 (MQ=255)
                                  gTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCtata                                                              >  2:2744087/1‑100 (MQ=255)
                                  gTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGAGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCctct                                                                 <  1:1337299/98‑1 (MQ=255)
                                   tttttAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTatat                                                             <  2:1943303/100‑1 (MQ=255)
                                    ttttAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTatat                                                             >  1:2299434/1‑100 (MQ=255)
                                    ttttAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTatat                                                             <  1:1737829/100‑1 (MQ=255)
                                    ttttAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATAtt                                                            >  1:504338/1‑100 (MQ=255)
                                     tttAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATAtt                                                            >  2:3029122/1‑99 (MQ=255)
                                     tttAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTa                                                           >  2:19745/1‑100 (MQ=255)
                                      ttAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAg                                                          >  2:454344/1‑100 (MQ=255)
                                       tAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGc                                                         <  2:2756169/100‑1 (MQ=255)
                                        aGCAAACCACTTTATTAAGCAATCTACTAATATCGATCGGCATAGGGAATTATGCAATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCg                                                        <  1:625501/100‑1 (MQ=255)
                                        aGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCg                                                        <  1:1152787/100‑1 (MQ=255)
                                         gCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGa                                                       >  1:1700851/1‑100 (MQ=255)
                                          cAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGCAAAACATCTTTTAACTCTGACACTTTCTTCTATATTAGCGAt                                                      >  1:1217173/1‑100 (MQ=255)
                                          cAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGAt                                                      <  1:763552/100‑1 (MQ=255)
                                           aaaCCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATa                                                     >  1:1612779/1‑100 (MQ=255)
                                           aaaCCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATa                                                     <  2:2702515/100‑1 (MQ=255)
                                            aaCCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATa                                                     >  1:107529/1‑100 (MQ=255)
                                              ccACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCt                                                  >  1:2270264/1‑100 (MQ=255)
                                                aCTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtgt                                                <  2:1829936/100‑1 (MQ=255)
                                                aCTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtgt                                                >  2:809971/1‑100 (MQ=255)
                                                aCTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtgt                                                <  1:1734679/100‑1 (MQ=255)
                                                  tttATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCtgttg                                              <  1:2923590/100‑1 (MQ=255)
                                                    tattaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTa                                            >  1:107508/1‑100 (MQ=255)
                                                    tattaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTa                                            <  1:3254987/100‑1 (MQ=255)
                                                      ttaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATAGCGGTTGTagc                                          <  1:868163/100‑3 (MQ=255)
                                                      ttaAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTAtc                                          >  2:2962842/1‑100 (MQ=255)
                                                       taAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTAtct                                         >  1:1455164/1‑100 (MQ=255)
                                                        aaGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTAtctc                                        >  2:2947781/1‑100 (MQ=255)
                                                         aGCAATCTACTAATAACGATCGGCATAGTGAATTATGCTATGAAAATCATCTTTTATCTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCa                                       <  1:3282876/100‑1 (MQ=255)
                                                         aGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATTCCTGTTGTAtcac                                        <  2:2300210/100‑3 (MQ=255)
                                                          gCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTAtctc                                        <  1:877073/99‑1 (MQ=255)
                                                          gCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCa                                       <  1:1063694/100‑1 (MQ=255)
                                                           cAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCa                                       <  1:1962547/99‑1 (MQ=255)
                                                           cAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATg                                     >  2:183418/1‑100 (MQ=255)
                                                            aaTCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGc                                    <  2:1197466/100‑1 (MQ=255)
                                                            aaTCTACACAGAACGATCTGGGTAGGGAATTATCCTATGAAAAAACATTTTTTAACTCTGACACTTTCCGCTATATTAGCTATACCTGTTGTATCTCa                                       <  1:3176496/98‑1 (MQ=255)
                                                             aTCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGc                                    <  1:2101787/100‑1 (MQ=255)
                                                             aTCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCa                                   <  2:2699602/100‑1 (MQ=255)
                                                                tactaATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGaa                                >  2:2672318/1‑100 (MQ=255)
                                                                  ctaATAACGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAAtt                              >  1:3244291/1‑100 (MQ=255)
                                                                   taataaCGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAAttt                             <  2:554087/100‑1 (MQ=255)
                                                                    aataaCGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAAttt                             >  1:576326/1‑100 (MQ=255)
                                                                    aataaCGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAAttt                             >  1:1580047/1‑100 (MQ=255)
                                                                    aataaCGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATAGCTGTTGTATCTCATGCAGAATTTa                            <  1:2672506/100‑1 (MQ=255)
                                                                      taaCGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTaaa                          >  1:2977700/1‑100 (MQ=255)
                                                                       aaCGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAg                         <  1:840261/100‑1 (MQ=255)
                                                                       aaCGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAg                         <  1:1642145/100‑1 (MQ=255)
                                                                        aCGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAgg                        >  1:1789414/1‑100 (MQ=255)
                                                                         cGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAgg                        <  2:2970302/100‑1 (MQ=255)
                                                                         cGATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGc                       <  2:2792397/100‑1 (MQ=255)
                                                                          naTCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCg                      >  2:443454/2‑100 (MQ=255)
                                                                          gATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGc                       >  2:2924992/1‑100 (MQ=255)
                                                                          gATCGGCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCg                      >  1:2484498/1‑100 (MQ=255)
                                                                             cGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGt                    <  2:767815/100‑1 (MQ=255)
                                                                              ggCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGttt                  >  1:2450706/1‑100 (MQ=255)
                                                                               gCATAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGtttt                 >  1:3174496/1‑100 (MQ=255)
                                                                                 aTAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGc               <  1:691248/100‑1 (MQ=255)
                                                                                  tAGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCa              >  1:1795612/1‑100 (MQ=255)
                                                                                   aGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAg             <  1:419408/100‑1 (MQ=255)
                                                                                   aGGGAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAACGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAg             >  1:1661344/1‑100 (MQ=255)
                                                                                    gggAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAg             <  2:3276713/100‑1 (MQ=255)
                                                                                    gggAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAg             <  2:1919431/100‑1 (MQ=255)
                                                                                    gggAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGa            >  1:430104/1‑100 (MQ=255)
                                                                                    gggAATTATGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGa            <  2:1648357/100‑1 (MQ=255)
                                                                                     ggAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGa            <  1:110825/100‑1 (MQ=255)
                                                                                      gAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGat           >  2:1571033/1‑100 (MQ=255)
                                                                                           aTGCTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGc     >  2:407866/1‑100 (MQ=255)
                                                                                              cTATGAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTg  >  2:2887535/1‑100 (MQ=255)
                                                                                               tATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTg  >  2:3271270/1‑100 (MQ=255)
                                                                                                   |                                                                                               
GGCGATATTAATTCCTTTGCATTAACAATATGTCGTTTTTAGCAAACCACTTTATTAAGCAATCTACTAATAACGATCGGCATAGGGAATTATGCTATGAAAAAACATCTTTTAACTCTGACACTTTCCTCTATATTAGCGATACCTGTTGTATCTCATGCAGAATTTAAAGGCGGTTTTGCAGATATCGGCGTG  >  NC_000913/2525031‑2525225

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