Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,101,039 (C)7→9 coding (585/720 nt) yggN ← putative EcfF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,0321.C98.1% 177.9 / ‑3.2 53coding (592/720 nt)yggNputative EcfF
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (22/30);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01
*NC_0009133,101,0322.C96.2% 172.5 / 0.0 53coding (592/720 nt)yggNputative EcfF
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (22/29);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

ATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑CCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTT  >  NC_000913/3100896‑3101160
                                                                                                                                         ||                                                                                                                                    
aTAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑ccc                                                                                                                                  <  2:405877/140‑1 (MQ=255)
 tAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cccc                                                                                                                                 >  2:1284966/1‑140 (MQ=255)
 tAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cccc                                                                                                                                 >  2:408402/1‑140 (MQ=255)
   gagGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cccccc                                                                                                                               >  2:886740/1‑140 (MQ=255)
   gagGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cccccc                                                                                                                               >  2:1281301/1‑140 (MQ=255)
   gagGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cccccc                                                                                                                               >  2:1057675/1‑140 (MQ=255)
    agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cccccc‑‑‑c                                                                                                                           >  2:1339128/1‑140 (MQ=255)
    agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cc                                                                                                                                   >  1:1027281/1‑135 (MQ=255)
    agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cc                                                                                                                                   <  2:296574/135‑1 (MQ=255)
    agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cc                                                                                                                                   >  1:296574/1‑135 (MQ=255)
    agGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cc                                                                                                                                   <  2:1027281/135‑1 (MQ=255)
         tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cccccc‑‑‑cccagg                                                                                                                      <  2:255782/140‑6 (MQ=255)
         tattTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑cccccc‑‑‑cccagg                                                                                                                      >  2:229304/1‑135 (MQ=255)
           ttttAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGATTGGATTGAGGATTGCAG‑‑‑cccccc‑‑‑cccaggct                                                                                                                    >  2:1179153/1‑133 (MQ=255)
                  ttCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGc                                                                                                               <  1:1057675/140‑1 (MQ=255)
                   tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCCC‑‑‑CAGGCTTCCCAGc                                                                                                               <  2:109346/140‑1 (MQ=255)
                   tCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGca                                                                                                              <  1:1284966/140‑1 (MQ=255)
                      cGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGC‑‑CCCCCC‑‑‑CAGGCTTCCCAGCACAtt                                                                                                          <  2:914595/140‑1 (MQ=255)
                       gACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACAtt                                                                                                          <  1:886740/140‑1 (MQ=255)
                       gACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACAtt                                                                                                          <  1:948244/140‑1 (MQ=255)
                       gACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACAtt                                                                                                          <  1:1339128/140‑1 (MQ=255)
                          cAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTg                                                                                                       >  2:877473/1‑140 (MQ=255)
                           aGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGt                                                                                                      <  1:296835/140‑1 (MQ=255)
                            gggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTaa                                                                                                    <  1:998269/140‑1 (MQ=255)
                             ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTaa                                                                                                    >  1:75499/1‑140 (MQ=255)
                             ggCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTaa                                                                                                    <  1:877473/140‑1 (MQ=255)
                                  tGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGt                                                                                               >  2:391348/1‑140 (MQ=255)
                                        tATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTAccgc                                                                                         >  2:825913/1‑140 (MQ=255)
                                          tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                       <  1:391348/140‑1 (MQ=255)
                                          tCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTAccgccg                                                                                       <  2:75499/140‑1 (MQ=255)
                                           cTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACcgccgc                                                                                      >  2:1028215/1‑140 (MQ=255)
                                               cAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCtttt                                                                                  <  2:1002103/140‑1 (MQ=255)
                                                 gagTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCa                                                                                >  2:1009443/1‑140 (MQ=255)
                                                  agTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAg                                                                               <  1:1009443/140‑1 (MQ=255)
                                                    tCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCa                                                                             >  2:938586/1‑140 (MQ=255)
                                                            cgGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTc                                                                     >  2:525586/1‑140 (MQ=255)
                                                             ggCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcg                                                                    <  1:754931/140‑1 (MQ=255)
                                                              gCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑CCCCCCCCACAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTcg                                                                    <  1:369707/140‑1 (MQ=255)
                                                               cTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACAGCCGCTTTTCAGCACCGCTTTcgcg                                                                  <  1:1281301/140‑1 (MQ=255)
                                                                 aCAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGcc                                                                <  1:1179153/140‑1 (MQ=255)
                                                                   aaaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCa                                                              <  1:525586/140‑1 (MQ=255)
                                                                    aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                             <  1:737904/140‑1 (MQ=255)
                                                                    aaCATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCAt                                                             <  1:408402/140‑1 (MQ=255)
                                                                           ttGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTa                                                      >  1:1037186/1‑140 (MQ=255)
                                                                            tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTaa                                                     >  2:65661/1‑140 (MQ=255)
                                                                            tGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTaa                                                     >  1:807979/1‑140 (MQ=255)
                                                                              ccAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATg                                                   <  1:825913/140‑1 (MQ=255)
                                                                                 aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                <  1:65661/140‑1 (MQ=255)
                                                                                 aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                >  1:446919/1‑140 (MQ=255)
                                                                                 aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                <  1:436852/140‑1 (MQ=255)
                                                                                 aaCTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTg                                                >  1:1067002/1‑140 (MQ=255)
                                                                                      ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTg                                           >  1:12405/1‑140 (MQ=255)
                                                                                      ctgGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTg                                           >  1:1079577/1‑140 (MQ=255)
                                                                                              cTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTc                                   >  1:876042/1‑140 (MQ=255)
                                                                                                         tttttCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGcc                        <  2:1079577/140‑1 (MQ=255)
                                                                                                          ttttCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCt                       <  2:1228949/140‑1 (MQ=255)
                                                                                                            ttCCACTCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGa                     <  2:807979/140‑1 (MQ=255)
                                                                                                                  tCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt               <  2:11715/140‑1 (MQ=255)
                                                                                                                  tCGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGAACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACt               >  2:1211217/1‑140 (MQ=255)
                                                                                                                   cGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGcc                        >  2:244186/1‑130 (MQ=255)
                                                                                                                   cGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGcc                        <  1:244186/130‑1 (MQ=255)
                                                                                                                   cGGTTTGGATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTa              >  2:2816/1‑140 (MQ=255)
                                                                                                                         ggATTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGCAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCTCCCATTTCATTAATGCTCTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGcc        >  2:946416/1‑139 (MQ=255)
                                                                                                                           aTTGAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTgg      <  1:911478/140‑1 (MQ=255)
                                                                                                                              gAGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCt   >  2:351396/1‑140 (MQ=255)
                                                                                                                               aGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCtt  >  1:1210557/1‑140 (MQ=255)
                                                                                                                               aGGATTGCAGCCCCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCtt  <  2:1067002/140‑1 (MQ=255)
                                                                                                                                         ||                                                                                                                                    
ATAGAGGATTATTTTAAATTCCCGACCAGGGCTTTGCGGCTATCTTCCAGAGTCACAACGCGGCTACAAACATCTTTGCCAAACTGCTGGAAATCTTTTTCCTGCTTTTTCCACTCGGTTTGGATTGAGGATTGCAG‑‑‑CCCCCC‑‑‑CAGGCTTCCCAGCACATTCTGTAATGGGTTACCGCCGCTTTTCAGCACCGCTTTCGCGCCCATTTCATTAATGCTGTCCTGTAAAATTCCACCCATTGCCTGATTCACTAATTGCTGGCCTT  >  NC_000913/3100896‑3101160

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