Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,020,511 C→G 6.2% intergenic (+22/‑29) ssnA → / → ygfM putative chlorohydrolase/aminohydrolase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,020,5110CG 6.2% 209.7 / 4.6 78intergenic (+22/‑29)ssnA/ygfMputative chlorohydrolase/aminohydrolase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (35/37);  new base G (3/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCC  >  NC_000913/3020363‑3020657
                                                                                                                                                     |                                                                                                                                                  
caGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGc                                                                                                                                                    >  1:305547/1‑149 (MQ=255)
caGCGTGATGCTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGc                                                                                                                                                    >  2:1045979/1‑149 (MQ=255)
 aGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCa                                                                                                                                                   >  1:537411/1‑149 (MQ=255)
  gCGTGATGGTCAATGGTGTGATGGGCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAg                                                                                                                                                  <  1:1094659/149‑1 (MQ=255)
   cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAgg                                                                                                                                                 <  2:845931/149‑1 (MQ=255)
   cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAgg                                                                                                                                                 <  1:408855/149‑1 (MQ=255)
   cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAgg                                                                                                                                                 <  2:215422/149‑1 (MQ=255)
     tGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGa                                                                                                                                               <  1:446036/149‑1 (MQ=255)
      gATTGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTATAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGaa                                                                                                                                              <  2:1067510/149‑1 (MQ=255)
      gATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGCATGGATGCGCTGGCATAAATGACAGATGCCCTCTACCCGCAGGGaa                                                                                                                                              >  1:176322/1‑149 (MQ=255)
         ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAagag                                                                                                                                           >  1:1130825/1‑149 (MQ=255)
          gTCAATGGTGTGATGGTCTATGAAGACCGGCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAgg                                                                                                                                          >  2:1087496/1‑149 (MQ=255)
           tCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCCCCGGGATGAggt                                                                                                                                         >  2:281898/1‑148 (MQ=255)
            cAATGGTGTGATTGTCTATGAAGACCGTCAGTATAACTTCTATTGCGAATCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGGCGGATGGATGCGCTGTCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgggc                                                                                                                                        <  2:334498/149‑2 (MQ=255)
            cAATGGTGTGATGGGCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCTGCGGGAAGAgggc                                                                                                                                        <  1:871499/149‑2 (MQ=255)
              aTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTCCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGACAGTATGTGGCGTCGGATCGATGCGCTGGCATAAATGACAGATGCCCTCTACCCGGAGGGAAGAGGGTTa                                                                                                                                      >  1:640313/1‑149 (MQ=255)
                      aTGGTCTATGAAGACCGTCAGATAAACTTCGATTGCTATTCCAATAATGAACAATCCAGAAAAGAAGCTGCCAGTATCTGGCGAAGTATGGATGACATGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAATGa                                                                                                                              >  2:828261/1‑149 (MQ=255)
                           cTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaagaa                                                                                                                         >  1:216030/1‑149 (MQ=255)
                              tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGATGAGGGTTAAATAAGGAAAGAAgat                                                                                                                      >  2:655498/1‑149 (MQ=255)
                              tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                                                                                                                      >  2:183734/1‑149 (MQ=255)
                              tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                                                                                                                      >  1:429227/1‑149 (MQ=255)
                              tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                                                                                                                      <  2:229417/149‑1 (MQ=255)
                                aaGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga                                                                                                                    <  1:631368/149‑1 (MQ=255)
                                aaGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga                                                                                                                    >  1:41882/1‑149 (MQ=255)
                                 aGACCGTCAGGTTAACTACGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                                   >  1:299595/1‑149 (MQ=255)
                                   aCCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTg                                                                                                                 <  1:576409/149‑1 (MQ=255)
                                    ccGTCAGTTTAAATTCGATTGCTATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGa                                                                                                                <  2:943999/149‑1 (MQ=255)
                                     cGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCTGTATGTGGCGTAGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGATAGAAGATGATTGaa                                                                                                               >  2:403516/1‑149 (MQ=255)
                                     cGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                               >  1:442166/1‑149 (MQ=255)
                                       tCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACa                                                                                                             <  1:520536/149‑1 (MQ=255)
                                           tttAACTTCGATTGCGAGTCCTTTTATGCAGAAGCCAGATAAGCCGCTGCCAGTATGGGGCGTCGGATGGATGCGCTGGCATAATTGACAGATGCCGTCTTCCCGCAGGGTTGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttt                                                                                                         <  1:160978/149‑1 (MQ=255)
                                             taccTTCGATTGCGATTCCATTTATGCCCAAGCCAGAAACGCCGCTGCAAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCAGCAGGGAAGAGGGTTGAATAAGGAAAGAATATGATTGAACAAttttt                                                                                                       <  2:809051/146‑1 (MQ=255)
                                             tAACTTCGATTGCGATTCCATTTATGCACAAGCCTGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttttt                                                                                                       <  1:1043014/149‑1 (MQ=255)
                                              aaCTTCGTTTGCTATTCCATTTATGCACAAGCCAGACAAGCCGCTGCCAGGATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTc                                                                                                      <  2:144784/149‑1 (MQ=255)
                                                cTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAg                                                                                                    <  1:127630/149‑1 (MQ=255)
                                                  tCGATTGCGATTCCATTTATTCACAAGCCAGAATAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTTGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTTAACAATTTTTCAGGc                                                                                                  <  2:569337/149‑1 (MQ=255)
                                                      ttGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAAGAAGGAAAGAAGATGATTGACCAATTTTTCAGGCCCGa                                                                                              >  2:1141281/1‑149 (MQ=255)
                                                            ttCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGt                                                                                        <  1:124587/149‑1 (MQ=255)
                                                             tCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTc                                                                                       <  1:358635/149‑1 (MQ=255)
                                                                aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                                    <  2:533645/149‑1 (MQ=255)
                                                                   tATGAACATGACATAAAAGCCGCTTCCAGTATGTGGCGGCGGTTGGTTGTGCTGGCATTAATGACAGATGCCCTCTTCCCGCAGGGAAGAGTGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCCGGCGCGACTCTGTCGAACAg                                                                                 <  2:1111303/149‑1 (MQ=255)
                                                                       cacaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAATTAAGGAAAGAAGATGTTTGAACAATTTTTCAGGCCCGACTCTGTCGAGCCGgcgc                                                                             >  1:563892/1‑149 (MQ=255)
                                                                        acaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCt                                                                            >  2:637188/1‑149 (MQ=255)
                                                                         caAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTg                                                                           <  1:393351/149‑1 (MQ=255)
                                                                          aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                          >  1:8829/1‑149 (MQ=255)
                                                                          aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                          >  2:1083058/1‑149 (MQ=255)
                                                                              cAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCACAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACt                                                                      >  1:536189/1‑149 (MQ=255)
                                                                                   aaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAgc                                                                 >  1:804071/1‑149 (MQ=255)
                                                                                   aaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAgc                                                                 >  1:804073/1‑149 (MQ=255)
                                                                                     cagggagCCACTATGTGGCGTCGGATGGATGCGCTGGCAAAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAACGAACGAAGATGATTGTATTATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAgagc                                                               <  2:706449/143‑1 (MQ=255)
                                                                                      ccGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGcgcc                                                              <  1:461437/149‑1 (MQ=255)
                                                                                        gCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGcgccgc                                                            >  2:570060/1‑149 (MQ=255)
                                                                                                 aTGTGGCGTCGTATGGATGTCCTGTATAATATGATAGATGCCTGCTTCCCCCATGGAAGAGGGTTAAATAAGGAAAGAAGATGATTCAACAAATTATCAGGCCCGACTCTGTCGAACACGCGCTGGAACTGAAGCGCCGCTACCAGGAt                                                   <  2:640313/149‑1 (MQ=255)
                                                                                                     ggCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAg                                               <  1:973941/149‑1 (MQ=255)
                                                                                                      gCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGc                                              >  2:485933/1‑149 (MQ=255)
                                                                                                          cGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTc                                          <  1:655498/149‑1 (MQ=255)
                                                                                                          cGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTc                                          <  1:62085/149‑1 (MQ=255)
                                                                                                           ggatggatGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGtat                                         >  1:495997/1‑147 (MQ=255)
                                                                                                             atggatgCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTgg                                       <  1:622932/149‑1 (MQ=255)
                                                                                                                 atgCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCg                                   >  1:887278/1‑149 (MQ=255)
                                                                                                                    cgcTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCg                                <  2:465005/149‑1 (MQ=255)
                                                                                                                      cTGGGATAAATGTCAGATGCCGTCTTCCCGCAGGGAAGTGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGATCTGTAGCGCCGCTAGCAGGTTGAAGCCGGCTGGTTCGCCggg                              <  2:299595/149‑1 (MQ=255)
                                                                                                                          cATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAAAAAGGAAAGAAGATGATAGAACAATTATTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgca                          <  1:1081744/149‑1 (MQ=255)
                                                                                                                                   aCAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTc                  <  2:216030/149‑1 (MQ=255)
                                                                                                                                   aCAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTCTCGAACAGGAGCTGGAACTGAAGCGCCGCTACCAGGATGACGCCGTCTGGTTCGCCGGGGGCACCTAACTc                  >  2:178127/1‑149 (MQ=255)
                                                                                                                                    cAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCa                 <  1:570060/149‑1 (MQ=255)
                                                                                                                                       ctGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATg              <  1:1045979/148‑1 (MQ=255)
                                                                                                                                       aTGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATg              >  1:1047720/1‑149 (MQ=255)
                                                                                                                                        tGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGc             >  1:290524/1‑149 (MQ=255)
                                                                                                                                         gCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCt            >  1:203741/1‑149 (MQ=255)
                                                                                                                                          cccTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAGTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGAGCCGCTCCCAGGAAGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGGTa           >  1:482792/1‑149 (MQ=255)
                                                                                                                                          cccTCTTCCCACAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTa           <  2:536189/149‑1 (MQ=255)
                                                                                                                                           ccTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTac          <  2:290524/149‑1 (MQ=255)
                                                                                                                                            ctctTCCCTGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaca         >  1:466844/1‑149 (MQ=255)
                                                                                                                                            ctctTCCCGAAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTACAGGCCCGACTCAGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGACTGGTTCGCCCGGGGCAGCAAACACAATGCTaca         >  2:786109/1‑149 (MQ=255)
                                                                                                                                             tctTCCCTGCGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTaccc        >  2:61762/1‑147 (MQ=255)
                                                                                                                                              ctTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACAcc       >  1:1167615/1‑149 (MQ=255)
                                                                                                                                                  ccGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAAcc   >  1:412646/1‑149 (MQ=255)
                                                                                                                                                  ccGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAATTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTAGTTCGCCTGGGGCAGCAAACTCAATGCTACACCAccc   >  2:370912/1‑149 (MQ=255)
                                                                                                                                                   cGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAAccc  >  1:635364/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAAATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCC  >  NC_000913/3020363‑3020657

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: