Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,536,373 C→G 8.0% intergenic (+30/+13) crr → / ← pdxK glucose‑specific enzyme IIA component of PTS/pyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,536,3730CG 8.0% 258.8 / 10.5 75intergenic (+30/+13)crr/pdxKglucose‑specific enzyme IIA component of PTS/pyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base G (2/4);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACTGA  >  NC_000913/2536227‑2536519
                                                                                                                                                  |                                                                                                                                                  
gCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATggg                                                                                                                                                  >  3:554170/1‑149 (MQ=255)
 ccAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGgc                                                                                                                                                 >  1:599204/1‑149 (MQ=255)
 ccAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGgc                                                                                                                                                 >  1:570002/1‑149 (MQ=255)
     gTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAACAATGGCGCCGATGGgcgcaa                                                                                                                                             <  2:494628/149‑3 (MQ=255)
     gTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCa                                                                                                                                             <  4:554170/149‑1 (MQ=255)
      tCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCAt                                                                                                                                            <  1:284232/149‑1 (MQ=255)
      tCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCAt                                                                                                                                            <  3:551154/149‑1 (MQ=255)
      tCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAt                                                                                                                                            >  3:203949/1‑149 (MQ=255)
           ccTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAtttttt                                                                                                                                       >  4:224880/1‑149 (MQ=255)
           ccTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAtttttt                                                                                                                                       >  4:132183/1‑149 (MQ=255)
           ccTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAtttttt                                                                                                                                       >  4:74449/1‑149 (MQ=255)
                   cGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTc                                                                                                                               <  3:611567/149‑1 (MQ=255)
                    ggTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTcc                                                                                                                              <  1:709997/149‑1 (MQ=255)
                     gttgttATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCg                                                                                                                             >  2:588942/1‑149 (MQ=255)
                      ttgttATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGc                                                                                                                            <  1:576213/149‑1 (MQ=255)
                        gttATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCa                                                                                                                          <  4:532676/149‑1 (MQ=255)
                              tcCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAgcggcg                                                                                                                    >  1:304857/1‑149 (MQ=255)
                              tcCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAgcggcg                                                                                                                    >  1:35352/1‑149 (MQ=255)
                               ccAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGcggcgg                                                                                                                   >  3:735871/1‑149 (MQ=255)
                                cNACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCggcggc                                                                                                                  >  3:195111/1‑149 (MQ=255)
                                 aaCATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCa                                                                                                                 >  1:652871/1‑149 (MQ=255)
                                       gACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCa                                                                                                           >  1:164172/1‑149 (MQ=255)
                                        aCGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCaa                                                                                                          <  2:628522/149‑1 (MQ=255)
                                        aCGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCaa                                                                                                          <  2:599204/149‑1 (MQ=255)
                                          gAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGc                                                                                                                                  >  3:377922/1‑123 (MQ=255)
                                          gAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGc                                                                                                                                  <  4:377922/123‑1 (MQ=255)
                                            aaTCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCa                                                                                                      >  4:282501/1‑149 (MQ=255)
                                              tCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATc                                                                                                    >  4:326916/1‑149 (MQ=255)
                                                aaaGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCaa                                                                                                          >  2:44887/1‑141 (MQ=255)
                                                aaaGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCaa                                                                                                          <  1:44887/141‑1 (MQ=255)
                                                  aGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGctc                                                                                                <  1:723899/149‑1 (MQ=255)
                                                  aGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGctc                                                                                                >  1:73717/1‑149 (MQ=255)
                                                      cTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAt                                                                                            <  1:248244/149‑1 (MQ=255)
                                                        gATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAtgc                                                                                          >  4:502888/1‑149 (MQ=255)
                                                         aTCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAtgct                                                                                         >  2:237486/1‑149 (MQ=255)
                                                         aTCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAtgct                                                                                         >  2:72069/1‑149 (MQ=255)
                                                              aCTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTgg                                                                                    <  2:304857/149‑1 (MQ=255)
                                                              aCTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTgg                                                                                    <  3:224880/149‑1 (MQ=255)
                                                                tGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGgt                                                                                  >  3:407020/1‑149 (MQ=255)
                                                                tGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGgt                                                                                  <  2:164172/149‑1 (MQ=255)
                                                                  tCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGgtgt                                                                                >  2:128531/1‑149 (MQ=255)
                                                                               ccGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTc                                                                   <  3:675097/149‑1 (MQ=255)
                                                                                cGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTcc                                                                  <  3:253871/149‑1 (MQ=255)
                                                                                 gTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCa                                                                 <  2:73717/149‑1 (MQ=255)
                                                                                       gAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTAcgc                                                           <  3:282501/149‑1 (MQ=255)
                                                                                        aaaCCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGAGTAGCGCATCACTTCCAGTAcgcg                                                          >  3:717950/1‑149 (MQ=255)
                                                                                            cccGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAAc                                                      >  3:257717/1‑149 (MQ=255)
                                                                                             ccGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAAcc                                                     <  3:37627/149‑1 (MQ=255)
                                                                                                  tATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCt                                                <  4:189829/149‑1 (MQ=255)
                                                                                                     ccGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCgg                                             >  2:190780/1‑149 (MQ=255)
                                                                                                       gCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGTCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCCGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTAGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTg                                           >  4:264930/1‑149 (MQ=255)
                                                                                                        cATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAACGCATCACATCCAGTACGCGCAACCCCGCTCGGTGc                                          >  2:526422/1‑149 (MQ=255)
                                                                                                          tCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCAc                                        <  4:2406/149‑1 (MQ=255)
                                                                                                          tCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCAc                                        <  3:123092/149‑1 (MQ=255)
                                                                                                            aagaagTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTg                                      >  3:144754/1‑149 (MQ=255)
                                                                                                              gaagTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCa                                    >  1:745184/1‑149 (MQ=255)
                                                                                                                agTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATc                                  <  4:195111/149‑1 (MQ=255)
                                                                                                                 gTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCg                                 >  1:712270/1‑149 (MQ=255)
                                                                                                                   aaTTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGt                               <  3:502888/149‑1 (MQ=255)
                                                                                                                      tCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTaa                            <  1:570537/149‑1 (MQ=255)
                                                                                                                       cTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAAc                           <  3:132183/149‑1 (MQ=255)
                                                                                                                             gCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTc                     >  1:535775/1‑149 (MQ=255)
                                                                                                                               aGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGATGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTccc                   <  4:735871/149‑1 (MQ=255)
                                                                                                                                   aaaaaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTc               <  4:203949/149‑1 (MQ=255)
                                                                                                                                   aaaaaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCACCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTc               <  1:65748/149‑1 (MQ=255)
                                                                                                                                     aaaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTTCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAg             >  2:678196/1‑149 (MQ=255)
                                                                                                                                      aaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGc            >  1:387958/1‑149 (MQ=255)
                                                                                                                                       aTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCa           >  3:50312/1‑149 (MQ=255)
                                                                                                                                       aTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCa           >  4:575508/1‑149 (MQ=255)
                                                                                                                                            gcCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCa      >  1:1702/1‑149 (MQ=255)
                                                                                                                                            gcCCATCGGCGCCATTTTTTTATGCTTCCGCCAGAGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCa      >  2:258209/1‑149 (MQ=255)
                                                                                                                                             cccATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCAc     <  2:387958/149‑1 (MQ=255)
                                                                                                                                              ccATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCgg                                <  4:357091/121‑1 (MQ=255)
                                                                                                                                              ccATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCgg                                >  3:357091/1‑121 (MQ=255)
                                                                                                                                              ccATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACt    <  1:326309/149‑1 (MQ=255)
                                                                                                                                                aTCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACtga  <  1:118794/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACTGA  >  NC_000913/2536227‑2536519

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