Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,605,542 T→A 6.6% G447G (GGA→GGT yjjI ← DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,605,5420TA 6.6% 561.2 / 24.1 182G447G (GGA→GGTyjjIDUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand):  ref base T (97/73);  new base A (4/8);  total (101/81)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 188 total aligned reads displayed.

AGTGTTGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGC  >  CP009273/4605395‑4605689
                                                                                                                                                   |                                                                                                                                                   
     tgcgtgcgGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGc                                                                                                                                                <  2:509725/148‑1 (MQ=255)
      gcgtgcgGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCt                                                                                                                                               >  1:1073102/1‑148 (MQ=255)
      gcgtgcgGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCt                                                                                                                                               <  1:3340856/148‑1 (MQ=255)
      gcgtgcgGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCt                                                                                                                                               >  2:2342432/1‑148 (MQ=255)
        gtgcgGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCttt                                                                                                                                             <  2:1151411/148‑1 (MQ=255)
           cgGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTaaa                                                                                                                                          >  1:114583/1‑148 (MQ=255)
              cttcttCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGc                                                                                                                                       >  1:384839/1‑148 (MQ=255)
               ttcttcGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCa                                                                                                                                      <  2:3354487/148‑1 (MQ=255)
               ttcttcGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCa                                                                                                                                      <  2:1800786/148‑1 (MQ=255)
               ttcttcGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCa                                                                                                                                      <  1:1837198/148‑1 (MQ=255)
                      cccAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGc                                                                                                                               >  2:2229308/1‑148 (MQ=255)
                       ccagccagGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCa                                                                                                                              <  1:290866/148‑1 (MQ=255)
                          gccagGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAg                                                                                                                           <  2:3124375/148‑1 (MQ=255)
                           ccagGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGc                                                                                                                          <  2:511767/148‑1 (MQ=255)
                           ccagGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGc                                                                                                                          >  1:2106534/1‑148 (MQ=255)
                             agGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTg                                                                                                                        >  2:2381798/1‑148 (MQ=255)
                               gtggtgTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTa                                                                                                                      <  2:449152/148‑1 (MQ=255)
                                tggtgTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTAc                                                                                                                     >  1:2249047/1‑148 (MQ=255)
                                        gTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCt                                                                                                             >  2:2600055/1‑148 (MQ=255)
                                          gCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGc                                                                                                           >  2:2900925/1‑148 (MQ=255)
                                           cGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCg                                                                                                          <  1:124555/148‑1 (MQ=255)
                                              gAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGt                                                                                                       <  2:16788/148‑1 (MQ=255)
                                              gAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGt                                                                                                       >  1:2204014/1‑148 (MQ=255)
                                              gAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGt                                                                                                       <  2:16920/148‑1 (MQ=255)
                                               aaCCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGtt                                                                                                      >  1:1199903/1‑148 (MQ=255)
                                                    tCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGtt                                                                                                 >  1:2002956/1‑148 (MQ=255)
                                                        cgcgATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGAt                                                                                             >  2:3536881/1‑148 (MQ=255)
                                                        cgcgATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGAt                                                                                             <  1:1152605/148‑1 (MQ=255)
                                                         gcgATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATg                                                                                            <  1:2320922/148‑1 (MQ=255)
                                                            atatTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGtt                                                                                         <  2:558301/148‑1 (MQ=255)
                                                            atatTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGtt                                                                                         <  1:1120149/148‑1 (MQ=255)
                                                            atatTTTTCTAAATCCGAAAAACGCACCATATAACCGAAAACGCGAAACAGATCGTTACCGCTGACATTGTCGGTAAATTCAAGCATTCCGGCATAAAAGGCACCGAGGCAAAGCTGTCCCAGTGCCTGCGGGTTAGGTATGATGGta                                                                                         >  1:1580412/1‑147 (MQ=255)
                                                             tatTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGttt                                                                                        <  2:1840252/148‑1 (MQ=255)
                                                              atTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTtc                                                                                       >  1:1062638/1‑148 (MQ=255)
                                                                ttttCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTtcgt                                                                                     <  1:177575/148‑1 (MQ=255)
                                                                  ttCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTtcgtcg                                                                                   <  2:2315984/148‑1 (MQ=255)
                                                                    cTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAg                                                                                 >  1:1951565/1‑148 (MQ=255)
                                                                    cTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAg                                                                                 >  1:421710/1‑148 (MQ=255)
                                                                    cTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAg                                                                                 <  1:1244123/148‑1 (MQ=255)
                                                                     tAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGc                                                                                >  1:1263742/1‑148 (MQ=255)
                                                                       aaTCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGt                                                                              <  1:3340674/148‑1 (MQ=255)
                                                                       aaTCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGt                                                                              <  2:2861552/148‑1 (MQ=255)
                                                                       aaTCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGt                                                                              <  2:1086720/148‑1 (MQ=255)
                                                                          ccGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAg                                                                           >  1:61042/1‑148 (MQ=255)
                                                                          ccGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGGTTGATGGTTTCGTCGAGCGTCAg                                                                           >  2:2663225/1‑148 (MQ=255)
                                                                           cGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGa                                                                          <  2:632284/148‑1 (MQ=255)
                                                                                aaCGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTc                                                                     >  2:1030373/1‑148 (MQ=255)
                                                                                 aCGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCg                                                                    <  1:3123546/148‑1 (MQ=255)
                                                                                 aCGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCg                                                                    <  1:2843705/148‑1 (MQ=255)
                                                                                  cGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGACTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGc                                                                   >  2:678544/1‑148 (MQ=255)
                                                                                    cACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTg                                                                 <  1:2210525/148‑1 (MQ=255)
                                                                                    cACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTg                                                                 <  1:1903320/148‑1 (MQ=255)
                                                                                      ccATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGAt                                                               >  2:2674260/1‑148 (MQ=255)
                                                                                      ccATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGAt                                                               >  1:3011511/1‑148 (MQ=255)
                                                                                       cATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATg                                                              <  2:2488391/148‑1 (MQ=255)
                                                                                         tataACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGcc                                                            >  1:3117147/1‑148 (MQ=255)
                                                                                         tataACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGcc                                                            >  2:2558272/1‑148 (MQ=255)
                                                                                          ataACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCg                                                           >  2:2006673/1‑148 (MQ=255)
                                                                                           taACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCgg                                                          >  2:3547499/1‑148 (MQ=255)
                                                                                           taACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCgg                                                          >  2:1521881/1‑148 (MQ=255)
                                                                                            aaCCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGa                                                         >  1:1189065/1‑148 (MQ=255)
                                                                                            aaCCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGa                                                         >  1:2890258/1‑148 (MQ=255)
                                                                                               cGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaata                                                      <  1:2495974/148‑1 (MQ=255)
                                                                                               cGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATACCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaata                                                      <  1:2510590/148‑1 (MQ=255)
                                                                                               cGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATACCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaata                                                      <  1:2248854/148‑1 (MQ=255)
                                                                                                ggTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaataa                                                     <  1:1669967/148‑1 (MQ=255)
                                                                                                 gTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaataat                                                    <  2:1950552/148‑1 (MQ=255)
                                                                                                 gTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATACCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaataat                                                    <  2:2989709/148‑1 (MQ=255)
                                                                                                    aCGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGaataataat                                                 >  1:2831437/1‑148 (MQ=255)
                                                                                                     cgcgAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAataataata                                                <  1:987784/148‑1 (MQ=255)
                                                                                                     cgcgAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAataataata                                                <  1:1465112/148‑1 (MQ=255)
                                                                                                     cgcgAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAataataata                                                >  1:1303510/1‑148 (MQ=255)
                                                                                                     cgcgAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAataataata                                                >  1:1101914/1‑148 (MQ=255)
                                                                                                      gcgAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAAtaataataa                                               >  1:3374039/1‑148 (MQ=255)
                                                                                                        gAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGc                                             >  1:871327/1‑148 (MQ=255)
                                                                                                         aaCCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCa                                            >  1:2468602/1‑148 (MQ=255)
                                                                                                           ccAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCatg                                          <  1:1521881/148‑1 (MQ=255)
                                                                                                           ccAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCatg                                          <  2:13364/148‑1 (MQ=255)
                                                                                                           ccAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCatg                                          >  2:536612/1‑148 (MQ=255)
                                                                                                            cAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATACCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGAGGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCatga                                         >  2:1895871/1‑148 (MQ=255)
                                                                                                               aTCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGc                                      >  2:980195/1‑148 (MQ=255)
                                                                                                               aTCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGc                                      >  1:2024775/1‑148 (MQ=255)
                                                                                                               aTCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGc                                      >  2:583399/1‑148 (MQ=255)
                                                                                                                tCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCg                                     >  2:897708/1‑148 (MQ=255)
                                                                                                                  gTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCggg                                   <  1:3491697/148‑1 (MQ=255)
                                                                                                                    tACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGc                                 >  1:3342503/1‑148 (MQ=255)
                                                                                                                       cGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGAc                              <  2:531711/148‑1 (MQ=255)
                                                                                                                           gACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGtt                          >  1:3169543/1‑148 (MQ=255)
                                                                                                                            aCATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGttt                         >  2:483477/1‑148 (MQ=255)
                                                                                                                               ttGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCa                      >  2:1233720/1‑148 (MQ=255)
                                                                                                                                tGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAg                     <  2:2106534/148‑1 (MQ=255)
                                                                                                                                tGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAg                     >  1:589595/1‑148 (MQ=255)
                                                                                                                                    ggTAAATTCACGCATACCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATgg                 <  2:1347410/148‑1 (MQ=255)
                                                                                                                                    ggTAAATTCACGCATACCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATgg                 <  2:1347369/148‑1 (MQ=255)
                                                                                                                                         aTTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGAt            <  2:114583/148‑1 (MQ=255)
                                                                                                                                         aTTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGAt            >  1:2027722/1‑148 (MQ=255)
                                                                                                                                         aTTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGAt            <  2:384883/148‑1 (MQ=255)
                                                                                                                                            cACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCgg         <  1:2600055/148‑1 (MQ=255)
                                                                                                                                                  ttCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTgg   >  2:3416471/1‑148 (MQ=255)
                                                                                                                                                   tCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGc  >  2:1593152/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGTGTTGCGTGCGGCTTCTTCGCCCAGCCAGGTGGTGTTGGTGCGTGAACCTTCGGCGCGATATTTTTCTAAATCCGACAAACGCACCATATAACCGGTAACGCGAACCAGATCGTTACCGCTGACATTGGCGGTAAATTCACGCATTCCGGCTTTAAAGGCACCGAGGCAAAGCTGTACCAGTGCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGC  >  CP009273/4605395‑4605689

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