Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 319,024 C→A 91.4% R326S (CGC→AGC)  ykgF → ferridoxin‑like LutB family protein; putative electron transport chain YkgEFG component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273319,0240CA91.4% 499.5 / 40.1 187R326S (CGC→AGC) ykgFferridoxin‑like LutB family protein; putative electron transport chain YkgEFG component
Reads supporting (aligned to +/- strand):  ref base C (8/8);  new base A (75/96);  total (83/104)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 188 total aligned reads displayed.

GCGCGAAGCTGGGCACGTTGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGG  >  CP009273/318879‑319168
                                                                                                                                                 |                                                                                                                                                
gcgcgAAGCTGGGCACGTTGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGc                                                                                                                                                <  2:2649011/148‑1 (MQ=255)
   cgAAGCTGGGCACGTTGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCat                                                                                                                                             <  2:1285145/148‑1 (MQ=255)
    gAAGCTGGGCACGTTGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCata                                                                                                                                            <  2:143116/148‑1 (MQ=255)
      aGCTGGGCACGTTGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATAtt                                                                                                                                          <  1:2561178/148‑1 (MQ=255)
        cTGGGCACGTTGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCCTATTgg                                                                                                                                        >  2:2060803/1‑148 (MQ=255)
        cTGGGCACGTTGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTgg                                                                                                                                        >  2:1009536/1‑148 (MQ=255)
          gggCACGTTGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCg                                                                                                                                      >  1:1035523/1‑148 (MQ=255)
            gCACGTTGATGGTCCTGAAGATTTTCATCAGGTTATTGTCGATAACGGGCGTTCTGAGGTGCAGGCCTCTGAAATTCTTTTTGTGCTGAGCTGTATTAGCTGCGTGGTTTGTATGTATACTTGTCCGGAATATAGCAATACTGGCGGt                                                                                                                                    <  2:553314/148‑1 (MQ=255)
               cGTTGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCAt                                                                                                                                 >  1:1025821/1‑148 (MQ=255)
                 ttGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATgg                                                                                                                               <  1:622598/148‑1 (MQ=255)
                  tGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGa                                                                                                                              >  1:2958631/1‑148 (MQ=255)
                   gATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGat                                                                                                                             <  2:2332228/148‑1 (MQ=255)
                   gATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGat                                                                                                                             <  2:148853/148‑1 (MQ=255)
                     tGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGatat                                                                                                                           >  2:1127282/1‑148 (MQ=255)
                     tGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGatat                                                                                                                           >  2:2434135/1‑148 (MQ=255)
                     tGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGatat                                                                                                                           >  1:227464/1‑148 (MQ=255)
                       gTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATgg                                                                                                                         <  2:1888490/148‑1 (MQ=255)
                        tCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGc                                                                                                                        >  2:769936/1‑148 (MQ=255)
                        tCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGc                                                                                                                        >  2:2979576/1‑148 (MQ=255)
                           tGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGctct                                                                                                                     >  1:1784547/1‑148 (MQ=255)
                            gAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTa                                                                                                                    >  1:920152/1‑148 (MQ=255)
                              agagTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTAtt                                                                                                                  >  2:2817349/1‑148 (MQ=255)
                              agagTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTAtt                                                                                                                  >  2:471032/1‑148 (MQ=255)
                               gagTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTAttt                                                                                                                 <  1:2256541/148‑1 (MQ=255)
                               gagTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTAttt                                                                                                                 >  2:2786237/1‑148 (MQ=255)
                                agTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTa                                                                                                                >  2:3138466/1‑148 (MQ=255)
                                 gTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTAt                                                                                                               <  2:2156727/148‑1 (MQ=255)
                                 gTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTAt                                                                                                               <  2:2741033/148‑1 (MQ=255)
                                  tttCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATc                                                                                                              >  2:3049193/1‑148 (MQ=255)
                                     cATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAg                                                                                                           <  1:1394165/148‑1 (MQ=255)
                                       tCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAggg                                                                                                         >  1:1428586/1‑148 (MQ=255)
                                        cTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGc                                                                                                        <  2:2149041/148‑1 (MQ=255)
                                        cTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGc                                                                                                        <  2:1456493/148‑1 (MQ=255)
                                          ggTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCa                                                                                                      >  2:2589038/1‑148 (MQ=255)
                                                  tCGAAAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGc                                                                                              >  2:218326/1‑148 (MQ=255)
                                                    gATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCgg                                                                                            <  1:1400836/148‑1 (MQ=255)
                                                    gATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCgg                                                                                            <  1:1557982/148‑1 (MQ=255)
                                                      tAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTg                                                                                          >  2:702525/1‑148 (MQ=255)
                                                      tAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTg                                                                                          >  2:702632/1‑148 (MQ=255)
                                                            gCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTtct                                                                                    >  1:2672193/1‑148 (MQ=255)
                                                            gCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTtct                                                                                    >  2:2357994/1‑148 (MQ=255)
                                                             cGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTtctc                                                                                   <  1:832444/148‑1 (MQ=255)
                                                                tCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGc                                                                                >  2:2353400/1‑148 (MQ=255)
                                                                tCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGc                                                                                >  2:3097669/1‑148 (MQ=255)
                                                                tCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGc                                                                                >  2:1661455/1‑148 (MQ=255)
                                                                  tGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTa                                                                              <  2:337627/148‑1 (MQ=255)
                                                                  tGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTa                                                                              <  2:336858/148‑1 (MQ=255)
                                                                     ggTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACtt                                                                           <  2:1928147/148‑1 (MQ=255)
                                                                      gTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTg                                                                          >  1:3060606/1‑148 (MQ=255)
                                                                            ggctctGAATTTCTGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCTTATTGGCGGTCTTGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTggcggat                                                                    <  1:3138697/146‑3 (MQ=255)
                                                                            gCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCt                                                                    >  2:2639131/1‑148 (MQ=255)
                                                                            gCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCt                                                                    <  1:3138620/148‑1 (MQ=255)
                                                                             ccTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCta                                                                   >  1:2985100/1‑148 (MQ=255)
                                                                             ccTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCta                                                                   <  1:1624078/148‑1 (MQ=255)
                                                                                ctGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaa                                                                <  2:2518605/148‑1 (MQ=255)
                                                                                ctGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaa                                                                <  1:1885157/148‑1 (MQ=255)
                                                                                  gAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaag                                                              <  1:521829/148‑1 (MQ=255)
                                                                                   aaTTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaaga                                                             >  2:830291/1‑148 (MQ=255)
                                                                                   aaTTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaaga                                                             <  2:1891573/148‑1 (MQ=255)
                                                                                    aTTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagat                                                            >  1:678959/1‑148 (MQ=255)
                                                                                     tttCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatt                                                           <  2:1714236/148‑1 (MQ=255)
                                                                                     tttCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatt                                                           <  2:39601/148‑1 (MQ=255)
                                                                                        cGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagatttta                                                        >  2:1985950/1‑148 (MQ=255)
                                                                                          ggATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaa                                                      <  2:35404/148‑1 (MQ=255)
                                                                                          ggATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaa                                                      >  2:1170829/1‑148 (MQ=255)
                                                                                           gATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaag                                                     <  1:1201525/148‑1 (MQ=255)
                                                                                           gATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaag                                                     <  2:2722739/148‑1 (MQ=255)
                                                                                            aagtgCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGTCCAATTGGTGAGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaaga                                                    >  2:1123091/3‑148 (MQ=255)
                                                                                            aTGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaaga                                                    <  1:2483452/148‑1 (MQ=255)
                                                                                             tgtgCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagat                                                   <  2:1341189/148‑1 (MQ=255)
                                                                                             tgtgCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagat                                                   <  2:1179714/148‑1 (MQ=255)
                                                                                             tgtgCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTAATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagat                                                   <  1:1528174/148‑1 (MQ=255)
                                                                                              gtgCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                  >  2:631548/1‑148 (MQ=255)
                                                                                              gtgCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTAtaaagattttaaagatt                                                  >  2:319742/1‑148 (MQ=255)
                                                                                                gctgcGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTa                                                <  2:2234011/148‑1 (MQ=255)
                                                                                                      cTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTAc                                          >  1:972436/1‑148 (MQ=255)
                                                                                                      cTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTAc                                          >  1:444184/1‑148 (MQ=255)
                                                                                                      cTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTAc                                          >  1:2305802/1‑148 (MQ=255)
                                                                                                       tGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACg                                         >  2:2257112/1‑148 (MQ=255)
                                                                                                        gTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGc                                        <  1:214707/148‑1 (MQ=255)
                                                                                                        gTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGc                                        >  1:1548497/1‑148 (MQ=255)
                                                                                                         tATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGcc                                       <  2:3111102/148‑1 (MQ=255)
                                                                                                           ttCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTg                                     <  1:1953982/148‑1 (MQ=255)
                                                                                                            tCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGc                                    <  2:2684686/148‑1 (MQ=255)
                                                                                                                tGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTGTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCtt                                <  1:1591452/148‑1 (MQ=255)
                                                                                                                tGCGGGGCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCtt                                <  2:423784/148‑1 (MQ=255)
                                                                                                                   ggggCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTAt                             <  1:2235688/148‑1 (MQ=255)
                                                                                                                   ggggCTTGTAAGCATTCTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTAt                             <  2:1836287/148‑1 (MQ=255)
                                                                                                                     ggCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGc                           >  1:746281/1‑148 (MQ=255)
                                                                                                                     ggCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGc                           >  1:378357/1‑148 (MQ=255)
                                                                                                                     ggCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGc                           >  1:1961240/1‑148 (MQ=255)
                                                                                                                      gCTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGca                          <  2:769252/148‑1 (MQ=255)
                                                                                                                       cTTGTATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGcac                         <  2:2960065/148‑1 (MQ=255)
                                                                                                                           tATGAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCt                     <  1:2786237/148‑1 (MQ=255)
                                                                                                                              gAATACTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTtgt                  >  1:1666362/1‑148 (MQ=255)
                                                                                                                                  aCTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACa              >  1:2652318/1‑148 (MQ=255)
                                                                                                                                   cTTGTCCGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACaa             >  1:2064601/1‑148 (MQ=255)
                                                                                                                                        ccGGCATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACgtgt        >  2:338467/1‑148 (MQ=255)
                                                                                                                                            cATATAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTcc    >  1:1380681/1‑148 (MQ=255)
                                                                                                                                             atatAGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCg   <  2:1051099/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCGCGAAGCTGGGCACGTTGATGGTCCTGAAGAGTTTCATCTGGTTATTGTCGATAACGGGCGTTCTGAGGTGCTGGCCTCTGAATTTCGGGATGTGCTGCGCTGTATTCGCTGCGGGGCTTGTATGAATACTTGTCCGGCATATCGCCATATTGGCGGTCATGGATATGGCTCTATTTATCCAGGGCCAATTGGTGCGGTGATTTCTCCGCTACTTGGCGGCTATAAAGATTTTAAAGATTTACCCTACGCCTGCTCTTTATGCACAGCTTGTGACAACGTGTGTCCGG  >  CP009273/318879‑319168

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: ‑