Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 188,226 T→C 4.7% intergenic (+31/‑116) tsf → / → pyrH translation elongation factor EF‑Ts/uridylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273188,2260TC 4.7% 247.8 / 2.9 86intergenic (+31/‑116)tsf/pyrHtranslation elongation factor EF‑Ts/uridylate kinase
Reads supporting (aligned to +/- strand):  ref base T (48/33);  new base C (2/2);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

CTGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTC  >  CP009273/188079‑188365
                                                                                                                                                    |                                                                                                                                           
ctgAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGcttc                                                                                                                                            >  1:182686/1‑149 (MQ=255)
   aaaGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCttt                                                                                                                                         >  1:541720/1‑149 (MQ=255)
    aaGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCGCCCTCAGGCGGCTCCtttt                                                                                                                                        <  1:1094476/149‑1 (MQ=255)
     agagCATAACGCTGAAGTGACTGGCTCCACCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTCGAGCCTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCttttt                                                                                                                                       <  1:51626/149‑1 (MQ=255)
     agagCAAAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCttttt                                                                                                                                       <  2:246194/149‑1 (MQ=255)
      gagCATTACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTtg                                                                                                                                      <  1:20170/149‑1 (MQ=255)
        gCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTtgtg                                                                                                                                    >  2:98248/1‑149 (MQ=255)
         cATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGc                                                                                                                                   <  1:618095/149‑1 (MQ=255)
         cATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCGCCCTCAGGCGGCTCCTTTTtgata                                                                                                                                   <  2:454996/149‑4 (MQ=38)
           tAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAGGGGATTGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGccc                                                                                                                                 >  1:1173720/1‑149 (MQ=255)
              cGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGGCTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTGAATTATCACAAAGGAGCCGCCTGAGGGCGCCTTCTTTTTGTGCCCATc                                                                                                                              >  1:373674/1‑149 (MQ=255)
                 tGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTg                                                                                                                           <  1:51022/149‑1 (MQ=255)
                    cgTGACTGGTTTCAGGCGCGTTGAAGGGGGTGAACGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTaa                                                                                                                        <  1:1188003/148‑1 (MQ=255)
                    aGTGACTGGCTTCATCCGCTTCCAAGTGGGTGAAGGCATCGATAAAGCTGAGACTGACTCTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGACTTCTTTTTGAGCCCATCTTGTaa                                                                                                                        <  2:785821/149‑1 (MQ=255)
                       gACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAAtt                                                                                                                     >  1:174943/1‑149 (MQ=255)
                        aCTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTc                                                                                                                    <  1:1141491/149‑1 (MQ=255)
                        aCTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTc                                                                                                                    >  2:206037/1‑149 (MQ=255)
                          tGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCACGCAGTGTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAg                                                                                                                  <  2:994125/149‑1 (MQ=255)
                          tGGCGTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAg                                                                                                                  <  1:1202711/149‑1 (MQ=255)
                           ggCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCGATCTTGTAAATTCAGc                                                                                                                 <  2:216498/149‑1 (MQ=255)
                           ggCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGc                                                                                                                 <  1:875007/149‑1 (MQ=255)
                            gCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCt                                                                                                                >  2:1003972/1‑149 (MQ=255)
                                  tCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGGTGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTGTTTGTGCCCATCTTGTAAATTCAGCTAACCCt                                                                                                          >  2:515688/1‑149 (MQ=255)
                                  tCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCt                                                                                                          >  1:1017567/1‑149 (MQ=255)
                                     gCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTtgt                                                                                                       <  1:448028/149‑1 (MQ=255)
                                      ctacGAAGGGGGTGAAGGCATCGAGAAAGTTGAGACTGTTTATGCATGAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTtgtg                                                                                                      <  2:871657/146‑1 (MQ=255)
                                      cTTCGAAGTGGGTGAAGGCATCGGGAAAGTTGAGAGTGACTTTGCAGCAGAAGTTGCTGCGAAGTCCAAGCAGTCTGAATTATCAAAAAGGTGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTtgtg                                                                                                      <  1:426108/149‑1 (MQ=255)
                                      cTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTACATTCAGCTAACCCTtgtg                                                                                                      >  2:460085/1‑149 (MQ=255)
                                      cTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCAGCTTGTAAATTCAGCTAACCCTggtg                                                                                                      >  2:277262/1‑149 (MQ=255)
                                                  tGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGa                                                                                          >  2:239011/1‑149 (MQ=255)
                                                  tGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGa                                                                                          >  1:1129213/1‑149 (MQ=255)
                                                      ggCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTACGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAg                                                                                      <  1:1019819/149‑1 (MQ=255)
                                                         aTCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGa                                                                                   >  1:1122570/1‑149 (MQ=255)
                                                          tCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCATGCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGAc                                                                                  >  2:421675/1‑149 (MQ=255)
                                                          tCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGAc                                                                                  <  1:436799/149‑1 (MQ=255)
                                                          tCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTCCGCTGAAAAGCGAc                                                                                  >  2:168959/1‑149 (MQ=255)
                                                          tCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGGCCAAGAAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGAc                                                                                  <  1:394119/149‑1 (MQ=255)
                                                           cGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTCTCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAACTCAGCTACCCCTTGTGGGGCTGCGCTGAAAAGCGACg                                                                                 <  1:803843/149‑1 (MQ=255)
                                                             agaAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTa                                                                               >  1:656970/1‑149 (MQ=255)
                                                             agaAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTa                                                                               >  1:794385/1‑149 (MQ=255)
                                                              gaAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTAc                                                                              >  1:494606/1‑149 (MQ=255)
                                                               aaaGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACa                                                                             >  1:501279/1‑149 (MQ=255)
                                                               aaaGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTGAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTGCTGCCCAGCCTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACa                                                                             >  2:223875/1‑149 (MQ=255)
                                                                     gagaCTGACTTGGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAACGTCg                                                                       <  1:192611/149‑1 (MQ=255)
                                                                       gactgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCt                                                                     >  1:460722/1‑149 (MQ=255)
                                                                           gactTTGCAGCAGAAGTTGCTGCGATGTCCAGGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTa                                                                 >  2:1001936/1‑149 (MQ=255)
                                                                                 gcagcaGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAAtt                                                           >  1:173906/1‑149 (MQ=255)
                                                                                 gcagcaGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCATGTGGGGCTGCGCTGAAAAGCGACGTACAATGTTGCTAGTATTAAtt                                                           >  1:916062/1‑149 (MQ=255)
                                                                                 gcagcaGAAGTTGCTGCGATGTCCAAGCAGGCTTAATGATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTGTTAAtt                                                           >  2:113055/1‑149 (MQ=255)
                                                                                    gcagAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCAt                                                        >  2:1008709/1‑149 (MQ=255)
                                                                                    gcagAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCAt                                                        >  1:771537/1‑149 (MQ=255)
                                                                                    gcagAAGTTGCGGCGATGTCCAAGCAGTCTTAATTAGCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCAt                                                        >  2:958281/1‑149 (MQ=255)
                                                                                      agAAGTTGCTGCCATGTCAAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTGTTGTGCCCATCTTGGAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAATAACGACGTACTATGTCGCTAGTATTAATTCAttt                                                      <  1:299498/149‑1 (MQ=255)
                                                                                          gTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAAt                                                  <  2:264197/149‑1 (MQ=255)
                                                                                              ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                              >  1:816413/1‑149 (MQ=255)
                                                                                              ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                              <  1:716941/149‑1 (MQ=255)
                                                                                                gctaTGACCAAGCAGTATTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGa                                            <  1:45615/146‑1 (MQ=255)
                                                                                                 cGATGTCCAAGCAGTCTTAATTATCAAAGAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGAc                                           >  2:729896/1‑149 (MQ=255)
                                                                                                 cGATGTCCAAGCAGTCTTAATTATCAAAATGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGGACAATGTCGCTAGTATTAATTCATTTCAATCGTTgcc                                           >  1:807223/1‑147 (MQ=255)
                                                                                                 cGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGAACAATGTCGCTAGTATTAATTCATTTCAATCGTTAAc                                           >  1:116827/1‑149 (MQ=255)
                                                                                                        cAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAATCTCa                                    >  1:93228/1‑149 (MQ=255)
                                                                                                        cAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATATTGTAAATTCAGCTAAACCTTTTGGGGCTGCGCTGAAAAGCGAGGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCa                                    >  1:1218591/1‑149 (MQ=255)
                                                                                                         aaGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCCATCGTTGACAGTCTCAg                                   >  1:448018/1‑149 (MQ=255)
                                                                                                         aaGCAGTCTTAATTATCAAAAAGGAGCCGCATGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAg                                   >  2:386898/1‑149 (MQ=255)
                                                                                                            ctgTCTTAATTATCAAAAAGGAGCCCCCTGAGGGCGGCGTCTCTTTGCGCCGATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaa                                <  2:261553/147‑1 (MQ=255)
                                                                                                            cAGTCTTAATTATCAAAAACGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaa                                <  2:26126/149‑1 (MQ=255)
                                                                                                             aGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaa                               <  1:981673/149‑1 (MQ=255)
                                                                                                                 ttAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATCAATTTCAATCGTTGACAGTCTCAGgcaaaaaa                           >  2:309125/1‑149 (MQ=255)
                                                                                                                   aaTTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATGAATCCATTTCAATCGTTGACAGTCTCAGGAAAGAAACa                         <  1:667979/149‑1 (MQ=255)
                                                                                                                      tATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATgg                      >  2:583168/1‑149 (MQ=255)
                                                                                                                       aTCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGc                     >  2:420374/1‑149 (MQ=255)
                                                                                                                        acAAAAAGAAGCCGCCCTCAGGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGGTGCGCCGAAAAGCGCAGTACAATGTCGCAAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCt                    >  1:562566/2‑149 (MQ=37)
                                                                                                                        acAAAAAGAAGCCGCCCTCAGGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCt                    >  1:1237133/2‑149 (MQ=38)
                                                                                                                         cAAAACGGAGCCGCCTGATGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTa                   <  1:792126/149‑1 (MQ=255)
                                                                                                                         cAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTa                   <  1:239011/149‑1 (MQ=255)
                                                                                                                          aaaaaGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACCGTCTCGGGAAACAAACCAGGCTAc                  >  2:188357/1‑149 (MQ=255)
                                                                                                                           aaaaGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAcc                 >  2:505736/1‑149 (MQ=255)
                                                                                                                             aaGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCaa               >  2:827700/1‑149 (MQ=255)
                                                                                                                             aaGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATGCATTCCAATCGTTGACAGTAGCAGGAAAGAAACATGGCAACCaa               >  1:124745/1‑149 (MQ=255)
                                                                                                                               ggAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATg             <  1:984081/149‑1 (MQ=255)
                                                                                                                                   ccgccTGAGGGCGGCTTCTTTTTGTGCCCATATTGAAAGCACTTCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCCATGCaaa         >  1:605800/1‑149 (MQ=255)
                                                                                                                                        cttaGGGGGGCTTCTTTTTGCGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCGTTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAAccc     <  2:451812/146‑1 (MQ=255)
                                                                                                                                        ccgAGGGCGTGTTCTTTTTGGGCCCATCTTGTAAATGCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTGGTTCCTTCCAATCTTTGACAGTCTCTGGAAAGAAACATTGCTACCAATGCAAAAccc     <  2:799464/147‑1 (MQ=255)
                                                                                                                                         tGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAATGAAACGTGGCTACCAATGCAAAACCCg    >  1:861922/1‑149 (MQ=255)
                                                                                                                                         tGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTCCCAATGCAAAACCCg    >  2:403476/1‑149 (MQ=255)
                                                                                                                                          gAGGGCGGCTGCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTGGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTACTTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGt   <  2:232036/149‑1 (MQ=255)
                                                                                                                                           aGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTAGTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAGCCCGTc  >  2:984504/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CTGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCGCCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTC  >  CP009273/188079‑188365

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