Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 371,999 T→A 5.5% I254N (ATT→AAT)  mhpC → 2‑hydroxy‑6‑ketonona‑2,4‑dienedioic acid hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913371,9990TA 5.5% 323.6 / 7.6 91I254N (ATT→AAT) mhpC2‑hydroxy‑6‑ketonona‑2,4‑dienedioic acid hydrolase
Reads supporting (aligned to +/- strand):  ref base T (40/46);  new base A (3/2);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

CTTCGTTAAGAGCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCG  >  NC_000913/371853‑372147
                                                                                                                                                  |                                                                                                                                                    
cTTCGTTAAGAGCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCAtt                                                                                                                                                     <  1:451910/148‑1 (MQ=255)
 ttCGTTAAGAGCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGc                                                                                                                                                   >  3:546969/1‑149 (MQ=255)
 ttCGTTAAGAGCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGc                                                                                                                                                   <  4:245414/149‑1 (MQ=255)
 ttCGTTAAGAGCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGATTTGc                                                                                                                                                   <  4:288994/149‑1 (MQ=255)
 ttCGTTAAGAGCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCACAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGc                                                                                                                                                   >  3:688630/1‑149 (MQ=255)
      tAAGAGCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGt                                                                                                                                               <  1:825399/148‑1 (MQ=255)
       aaGAGCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGtt                                                                                                                                              <  1:138918/148‑1 (MQ=255)
       aaGAGCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCAATGCCggaca                                                                                                                                             <  4:106787/149‑4 (MQ=255)
           gCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGaa                                                                                                                                         <  4:406500/149‑1 (MQ=255)
           gCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGa                                                                                                                                          <  1:905294/148‑1 (MQ=255)
              tGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAAc                                                                                                                                        <  3:133411/147‑1 (MQ=255)
              tGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAAc                                                                                                                                        >  4:133411/1‑147 (MQ=255)
              tGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTg                                                                                                                                      >  4:292664/1‑149 (MQ=255)
               ggAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGc                                                                                                                                     >  3:788648/1‑149 (MQ=255)
                gAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCa                                                                                                                                    >  4:394580/1‑149 (MQ=255)
                       aTCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCtt                                                                                                                              >  1:322645/1‑148 (MQ=255)
                        tCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTcc                                                                                                                            >  4:565151/1‑149 (MQ=255)
                        tCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTc                                                                                                                             >  2:30764/1‑148 (MQ=255)
                         ccGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCAATGCCggac                                                                                                                                              <  2:911419/130‑3 (MQ=255)
                         ccGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCAATGCCggac                                                                                                                                              >  1:911419/1‑128 (MQ=255)
                             aaCAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGAc                                                                                                                       <  3:489104/149‑1 (MQ=255)
                                 gTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACtgtg                                                                                                                   >  3:106044/1‑149 (MQ=255)
                                  ttCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACtgtg                                                                                                                   <  2:683578/148‑1 (MQ=255)
                                     ccGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCa                                                                                                               >  4:826617/1‑149 (MQ=255)
                                        gATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACt                                                                                                             <  2:207877/148‑1 (MQ=255)
                                          ttttGGCCCACGTCTGGCGGAAATCAAAGCGCAAACGCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTggg                                                                                                          >  3:761842/1‑149 (MQ=255)
                                          ttttGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTGTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTggg                                                                                                          >  3:638206/1‑149 (MQ=255)
                                               gCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACtgt                                                                                                                    <  4:876301/134‑1 (MQ=255)
                                               gCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACtgt                                                                                                                    >  3:876301/1‑134 (MQ=255)
                                                 ccACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGt                                                                                                    >  2:115061/1‑148 (MQ=255)
                                                 ccACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTg                                                                                                   >  3:394700/1‑149 (MQ=255)
                                                  cACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTg                                                                                                   <  2:278308/148‑1 (MQ=255)
                                                    cGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTggg                                                                                                 >  1:654041/1‑148 (MQ=255)
                                                     gTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGaa                                                                                               <  4:394700/149‑1 (MQ=255)
                                                     gTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGaa                                                                                               <  4:466311/149‑1 (MQ=255)
                                                     gTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGa                                                                                                <  2:482521/148‑1 (MQ=255)
                                                          gCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGc                                                                                          <  4:761842/149‑1 (MQ=255)
                                                             gAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCg                                                                                        <  2:419654/148‑1 (MQ=255)
                                                                    taGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGAGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCATTGGGCGCAGTGGGAACGTGCCGACGCTTtc                                                                                >  4:213889/2‑149 (MQ=255)
                                                                    aaGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCttt                                                                                 >  1:326113/1‑148 (MQ=255)
                                                                        gcAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTtcaatc                                                                            >  4:513336/1‑149 (MQ=255)
                                                                             cccTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACt                                                                        <  2:811503/148‑1 (MQ=255)
                                                                              ccTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTgg                                                                      >  4:335468/1‑149 (MQ=255)
                                                                                tGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTACATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGt                                                                     >  1:262867/1‑148 (MQ=255)
                                                                                  aTTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACtg                                                                                                                     <  3:847318/98‑1 (MQ=255)
                                                                                  aTTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACtg                                                                                                                     >  4:847318/1‑98 (MQ=255)
                                                                                   ttGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCt                                                                  >  1:114739/1‑148 (MQ=255)
                                                                                     gTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGaa                                                               <  3:394580/149‑1 (MQ=255)
                                                                                      tCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGaa                                                               <  1:612047/148‑1 (MQ=255)
                                                                                          ggggCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTcc                                                          >  4:847326/1‑149 (MQ=255)
                                                                                           gggCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTcc                                                          >  2:558514/1‑148 (MQ=255)
                                                                                             gcgcAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCg                                                       >  3:175682/1‑149 (MQ=255)
                                                                                              cgcAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACt                                                                                                             >  2:299903/1‑94 (MQ=255)
                                                                                              cgcAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACt                                                                                                             <  1:299903/94‑1 (MQ=255)
                                                                                                  aCGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGc                                                  <  4:369816/149‑1 (MQ=255)
                                                                                                    gACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGccc                                                >  3:410502/1‑149 (MQ=255)
                                                                                                     aCCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGccc                                                >  2:805545/1‑148 (MQ=255)
                                                                                                      ccGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCtt                                              >  4:823521/1‑149 (MQ=255)
                                                                                                      ccGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCt                                               <  1:805545/148‑1 (MQ=255)
                                                                                                      ccGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCt                                               <  2:548091/148‑1 (MQ=255)
                                                                                                         cTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAg                                           <  3:71052/149‑1 (MQ=255)
                                                                                                           ttGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAgg                                          >  2:664073/1‑148 (MQ=255)
                                                                                                             gtgCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAAt                                       <  3:335468/149‑1 (MQ=255)
                                                                                                              tgCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAAt                                       <  1:393457/148‑1 (MQ=255)
                                                                                                                 cGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTc                                   <  4:410502/149‑1 (MQ=255)
                                                                                                                   atggatgCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTc                                   <  3:561313/147‑1 (MQ=255)
                                                                                                                   atggatgCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTc                                   >  4:561313/1‑147 (MQ=255)
                                                                                                                         gCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTtcaa                                                                              <  1:27625/98‑1 (MQ=255)
                                                                                                                         gCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTtcaa                                                                              >  2:27625/1‑98 (MQ=255)
                                                                                                                         gCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGaa                           <  3:88345/149‑1 (MQ=255)
                                                                                                                         gCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGaa                           <  3:513336/149‑1 (MQ=255)
                                                                                                                         gCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGaa                           <  4:788648/149‑1 (MQ=255)
                                                                                                                         gCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGa                            <  2:758939/148‑1 (MQ=255)
                                                                                                                         gCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGa                            <  2:114739/148‑1 (MQ=255)
                                                                                                                          cGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAg                          <  4:106044/149‑1 (MQ=255)
                                                                                                                             gtctgcgtctgcTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCAt                       >  3:454005/1‑149 (MQ=255)
                                                                                                                                 gcgtctgcTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATActc                   <  3:823521/149‑1 (MQ=255)
                                                                                                                                  cgtctgcTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTaa                                            >  1:116683/1‑123 (MQ=255)
                                                                                                                                  cgtctgcTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTaa                                            <  2:116683/123‑1 (MQ=255)
                                                                                                                                     ctgctgTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGaa                                                               >  2:353662/1‑101 (MQ=255)
                                                                                                                                     ctgctgTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGaa                                                               <  1:353662/101‑1 (MQ=255)
                                                                                                                                     ctgctgTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTg                <  2:322645/148‑1 (MQ=255)
                                                                                                                                     ctgctgTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGa               <  4:384667/149‑1 (MQ=255)
                                                                                                                                      tgctgTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAg              <  4:454005/149‑1 (MQ=255)
                                                                                                                                       gctgTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAg              >  2:910065/1‑148 (MQ=255)
                                                                                                                                         gaaccGGCAATGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCa            >  1:734267/4‑148 (MQ=255)
                                                                                                                                          gTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCaa           <  2:262867/148‑1 (MQ=255)
                                                                                                                                          aaccGGCAATGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCaa           >  1:524511/3‑148 (MQ=255)
                                                                                                                                             cGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTgg                                                                                                                  >  3:183080/1‑42 (MQ=255)
                                                                                                                                             cGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTgg                                                                                                                  <  4:183080/42‑1 (MQ=255)
                                                                                                                                              ggCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTgg       >  1:43507/1‑148 (MQ=255)
                                                                                                                                                 aTTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTggcggc   <  4:54458/149‑1 (MQ=255)
                                                                                                                                                  ttGCCGGTTCTGAACTGCATATCTTCCGCTACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGgcggcg  >  3:737803/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTCGTTAAGAGCCTGGAAGCTAATCCGAAACAGTTCCCGGATTTTGGCCCACGTCTGGCGGAAATCAAAGCGCAAACCCTGATTGTCTGGGGGCGCAACGACCGCTTTGTGCCGATGGATGCGGGTCTGCGTCTGCTGTCCGGCATTGCCGGTTCTGAACTGCATATCTTCCGCGACTGTGGTCACTGGGCGCAGTGGGAACATGCCGACGCTTTCAATCAACTGGTGCTGAATTTCCTCGCACGCCCTTAAGGAATGGTCATGACGAAGCATACTCTTGAGCAACTGGCGGCG  >  NC_000913/371853‑372147

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: