Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,536,373 C→G 6.7% 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 6.7% 167.7 / 4.5 71intergenic (+30/+13)crr/pdxKglucose‑specific enzyme IIA component of PTS/pyridoxal‑pyridoxamine kinase/hydroxymethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base G (2/3);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

AGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACTGAT  >  NC_000913/2536226‑2536520
                                                                                                                                                   |                                                                                                                                                   
aggcAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCGCGGGTATCCGCGTCAAGAAGTAATTCTTGCCGCCGTGTAAAATGGCGCGCATCg                                                                                                                                                    >  2:875431/4‑149 (MQ=255)
 gCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCgg                                                                                                                                                   <  1:288639/149‑1 (MQ=255)
          aCCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAAGCGCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCATAAAAAAATGGCGCCGATGGGCGCCAtttt                                                                                                                                          <  1:476862/149‑1 (MQ=255)
          aCCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGAACAAACTGTCCGGTAGCGTAACCGTGGGTGAAACACCGGATATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCAtgtt                                                                                                                                          >  2:183239/1‑149 (MQ=255)
           cccTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCAttttt                                                                                                                                         <  2:428818/149‑1 (MQ=255)
           cccTGACTCCGGATCTTATCTCCTACTTTGACGACAGCAATGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGTAGTAATTCTTGCCGCAGTGAAAACTGGCGCCCATCGGCGCCAttttt                                                                                                                                         <  1:204482/149‑1 (MQ=255)
             cTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTTGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCAGCGGCGCCAttttttt                                                                                                                                       <  1:594324/149‑1 (MQ=255)
             cTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCGATGGGCGCCAtttttca                                                                                                                                       <  1:379725/149‑3 (MQ=255)
                    cGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAATATGTAATTCTTGCCGCAGTGATAAATGGCGCGACTCGGCGCCATTTTTTTATACTTc                                                                                                                                >  1:48604/1‑149 (MQ=255)
                         gttATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCa                                                                                                                           >  2:718748/1‑149 (MQ=255)
                          ttATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAg                                                                                                                          >  2:24220/1‑149 (MQ=255)
                             tctcCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGTGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAgcgg                                                                                                                       >  1:338272/1‑149 (MQ=255)
                             tctcCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAgcgg                                                                                                                       >  2:877865/1‑149 (MQ=255)
                                ccAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGcggcgg                                                                                                                    >  2:10960/1‑149 (MQ=255)
                                 cAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCCTCGGCGCCAGTTTTTTATGCTTCCGCCAGCggcggc                                                                                                                   >  1:350047/1‑149 (MQ=255)
                                  aaCATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTGTGCTTCCGCCAGCGGCGGCa                                                                                                                  >  2:760105/1‑149 (MQ=255)
                                  aaCATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCa                                                                                                                  >  2:661033/1‑149 (MQ=255)
                                   aCATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCaa                                                                                                                 >  2:793932/1‑149 (MQ=255)
                                    cATGGACGAAATCAAAGAACTGATCAAACTGTCCTGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCACATAGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCaaa                                                                                                                <  1:557426/149‑1 (MQ=255)
                                         aCGAAATCAAAGAACAGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTATATGCTTCCGCCAGCGGCGGCAAAATCaa                                                                                                           <  2:911811/149‑1 (MQ=255)
                                          cGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAAt                                                                                                          <  2:986443/149‑1 (MQ=255)
                                          cGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTATTATGCTTCCGCCAGCGGCGGCAAAATCaaa                                                                                                          >  1:596403/1‑148 (MQ=255)
                                           gAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCTGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCGATTATTTTATGATTCCGCCAGCTGCGGCAAAATCAAtt                                                                                                         >  2:568886/1‑149 (MQ=255)
                                               tCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATc                                                                                                     <  1:29324/149‑1 (MQ=255)
                                                     aaCTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCtctc                                                                                               <  1:23570/149‑1 (MQ=255)
                                                       cTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCAt                                                                                             <  1:392579/149‑1 (MQ=255)
                                                            cAAACTGTCCGGTATCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGAGCCCATAGGCGCCATTTTTTTATGCATCCGCCAGCCGCGGCAAAATAAATTCAACTCAATCAtgctgc                                                                                        >  2:102816/1‑149 (MQ=255)
                                                             aaaCTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATgctgct                                                                                       <  2:975355/149‑1 (MQ=255)
                                                              aaCTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGctgctg                                                                                      >  1:647651/1‑149 (MQ=255)
                                                               ccTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTgg                                                                                     <  1:147565/148‑1 (MQ=255)
                                                                cTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTggg                                                                                    <  1:69200/149‑1 (MQ=255)
                                                                 tGTCCGGTAGCGTACCCGGGGGGGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGgt                                                                                   <  2:602481/149‑1 (MQ=255)
                                                                   tCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGTCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGgtgt                                                                                 <  2:274925/149‑1 (MQ=255)
                                                                      gggaGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAgc                                                                              <  2:814227/146‑1 (MQ=255)
                                                                      ggTAGCGTTAACGTTGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATTGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAgc                                                                              <  1:971248/149‑1 (MQ=255)
                                                                       gTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTTATTCTTGCCGCAGTGAAAAATGGCGCCCTTCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTAATCGCTCTCATGCTGCTGGGTGTAgcg                                                                             >  1:480859/1‑149 (MQ=255)
                                                                       gTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTGTATGCTTCCGCCAGCGGCGGCAAAAACCATTGAAATCTGTAAGGGTGCTGGGTGTaccc                                                                             >  2:863776/1‑146 (MQ=255)
                                                                        tAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAgcgc                                                                            <  2:242439/149‑1 (MQ=255)
                                                                           cGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTTCCGCAGTGAAAAATGTCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGAGCATc                                                                         >  2:105942/1‑149 (MQ=255)
                                                                           cGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATc                                                                         <  1:867740/149‑1 (MQ=255)
                                                                            gTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATcc                                                                        >  1:1008810/1‑148 (MQ=255)
                                                                               aCCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACtt                                                                     >  1:654636/1‑149 (MQ=255)
                                                                                  gTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCa                                                                  >  1:534132/1‑149 (MQ=255)
                                                                                       tGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTAcg                                                             >  1:248983/1‑149 (MQ=255)
                                                                                        gAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCAGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTAcgc                                                            <  1:978129/149‑1 (MQ=255)
                                                                                         aacccccGGTTATCCGCATCAAGAAGTTATTCTTGCCGCAGTGAAAAATGGCGCCGATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCGACATCAAGTCATCGCTCTCATGCTGCTTGGTGTAGAGCATCACTGCGAGTAcgcg                                                           >  2:676048/4‑149 (MQ=255)
                                                                                           aCCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCa                                                         >  1:147629/1‑149 (MQ=255)
                                                                                            ccccGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCaa                                                        >  1:378347/1‑149 (MQ=255)
                                                                                              ccGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCACTACGCGCAAcc                                                      >  1:592898/1‑149 (MQ=255)
                                                                                                     tCCGCTTCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCGCATCGGCGCCATTTGTTTATGCTTCCGCCAGCGGCGGCAAACTCTATTCAGCGCTCTCATGCTGCTGGGTGTAGCGCAGCACTTCCCGTACGCGCACCCCCGCtag                                               >  2:191426/1‑147 (MQ=255)
                                                                                                     tCCGCATCATGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATAGGAGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCCTTCTGCTGGGTGTAGCTCAGCACTTCCAGTTAGCGCCCCCCCgaggg                                               >  2:187926/1‑145 (MQ=255)
                                                                                                     tCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCg                                               <  1:180561/149‑1 (MQ=255)
                                                                                                      ccGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCgg                                              >  1:878214/1‑149 (MQ=255)
                                                                                                        gCATCAAGAAGTAATTCTTTCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTg                                            <  2:272596/149‑1 (MQ=255)
                                                                                                        gCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTg                                            >  2:383368/1‑149 (MQ=255)
                                                                                                            cAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCGCGGTGCACt                                        <  2:68853/149‑1 (MQ=255)
                                                                                                               gaagTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCa                                     <  1:24220/149‑1 (MQ=255)
                                                                                                                  gTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCg                                  <  1:797207/149‑1 (MQ=255)
                                                                                                                    aaTTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGt                                >  2:307640/1‑149 (MQ=255)
                                                                                                                        cTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAAc                            <  2:259062/149‑1 (MQ=255)
                                                                                                                              gCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTc                      <  1:106696/149‑1 (MQ=255)
                                                                                                                                 taaaaaaaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCTTGCTGCTGGGTGTAGCGCATCCCTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCt                   >  2:202556/4‑149 (MQ=255)
                                                                                                                                  aaaaaaaTGGCGCCGATGGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCtt                  >  1:780114/3‑149 (MQ=255)
                                                                                                                                    aaataTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCGGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTATCTCCTTCCCtttg                <  1:102816/149‑2 (MQ=255)
                                                                                                                                    aaaaaTGGCGCCCATCGGCGCCAGTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCTCTTCCAGTACGCGCAACCCCGCGCGGTGCACTGCATCGGTTAACGCCTTCACTTTc                >  1:590991/1‑149 (MQ=255)
                                                                                                                                     aaaaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCAGTGCATCGGTTAACGCCTTCCCTTTCa               <  1:645914/149‑1 (MQ=255)
                                                                                                                                       aaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGTTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATTGGTTAACGCCTTCCCTTTCAGc             >  1:349484/1‑149 (MQ=255)
                                                                                                                                       aaTGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAAACATTTCATCGCTCTCATGCTGCTGGGTGTAGCGCATAACTGCCACTACGAGCAGCCCCGCTCGGTGCACTGAATCGGTAAACGCGTTCCCAGTCAgg             >  1:85269/1‑148 (MQ=255)
                                                                                                                                         tGGCGCCCATCGGCGCCATTTTTTTGTGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCaa           <  1:760105/149‑1 (MQ=255)
                                                                                                                                         tGGCGCCAATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCaa           <  1:1010440/149‑1 (MQ=255)
                                                                                                                                                  tCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACtgat  >  2:888185/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGGACGAAATCAAAGAACTGATCAAACTGTCCGGTAGCGTAACCGTGGGTGAAACCCCGGTTATCCGCATCAAGAAGTAATTCTTGCCGCAGTGAAAAATGGCGCCCATCGGCGCCATTTTTTTATGCTTCCGCCAGCGGCGGCAAAATCAATTCATCGCTCTCATGCTGCTGGGTGTAGCGCATCACTTCCAGTACGCGCAACCCCGCTCGGTGCACTGCATCGGTTAACGCCTTCCCTTTCAGCAAGCCACTGAT  >  NC_000913/2536226‑2536520

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