Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,020,049 C→G 8.2% P55A (CCG→GCG)  fre → NAD(P)H‑flavin reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,020,0490CG 8.2% 759.8 / 46.4 205P55A (CCG→GCG) freNAD(P)H‑flavin reductase
Reads supporting (aligned to +/- strand):  ref base C (85/103);  new base G (6/11);  total (91/114)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

GTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAG  >  CP009273/4019903‑4020195
                                                                                                                                                  |                                                                                                                                                  
gTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGcc                                                                                                                                                   <  1:779193/148‑1 (MQ=255)
  aaaGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCgg                                                                                                                                                 <  2:862426/148‑1 (MQ=255)
    aGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGAt                                                                                                                                               >  1:1784637/1‑148 (MQ=255)
         ccTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGaaaa                                                                                                                                          <  1:1112433/148‑1 (MQ=255)
          cTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGaaaaa                                                                                                                                         >  2:209477/1‑148 (MQ=255)
            cGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGTTGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAgg                                                                                                                                       <  1:1557324/148‑1 (MQ=255)
                 gAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGCCGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTa                                                                                                                                  <  2:155408/148‑1 (MQ=255)
                 gAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTa                                                                                                                                  <  2:3080793/148‑1 (MQ=255)
                 gAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTa                                                                                                                                  >  2:1340034/1‑148 (MQ=255)
                  aaGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTAt                                                                                                                                 <  1:2895462/148‑1 (MQ=255)
                  aaGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTAt                                                                                                                                 <  1:774713/148‑1 (MQ=255)
                  aaGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTAt                                                                                                                                 <  1:590089/148‑1 (MQ=255)
                    gCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCg                                                                                                                               <  1:2845605/148‑1 (MQ=255)
                    gCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCg                                                                                                                               <  2:1775611/148‑1 (MQ=255)
                    gCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCg                                                                                                                               >  1:1626833/1‑148 (MQ=255)
                     cTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGa                                                                                                                              >  1:2768332/1‑148 (MQ=255)
                     cTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGa                                                                                                                              >  1:1759137/1‑148 (MQ=255)
                     cTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGa                                                                                                                              >  1:3517299/1‑148 (MQ=255)
                      tATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAg                                                                                                                             <  1:1018315/148‑1 (MQ=255)
                          aCGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGc                                                                                                                         <  2:940014/148‑1 (MQ=255)
                           cGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCa                                                                                                                        <  1:959625/148‑1 (MQ=255)
                           cGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCa                                                                                                                        <  2:2259830/148‑1 (MQ=255)
                             gATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCata                                                                                                                      >  1:721011/1‑148 (MQ=255)
                             gATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCata                                                                                                                      <  1:2996141/148‑1 (MQ=255)
                              aTACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATAGAGCTGCatat                                                                                                                     >  1:1268293/1‑148 (MQ=255)
                                aCCGTATATCGTGTCCGCATCGTGCCCGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTg                                                                                                                   <  2:399335/148‑1 (MQ=255)
                                aCCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCATATTg                                                                                                                   <  1:383347/148‑1 (MQ=255)
                                aCCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTg                                                                                                                   >  2:616421/1‑148 (MQ=255)
                                aCCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTg                                                                                                                   <  2:398858/148‑1 (MQ=255)
                                    tataTCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATAGAGCTACATATTGgcgc                                                                                                               >  1:2466590/1‑148 (MQ=255)
                                     atatCGTGTCCGCATCGTGCCAGACGCGACCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCt                                                                                                              >  2:2034279/1‑148 (MQ=255)
                                     atatCGTGTCCGCATCGTGCCAGACGCGACCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCt                                                                                                              >  2:236976/1‑148 (MQ=255)
                                      tttCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCTTTCTCAATGGCTTCGACGCCGGATGAAAAAGGTTTTATCGAGCTGCATATTGGCGCtt                                                                                                             <  2:300266/146‑1 (MQ=255)
                                      tatCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCtt                                                                                                             <  2:3448091/148‑1 (MQ=255)
                                      tatCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCtt                                                                                                             <  1:2395579/148‑1 (MQ=255)
                                       atCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTc                                                                                                            <  1:838164/148‑1 (MQ=255)
                                        tCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCt                                                                                                           <  1:3525450/148‑1 (MQ=255)
                                         cGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCAGATTGGCGCTTCTg                                                                                                          >  2:1977800/1‑148 (MQ=255)
                                         cGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGGAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTg                                                                                                          >  1:3447299/1‑148 (MQ=255)
                                           tgtCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGTTGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGaa                                                                                                        >  1:2423107/1‑148 (MQ=255)
                                                 cATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAAc                                                                                                  <  1:2459688/148‑1 (MQ=255)
                                                 cATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAAc                                                                                                  <  1:260371/148‑1 (MQ=255)
                                                    cGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCtt                                                                                               <  1:212388/148‑1 (MQ=255)
                                                           gACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTTGCGCTTCTGAAATCAACCTTTAAGCGa                                                                                        >  2:2104717/1‑148 (MQ=255)
                                                               cgGCCTTTTCTTTTCGTGCTGGTCAGTATTTGAGGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGc                                                                                    <  2:2700528/148‑1 (MQ=255)
                                                                ggCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCa                                                                                   <  1:467261/148‑1 (MQ=255)
                                                                 gCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAg                                                                                  <  2:2289672/148‑1 (MQ=255)
                                                                   cttttcttttCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTc                                                                                <  2:2411598/148‑1 (MQ=255)
                                                                     tttcttttcGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCAt                                                                              <  2:331488/148‑1 (MQ=255)
                                                                     tttcttttcGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCAt                                                                              >  2:1729323/1‑148 (MQ=255)
                                                                        cttttcGTGCTTGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGa                                                                           <  2:1256892/148‑1 (MQ=255)
                                                                        cttttcGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGa                                                                           <  1:1824522/148‑1 (MQ=255)
                                                                         ttttcGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGAc                                                                          >  2:2787802/1‑148 (MQ=255)
                                                                               tGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATc                                                                    <  2:1784645/148‑1 (MQ=255)
                                                                                    gTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCaa                                                               >  2:3076038/1‑148 (MQ=255)
                                                                                    gTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCaa                                                               >  1:2873676/1‑148 (MQ=255)
                                                                                        gTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGat                                                           <  1:859237/148‑1 (MQ=255)
                                                                                         tATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAACAAGGGTTNATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGatc                                                          <  1:1806348/148‑1 (MQ=255)
                                                                                          aTTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGatca                                                         <  1:1086587/148‑1 (MQ=255)
                                                                                          aTTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGatca                                                         <  1:1023908/148‑1 (MQ=255)
                                                                                          aTTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGatca                                                         <  1:2559769/148‑1 (MQ=255)
                                                                                          aTTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGatca                                                         <  1:747254/148‑1 (MQ=255)
                                                                                                 ggTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATc                                                  >  1:3428729/1‑148 (MQ=255)
                                                                                                    aGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTg                                               >  1:1685991/1‑148 (MQ=255)
                                                                                                      tGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGt                                             >  1:2035293/1‑148 (MQ=255)
                                                                                                      tGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGt                                             >  1:2328486/1‑148 (MQ=255)
                                                                                                       gatggatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTc                                            >  1:2657762/1‑148 (MQ=255)
                                                                                                         tggatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGa                                          <  2:2757331/148‑1 (MQ=255)
                                                                                                         tggatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGa                                          <  1:2948625/148‑1 (MQ=255)
                                                                                                          ggatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGAc                                         <  1:1064094/148‑1 (MQ=255)
                                                                                                          ggatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGAc                                         >  2:86597/1‑148 (MQ=255)
                                                                                                           gatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACa                                        <  1:616421/148‑1 (MQ=255)
                                                                                                           gatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACa                                        <  2:378079/148‑1 (MQ=255)
                                                                                                           gatgAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACa                                        <  1:2062217/148‑1 (MQ=255)
                                                                                                               aGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTcc                                    <  2:1535501/148‑1 (MQ=255)
                                                                                                                 cgcgACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTcccc                                  >  1:3442178/1‑148 (MQ=255)
                                                                                                                 cgcgACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTcccc                                  <  1:1483042/148‑1 (MQ=255)
                                                                                                                     aCAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGAGGTCGACATTCCCCACgg                              >  2:130930/1‑148 (MQ=255)
                                                                                                                      cAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGga                             >  2:2535471/1‑148 (MQ=255)
                                                                                                                       aaaCGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGgag                            >  1:2024644/1‑148 (MQ=255)
                                                                                                                        aaCGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGgaga                           <  1:1434408/148‑1 (MQ=255)
                                                                                                                        aaCGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGgaga                           >  1:3241862/1‑148 (MQ=255)
                                                                                                                        aaCGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGgaga                           >  1:2324030/1‑148 (MQ=255)
                                                                                                                          cGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAg                         >  2:796704/1‑148 (MQ=255)
                                                                                                                            tCCGTTCTCAATGGCTTCGACGCCGGATGGAAAAGGGTTTATCGAGCTGCAGATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCg                       >  2:3444724/1‑148 (MQ=255)
                                                                                                                            tCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCg                       >  1:1378823/1‑148 (MQ=255)
                                                                                                                               gTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTgg                    <  1:1643589/148‑1 (MQ=255)
                                                                                                                                ttCTCAATGGCTTCGACGGCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGc                   >  1:319098/1‑148 (MQ=255)
                                                                                                                                ttCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGc                   <  1:2137031/148‑1 (MQ=255)
                                                                                                                                     aaTGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTgcgc              >  2:2460934/1‑148 (MQ=255)
                                                                                                                                      aTGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGcgcg             <  1:1134908/148‑1 (MQ=255)
                                                                                                                                        ggCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCgat           >  2:3176555/1‑148 (MQ=255)
                                                                                                                                        ggCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCgat           >  2:3062368/1‑148 (MQ=255)
                                                                                                                                         gCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCgatg          >  1:2095115/1‑148 (MQ=255)
                                                                                                                                          cTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCgatga         >  2:2415146/1‑148 (MQ=255)
                                                                                                                                          cTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCgatga         <  1:1047102/148‑1 (MQ=255)
                                                                                                                                            tCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGACGCGTGGCTGCGCGatgatg       >  2:1009481/1‑148 (MQ=255)
                                                                                                                                            tCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGatgatg       >  1:2147807/1‑148 (MQ=255)
                                                                                                                                                cGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAaga   <  2:2436396/148‑1 (MQ=255)
                                                                                                                                                 gCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAagag  <  2:2445928/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCAGACGCGGCCTTTTCTTTTCGTGCTGGTCAGTATTTGATGGTAGTGATGGATGAGCGCGACAAACGTCCGTTCTCAATGGCTTCGACGCCGGATGAAAAAGGGTTTATCGAGCTGCATATTGGCGCTTCTGAAATCAACCTTTACGCGAAAGCAGTCATGGACCGCATCCTCAAAGATCATCAAATCGTGGTCGACATTCCCCACGGAGAAGCGTGGCTGCGCGATGATGAAGAG  >  CP009273/4019903‑4020195

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