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

CAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCC  >  NC_000913/3020363‑3020657
                                                                                                                                                     |                                                                                                                                                  
caGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGc                                                                                                                                                    >  1:305547/1‑149 (MQ=255)
caGCGTGATGCTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGc                                                                                                                                                    >  2:1045979/1‑149 (MQ=255)
 aGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCa                                                                                                                                                   >  1:537411/1‑149 (MQ=255)
  gCGTGATGGTCAATGGTGTGATGGGCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAg                                                                                                                                                  <  1:1094659/149‑1 (MQ=255)
   cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAgg                                                                                                                                                 <  2:845931/149‑1 (MQ=255)
   cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAgg                                                                                                                                                 <  1:408855/149‑1 (MQ=255)
   cGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAgg                                                                                                                                                 <  2:215422/149‑1 (MQ=255)
     tGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGa                                                                                                                                               <  1:446036/149‑1 (MQ=255)
      gATTGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTATAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGaa                                                                                                                                              <  2:1067510/149‑1 (MQ=255)
      gATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGCATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTACCCGCAGGGaa                                                                                                                                              >  1:176322/1‑149 (MQ=255)
         ggTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAagag                                                                                                                                           >  1:1130825/1‑149 (MQ=255)
          gTCAATGGTGTGATGGTCTATGAAGACCGGCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAgg                                                                                                                                          >  2:1087496/1‑149 (MQ=255)
           tCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCCCCGGGATGAggt                                                                                                                                         >  2:281898/1‑148 (MQ=255)
            cAATGGTGTGATTGTCTATGAAGACCGTCAGTATAACTTCTATTGCGAATCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGGCGGATGGATGCGCTGTCATAA‑ATGACAGATGCCCTCTTCCCTGCGGGAAGAgggc                                                                                                                                        <  2:334498/149‑2 (MQ=255)
            cAATGGTGTGATGGGCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCTGCGGGAAGAgggc                                                                                                                                        <  1:871499/149‑2 (MQ=255)
              aTGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTCCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGACAGTATGTGGCGTCGGATCGATGCGCTGGCATAA‑ATGACAGATGCCCTCTACCCGGAGGGAAGAGGGTTa                                                                                                                                      >  1:640313/1‑149 (MQ=255)
                      aTGGTCTATGAAGACCGTCAGATAAACTTCGATTGCTATTCCAATAATGAACAATCCAGAAAAGAAGCTGCCAGTATCTGGCGAAGTATGGATGACATGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAATGa                                                                                                                              >  2:828261/1‑149 (MQ=255)
                           cTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAaagaa                                                                                                                         >  1:216030/1‑149 (MQ=255)
                              tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGATGAGGGTTAAATAAGGAAAGAAgat                                                                                                                      >  2:655498/1‑149 (MQ=255)
                              tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                                                                                                                      >  2:183734/1‑149 (MQ=255)
                              tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                                                                                                                      >  1:429227/1‑149 (MQ=255)
                              tGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgat                                                                                                                      <  2:229417/149‑1 (MQ=255)
                                aaGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga                                                                                                                    <  1:631368/149‑1 (MQ=255)
                                aaGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatga                                                                                                                    >  1:41882/1‑149 (MQ=255)
                                 aGACCGTCAGGTTAACTACGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAgatgat                                                                                                                   >  1:299595/1‑149 (MQ=255)
                                   aCCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTg                                                                                                                 <  1:576409/149‑1 (MQ=255)
                                    ccGTCAGTTTAAATTCGATTGCTATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGa                                                                                                                <  2:943999/149‑1 (MQ=255)
                                     cGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCTGTATGTGGCGTAGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGATAGAAGATGATTGaa                                                                                                               >  2:403516/1‑149 (MQ=255)
                                     cGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGaa                                                                                                               >  1:442166/1‑149 (MQ=255)
                                       tCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACa                                                                                                             <  1:520536/149‑1 (MQ=255)
                                           tttAACTTCGATTGCGAGTCCTTTTATGCAGAAGCCAGATAAGCCGCTGCCAGTATGGGGCGTCGGATGGATGCGCTGGCATAA‑TTGACAGATGCCGTCTTCCCGCAGGGTTGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttt                                                                                                         <  1:160978/149‑1 (MQ=255)
                                             taccTTCGATTGCGATTCCATTTATGCCCAAGCCAGAAACGCCGCTGCAAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCAGCAGGGAAGAGGGTTGAATAAGGAAAGAATATGATTGAACAAttttt                                                                                                       <  2:809051/146‑1 (MQ=255)
                                             tAACTTCGATTGCGATTCCATTTATGCACAAGCCTGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAAttttt                                                                                                       <  1:1043014/149‑1 (MQ=255)
                                              aaCTTCGTTTGCTATTCCATTTATGCACAAGCCAGACAAGCCGCTGCCAGGATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTc                                                                                                      <  2:144784/149‑1 (MQ=255)
                                                cTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAg                                                                                                    <  1:127630/149‑1 (MQ=255)
                                                  tCGATTGCGATTCCATTTATTCACAAGCCAGAATAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTTGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTTAACAATTTTTCAGGc                                                                                                  <  2:569337/149‑1 (MQ=255)
                                                      ttGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAAGAAGGAAAGAAGATGATTGACCAATTTTTCAGGCCCGa                                                                                              >  2:1141281/1‑149 (MQ=255)
                                                            ttCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGt                                                                                        <  1:124587/149‑1 (MQ=255)
                                                             tCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTc                                                                                       <  1:358635/149‑1 (MQ=255)
                                                                aTTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGaa                                                                                    <  2:533645/149‑1 (MQ=255)
                                                                   tATGAACATGACATAAAAGCCGCTTCCAGTATGTGGCGGCGGTTGGTTGTGCTGGCATTA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGTGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCCGGCGCGACTCTGTCGAACAg                                                                                 <  2:1111303/149‑1 (MQ=255)
                                                                       cacaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAATTAAGGAAAGAAGATGTTTGAACAATTTTTCAGGCCCGACTCTGTCGAGCCGgcgc                                                                             >  1:563892/1‑149 (MQ=255)
                                                                        acaAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCt                                                                            >  2:637188/1‑149 (MQ=255)
                                                                         caAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTg                                                                           <  1:393351/149‑1 (MQ=255)
                                                                          aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                          >  1:8829/1‑149 (MQ=255)
                                                                          aaGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTgg                                                                          >  2:1083058/1‑149 (MQ=255)
                                                                              cAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCACAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACt                                                                      >  1:536189/1‑149 (MQ=255)
                                                                                   aaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAgc                                                                 >  1:804071/1‑149 (MQ=255)
                                                                                   aaGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAgc                                                                 >  1:804073/1‑149 (MQ=255)
                                                                                     cagggagCCACTATGTGGCGTCGGATGGATGCGCTGGCAAAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAACGAACGAAGATGATTGTATTATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAgagc                                                               <  2:706449/143‑1 (MQ=255)
                                                                                      ccGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGcgcc                                                              <  1:461437/149‑1 (MQ=255)
                                                                                        gCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGcgccgc                                                            >  2:570060/1‑149 (MQ=255)
                                                                                                 aTGTGGCGTCGTATGGATGTCCT‑GTATAATATGATAGATGCCTGCTTCCCCCATGGAAGAGGGTTAAATAAGGAAAGAAGATGATTCAACAAATTATCAGGCCCGACTCTGTCGAACACGCGCTGGAACTGAAGCGCCGCTACCAGGAt                                                   <  2:640313/149‑1 (MQ=255)
                                                                                                     ggCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAg                                               <  1:973941/149‑1 (MQ=255)
                                                                                                      gCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGc                                              >  2:485933/1‑149 (MQ=255)
                                                                                                          cGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTc                                          <  1:655498/149‑1 (MQ=255)
                                                                                                          cGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTc                                          <  1:62085/149‑1 (MQ=255)
                                                                                                           ggatggatGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGtat                                         >  1:495997/1‑147 (MQ=255)
                                                                                                             atggatgCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTgg                                       <  1:622932/149‑1 (MQ=255)
                                                                                                                 atgCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCg                                   >  1:887278/1‑149 (MQ=255)
                                                                                                                    cgcTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCg                                <  2:465005/149‑1 (MQ=255)
                                                                                                                      cTGGGATAA‑ATGTCAGATGCCGTCTTCCCGCAGGGAAGTGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGATCTGTAGCGCCGCTAGCAGGTTGAAGCCGGCTGGTTCGCCggg                              <  2:299595/149‑1 (MQ=255)
                                                                                                                          cATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAAAAAGGAAAGAAGATGATAGAACAATTATTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGgca                          <  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)
                                                                                                                                                     |                                                                                                                                                  
CAGCGTGATGGTCAATGGTGTGATGGTCTATGAAGACCGTCAGTTTAACTTCGATTGCGATTCCATTTATGCACAAGCCAGAAAAGCCGCTGCCAGTATGTGGCGTCGGATGGATGCGCTGGCATAA‑ATGACAGATGCCCTCTTCCCGCAGGGAAGAGGGTTAAATAAGGAAAGAAGATGATTGAACAATTTTTCAGGCCCGACTCTGTCGAACAGGCGCTGGAACTGAAGCGCCGCTACCAGGATGAAGCCGTCTGGTTCGCCGGGGGCAGCAAACTCAATGCTACACCAACCC  >  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: ‑