Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,346,792 C→G 5.6% P72R (CCG→CGG)  ptsN → sugar‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,346,7920CG 5.6% 247.4 / 5.9 71P72R (CCG→CGG) ptsNsugar‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (34/33);  new base G (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

TCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCGGCAGACCAAACTAAAACGCAC  >  NC_000913/3346645‑3346937
                                                                                                                                                   |                                                                                                                                                 
tCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCg                                                                                                                                                  >  3:160470/1‑149 (MQ=255)
  cACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCa                                                                                                                                                >  2:247745/1‑149 (MQ=255)
    cTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATg                                                                                                                                              <  2:485437/149‑1 (MQ=255)
      gTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGc                                                                                                                                            <  3:148563/149‑1 (MQ=255)
         agagCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCaaa                                                                                                                                         <  3:358688/149‑1 (MQ=255)
          gagCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATGCGGAATGGCaata                                                                                                                                        <  3:456175/149‑3 (MQ=255)
               aaaaaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTACCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGc                                                                                                                                                 <  1:284434/135‑1 (MQ=255)
               aaaaaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTACCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGc                                                                                                                                                 >  2:284434/1‑135 (MQ=255)
                 aaaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGgaaga                                                                                                                                 >  3:246795/1‑149 (MQ=255)
                  aaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCa                                                                                                                                                >  1:739227/1‑133 (MQ=255)
                  aaCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCa                                                                                                                                                <  2:739227/133‑1 (MQ=255)
                   aCGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGaagaag                                                                                                                               <  4:121876/149‑1 (MQ=255)
                      cgcCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATa                                                                                                                            <  4:674005/149‑1 (MQ=255)
                                   atcaGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATGCGGAATGGCaat                                                                                                                                         >  1:300797/1‑122 (MQ=255)
                                   atcaGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATGCGGAATGGCaat                                                                                                                                         <  2:300797/123‑2 (MQ=255)
                                    tcaGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGtt                                                                                                              <  3:335644/149‑1 (MQ=255)
                                    tcaGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGtt                                                                                                              >  1:513739/1‑149 (MQ=255)
                                    tcaGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGtt                                                                                                              >  3:182170/1‑149 (MQ=255)
                                     caGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTg                                                                                                             <  4:366108/149‑1 (MQ=255)
                                            cTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGtttt                                                                                                      <  4:627214/149‑1 (MQ=255)
                                                 ggcgAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGgaa                                                                                                                                   >  1:635644/1‑115 (MQ=255)
                                                 ggcgAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGgaa                                                                                                                                   <  2:635644/115‑1 (MQ=255)
                                                 ggcgAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTc                                                                                                 <  2:537172/149‑1 (MQ=255)
                                                 ggcgAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTc                                                                                                 >  4:589399/1‑149 (MQ=255)
                                                  gcgAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCa                                                                                                >  1:87832/1‑149 (MQ=255)
                                                    gAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGc                                                                                              >  1:189502/1‑149 (MQ=255)
                                                     aaaCAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCt                                                                                             >  1:230846/1‑149 (MQ=255)
                                                       acaacTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCg                                                                                           >  1:380536/1‑149 (MQ=255)
                                                          acTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGaaa                                                                                        <  1:551708/149‑1 (MQ=255)
                                                            ttATCCTGCCGCCTCAGGTGGTTGTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGCCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACg                                                                                      >  4:724481/1‑149 (MQ=255)
                                                            ttAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACg                                                                                      >  2:509930/1‑149 (MQ=255)
                                                             tAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGc                                                                                     <  1:478314/149‑1 (MQ=255)
                                                                 cTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGtt                                                                                                  >  1:611415/1‑132 (MQ=255)
                                                                 cTGCCGCCTAAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGtt                                                                                                  <  2:611415/132‑1 (MQ=255)
                                                                    ccgccTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGc                                                                              >  1:630366/1‑149 (MQ=255)
                                                                    ccgccTCAGGCGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGc                                                                              >  3:48563/1‑149 (MQ=255)
                                                                          cAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGa                                                                        <  4:182170/149‑1 (MQ=255)
                                                                             gtggtTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTg                                                                                                             >  1:193188/1‑109 (MQ=255)
                                                                             gtggtTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTg                                                                                                             <  2:193188/109‑1 (MQ=255)
                                                                                   tttGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCGGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGa                                                               >  3:472706/1‑149 (MQ=255)
                                                                                      gAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGAcaa                                                            <  4:535874/149‑1 (MQ=255)
                                                                                       aaGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGAcaac                                                           <  3:86699/149‑1 (MQ=255)
                                                                                             tCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCg                                                     >  1:92118/1‑149 (MQ=255)
                                                                                              ccTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATGCTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCgg                                                    >  3:354450/1‑149 (MQ=255)
                                                                                                   cgcgGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGac                                               <  2:694744/149‑1 (MQ=255)
                                                                                                   cgcgGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGac                                               <  2:513739/149‑1 (MQ=255)
                                                                                                    gcgGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGact                                              <  1:361912/149‑1 (MQ=255)
                                                                                                      gggAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGactta                                            >  2:296272/1‑149 (MQ=255)
                                                                                                        gAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGacttact                                          >  4:706650/1‑149 (MQ=255)
                                                                                                         aaaaaaTGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGacttactt                                         <  1:414007/149‑1 (MQ=255)
                                                                                                            aaaTGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACttttt                                      <  3:248480/149‑1 (MQ=255)
                                                                                                            aaaTGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACttttt                                      <  1:668080/149‑1 (MQ=255)
                                                                                                             aaTGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTg                                     >  1:364712/1‑149 (MQ=255)
                                                                                                             aaTGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTg                                     >  2:496979/1‑149 (MQ=255)
                                                                                                                gggCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGccc                                  <  1:247745/149‑1 (MQ=255)
                                                                                                                  gCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCctg                                >  4:569841/1‑149 (MQ=255)
                                                                                                                  gCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGATTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCctg                                >  2:692921/1‑149 (MQ=255)
                                                                                                                         cGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGc                         >  2:369134/1‑149 (MQ=255)
                                                                                                                          ggTATCGGCAATGGTATTGCTATTCCGCAAGGCAAACTGGAAGAAGCTACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAACCGCCTATAGCTTTCGATGCAATCGCCAGCCAACCGGTAGACTTACTATTAGCTGTCATGGTGcc                        <  2:519369/149‑1 (MQ=255)
                                                                                                                           gTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCg                       >  1:178899/1‑149 (MQ=255)
                                                                                                                              tCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCttt                                                                           >  3:345494/1‑94 (MQ=255)
                                                                                                                              tCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCttt                                                                           <  4:345494/94‑1 (MQ=255)
                                                                                                                               cGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCGGCAg                   >  3:2614/1‑149 (MQ=255)
                                                                                                                                 gCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCGGCAGAc                 >  1:467852/1‑149 (MQ=255)
                                                                                                                                     tGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCGGCAGACCaaa             <  4:149903/149‑1 (MQ=255)
                                                                                                                                     tGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCGGCAGACCaaa             >  1:285307/1‑149 (MQ=255)
                                                                                                                                     tGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCGGCAGACCaaa             >  3:669035/1‑149 (MQ=255)
                                                                                                                                           tGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCGGCAGACCAAACTaaaa       <  1:750727/149‑1 (MQ=255)
                                                                                                                                            gCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCGGCAGACCAAACTAAAAc      <  3:262621/149‑1 (MQ=255)
                                                                                                                                                ttCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCGGCAGACCAAACTAAAACGCAc  <  4:2614/149‑1 (MQ=255)
                                                                                                                                                ttCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCGGCAGACCAAACTAAAACGCAc  <  3:706650/149‑1 (MQ=255)
                                                                                                                                                   cGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCAt                                                                  <  4:327772/82‑1 (MQ=255)
                                                                                                                                                   cGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCAt                                                                  >  3:327772/1‑82 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCCACTGTCAGAGCAAAAAACGCGCCCTGGAAATCATCAGCGAGCTGGCGGCGAAACAACTTAGCCTGCCGCCTCAGGTGGTTTTTGAAGCTATCCTGACGCGGGAAAAAATGGGCAGTACCGGTATCGGCAATGGTATTGCCATTCCGCATGGCAAACTGGAAGAAGATACTCTGCGCGCCGTTGGCGTTTTCGTTCAGCTCGAAACGCCTATAGCTTTCGATGCCATCGACAACCAACCGGTGGACTTACTTTTTGCCCTGCTGGTGCCGGCAGACCAAACTAAAACGCAC  >  NC_000913/3346645‑3346937

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