Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,883 +C coding (667/687 nt) rph ← truncated RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,8831.C95.1% 126.5 / ‑2.3 42coding (667/687 nt)rphtruncated RNase PH
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (24/15);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815633‑3816132
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
gATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑Tc‑ccc                                                                                                                                                                                                                                                            >  1:137349/1‑251 (MQ=255)
         cGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACGTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCa                                                                                                                                                                                                                                                   <  1:297539/251‑1 (MQ=255)
                              tGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAACCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCCGCGCCGCCTTCTGCGTAGCTAAAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGTTGAGTCGct                                                                                                                                                                                                                              >  2:79386/1‑251 (MQ=255)
                                  gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTCCAGACAAAAAAAAGCCGCCTCATCAGTCGCCTTAAACATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCAAACAAGATGAGTAGCTATTc                                                                                                                                                                                                                          >  2:53393/1‑251 (MQ=255)
                                                           gtCAGCGAAAACACGCCAAACTTTAACACCTTCTTGCAAAGAGAAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTAATTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTGTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCACGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCGCGCCTTCTGcc                                                                                                                                                                                                 <  2:271432/250‑1 (MQ=255)
                                                              cgCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCACATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACggatc                                                                                                                                                                                                           >  2:241217/2‑235 (MQ=255)
                                                              aGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                                                                                                                                                                                                           <  1:241217/238‑1 (MQ=255)
                                                                cGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGctct                                                                                                                                                                                                                            >  1:334370/1‑219 (MQ=255)
                                                                cGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGctct                                                                                                                                                                                                                            <  2:334370/219‑1 (MQ=255)
                                                                 gtaatcgcGCCAAACTTTAACACCTTCTTGCTAAGAGTAAACTCAATAAACTGGCGCTGATTTGGTTTCATGCCGTCGCTCCTCATCTTACGTTTCTACAGACAAAAATAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATGCC‑CCC‑TCGGGCCAGAGACAACAAGATGAGTAGCTCGTCATGGGTGTACGGCTCGCCTTCTGCCGTcccc                                                                                                                                                                                           <  2:344340/244‑1 (MQ=255)
                                                                   aaaCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTGGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGCGACTCCTCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGAAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg                                                                                                                                                                                         >  1:49580/1‑251 (MQ=255)
                                                                      cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGAGAACGGCTCGCCTTCTGCCGTCCCCTGCAc                                                                                                                                                                                      >  1:298497/1‑251 (MQ=255)
                                                                            aaaCTTTAACCCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTGCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAat                                                                                                                                                                                <  2:137349/251‑1 (MQ=255)
                                                                                     cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCGACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAGCAGTTTGCCAGCTCCGCCTGCTGCGTCGCTACAATGGATGCGAT‑TC‑CCCCTCGGGCCAGATCCAACAAGATGAGTAGCTCTGCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGccc                                                                                                                                                                       <  1:326626/251‑1 (MQ=255)
                                                                                                    cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCAc                                                                                                                                                        >  2:279461/1‑251 (MQ=255)
                                                                                                                 aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACCAGATGAGTAGCTCTTCATg                                                                                                                                                                                                                       >  2:40100/1‑175 (MQ=255)
                                                                                                                 aaCTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATg                                                                                                                                                                                                                       <  1:40100/175‑1 (MQ=255)
                                                                                                                     ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGt                                                                                                                                         <  1:53393/249‑1 (MQ=255)
                                                                                                                         cTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGc                                                                                                                                   >  1:274936/1‑251 (MQ=255)
                                                                                                                                  tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTc                                                                                                                          >  1:9182/1‑251 (MQ=255)
                                                                                                                                        gCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGCCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGCGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTTCCGTCCCCTGCACTTCCATGATGCGCCCGTCTTCGGTCATCACTACGTTCCTGTCGGTCTCCGTGGCCGAGTCTTccac                                                                                                                    >  2:248171/1‑248 (MQ=255)
                                                                                                                                             cgatcctcATCTTACTTTTCGACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAACTCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGCTTCGAT‑TC‑CCCATCGGGCCAGAGCCAACAAGCTGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGAGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGCGTCTTCAAcctctt                                                                                                               >  2:94396/4‑246 (MQ=255)
                                                                                                                                                        ttACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCCGATCGCa                                                                                                     >  1:194597/1‑250 (MQ=255)
                                                                                                                                                          aCTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCaaa                                                                                                   <  1:279003/250‑1 (MQ=255)
                                                                                                                                                            ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGt                                                                                                                                                                                               <  2:80556/156‑1 (MQ=255)
                                                                                                                                                            ttttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGt                                                                                                                                                                                               >  1:80556/1‑156 (MQ=255)
                                                                                                                                                             tttCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGc                                                                                               >  1:32098/1‑251 (MQ=255)
                                                                                                                                                                         aaaaaaaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTcac                                                                                   >  1:129551/1‑251 (MQ=255)
                                                                                                                                                                              aaGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTc                                                                              >  1:210460/1‑251 (MQ=255)
                                                                                                                                                                               aGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTcc                                                                             <  1:67753/251‑1 (MQ=255)
                                                                                                                                                                                 gCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATAGCAAACCGCTTCGCCGTTCACAATTcctc                                                                           >  2:296951/1‑249 (MQ=255)
                                                                                                                                                                                       catcaGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGaa                                                                      >  1:82488/1‑250 (MQ=255)
                                                                                                                                                                                            gTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCTCTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCg                                                                <  1:323728/251‑1 (MQ=255)
                                                                                                                                                                                              cGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCgg                                                               <  1:279461/250‑1 (MQ=255)
                                                                                                                                                                                                            gTTTGACAGAGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGATCATCACTACGTTCATATCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGCCAGCAACTGCGTCTACCATCCCTTTcc                                                >  2:15291/1‑250 (MQ=255)
                                                                                                                                                                                                               tGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATCCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGAAAACCGCTTCGCCGTTCACAATTCCGACAGACACTGCGGCCAACATCCCTTTCATCg                                             >  2:43144/1‑251 (MQ=255)
                                                                                                                                                                                                                gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCgg                                            >  1:90671/1‑251 (MQ=255)
                                                                                                                                                                                                                gCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCg                                             <  1:255272/250‑1 (MQ=255)
                                                                                                                                                                                                                       ccgccTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGtt                                     >  1:203847/1‑251 (MQ=255)
                                                                                                                                                                                                                         gccTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCCCTCGCGCCAGCGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCCATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGCCTCTCCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCCTCGCCGTTCCCAATTCCGCCAGCAACTCCGGCTTCCATCCCTTTCATCGCGTTGGtttt                                   >  1:28829/1‑251 (MQ=255)
                                                                                                                                                                                                                                gCGTAGCTACAATGGACTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCtg                             <  2:82488/250‑2 (MQ=255)
                                                                                                                                                                                                                                  gTCGCTACAATGGATTCGAT‑TC‑CCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCACACCGCTTCGCCGTTCACAATTCAGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGcc                          >  2:131331/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aCAATGGTGTCGAT‑TCTCCC‑TCGGGCCAGATCCAACAAGATGATTAGCGCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCTCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCGTTGCGGCAGTGTCTTCAACGTATTCCAGATGGGAAACCGCTTCGCCGTTCACAATTCCTACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTc                    <  1:234663/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       cc‑ccc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAGCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAg     <  1:66950/250‑1 (MQ=255)
                                                                                                                                                                                                                                                       cc‑ccc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGt       >  1:257154/2‑249 (MQ=255)
                                                                                                                                                                                                                                                           cc‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAgcgc  >  1:102342/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         
GATCGCGCCCGGTATTAAACAGCCCGGCGTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGAT‑TC‑CCC‑TCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGC  >  NC_000913/3815633‑3816132

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑