Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 373,857 G→A 6.4% intergenic (‑33/+61) frmB ← / ← frmA S‑formylglutathione hydrolase/alcohol dehydrogenase class III/glutathione‑dependent formaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273373,8570GA 6.4% 183.3 / 3.2 62intergenic (‑33/+61)frmB/frmAS‑formylglutathione hydrolase/alcohol dehydrogenase class III/glutathione‑dependent formaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base A (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

GGCGCAACAACAATAATGTTGTGCTCAGCCGCGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACAAACGGTTCCAGATCGATATCACCTTTCATCGCATCTTCAACCATGCCCGGTAACTGGGAACGACCTTTCACGCCGCCAAACGCGGAACCTTTCCATACGCGACCAGTGACC  >  CP009273/373589‑374119
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                      
ggCGCAACAACAATAATGTGGTGCTCAGCCGCGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGccc                                                                                                                                                                                                                                                                  >  1:649547/1‑275 (MQ=255)
     aacaacaaGAATGTTGTGCTCAGCCGCGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGTGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                                                                                                                                             <  1:149948/275‑1 (MQ=255)
     aaaaaaaaaaatGTTGTGCTCAGCAGAGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                                                                                                                                             <  2:649540/266‑1 (MQ=255)
                               cgTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCTGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCATTCATTATCTCAGGGTTAATCCCCCGAACGGTGTCCCGGGGAGAGTATAAAAAACCCCCGCCCGAAgg                                                                                                                                                                                                                                   >  1:781037/1‑274 (MQ=255)
                               cgTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGANATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAg                                                                                                                                                                                                                                    >  1:589733/1‑274 (MQ=255)
                                   aCGCTGCATCCCCGATTTAGTAATGTAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTTACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGa                                                                                                                                                                                                                                <  1:239975/274‑1 (MQ=255)
                                    cGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACg                                                                                                                                                                                                                              >  1:1435/1‑275 (MQ=255)
                                      cTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTGTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCACGGGGAGAGTCTTAAGAACCCCCGCCTGAAGGGGAcc                                                                                                                                                                                                                              >  2:239971/1‑272 (MQ=255)
                                                     aGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCg                                                                                                                                                                                                                                                                 >  2:774099/1‑223 (MQ=255)
                                                         aTGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTCTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                                                                                                                                             <  1:774108/223‑1 (MQ=255)
                                                                                   gCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCg                                                                                                                                                                                                                                                                 >  2:200539/1‑193 (MQ=255)
                                                                                   gCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATGTCAGTGATTGGGAATAATGCCGATAGATATTTTGCCAGCCGCCAAAGCTTACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCGGGGGAGAATATTAAGAACCACCGACTGACGGGGAGGGGGGAAACCTTGGGGAAAGAAGTACAGAATTACGGTTCGaagagat                                                                                                                                                                                >  2:192523/1‑269 (MQ=255)
                                                                                       gAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                                                                                                                                             <  1:200542/193‑1 (MQ=255)
                                                                                                     cacaGGCAATGTTTCATTCGAGGCTTTTGGGGGGAAATAGACGCCGACATTCATTTCACATTTCAGTGATGGGGAATAATGACGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTAACTCAGGGTTAATCGCCCGGACGGCTGCACGGGGAGAGTATTAAGAACAACCGCCTGAAGTGGAAGGGGTAAACCTGACGGAAATAAGGAACGCAATACGGTTCGAATGGATTTGCCTTCATGCATCaatg                                                                                                                                                             >  2:255053/1‑272 (MQ=255)
                                                                                                          gCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAg                                                                                                                                                                                                                                    >  1:453692/1‑199 (MQ=255)
                                                                                                          gCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAAACAGTAACGAATTACGGTTCGAATGGATTGGCCTGCCTGCATCAGGTcaacg                                                                                                                                                        >  2:253216/1‑271 (MQ=255)
                                                                                                              tttttCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGa                                                                                                                                                                                                                                <  2:453687/199‑1 (MQ=255)
                                                                                                                    tntagcgGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTTGGAATAATGCCGATACATATGTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCAt                                                                                                                                                 <  2:1435/268‑1 (MQ=255)
                                                                                                                     tttgcgcgttttGGTGGCGAATAAAANACGTCATTCAGTTCCCATTTCAGGTATTGGGAATAATTCCGTAACATATTTTGCCAGCCGCCAAAGCTGACATGTTGTTCAATGAGGTCCATCAGTCAGTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGGGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAAt                                                                                                                                             <  2:553379/267‑1 (MQ=255)
                                                                                                                             ttttGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCa                                                                                                                                     <  1:253220/275‑1 (MQ=255)
                                                                                                                                      gagaTAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATTAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATgg                                                                                                                            >  2:287527/1‑275 (MQ=255)
                                                                                                                                             aCGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCg                                                                                                                                                                                                                                                                 >  2:559190/1‑135 (MQ=255)
                                                                                                                                                 cGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                                                                                                                                             <  1:559195/135‑1 (MQ=255)
                                                                                                                                                    aattaatgTCACAGTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATATCAGGGTTAATCTCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACaaa                                                                                                               <  1:354934/273‑1 (MQ=255)
                                                                                                                                                     actaagtgaaagtttcaGGTAGTGGGAATAATGCCGATACATATTTTGCCAGGCGCACAAGCTGACATGTGTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACaaa                                                                                                               <  2:794503/261‑1 (MQ=255)
                                                                                                                                                                       tGATTGGGAATAATGCCGATACATCTTTGTCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCg                                                                                                                                                                                                                                                                 >  1:555347/1‑109 (MQ=255)
                                                                                                                                                                        gATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACAAACGGTTCCAGATCGATATCAcc                                                                                          <  1:192526/275‑1 (MQ=255)
                                                                                                                                                                           tGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                                                                                                                                             <  2:555342/109‑1 (MQ=255)
                                                                                                                                                                            gggAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGACTGccc                                                                                                                                                                                                                                                                  >  2:254743/1‑103 (MQ=255)
                                                                                                                                                                            gggAATAATGACGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATCCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGa                                                                                                                   >  1:794512/1‑246 (MQ=255)
                                                                                                                                                                                ataatGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGACTGCCCgggg                                                                                                                                                                                                                                                              <  1:254747/103‑1 (MQ=255)
                                                                                                                                                                                    tGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAg                                                                                                                                                                                                                                    >  1:145524/1‑125 (MQ=255)
                                                                                                                                                                                       cGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCg                                                                                                                                                                                                                                                                 >  1:268018/1‑93 (MQ=255)
                                                                                                                                                                                        gATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGa                                                                                                                                                                                                                                <  2:145524/125‑1 (MQ=255)
                                                                                                                                                                                         aTACATATTTTGGCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTACATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGGGTATTAAGAAccc                                                                                                                                                                                                                                              >  2:613791/1‑110 (MQ=255)
                                                                                                                                                                                          tACATATTTTGCCAGCCGCCAAAGCTGACATGTTATTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGACTGCCCg                                                                                                                                                                                                                                                                 >  1:552674/1‑90 (MQ=255)
                                                                                                                                                                                           aCATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                                                                                                                                             <  2:268014/93‑1 (MQ=255)
                                                                                                                                                                                             atatTTTGCCAGCCGCCTAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTGGATCCCCCGGCCGGCCGCCCGGGGAGAGTATTAAGAAGCCCCGc                                                                                                                                                                                                                                          <  1:613796/110‑1 (MQ=255)
                                                                                                                                                                                              tatTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGACTGCCCGGGga                                                                                                                                                                                                                                                             <  2:552669/90‑1 (MQ=255)
                                                                                                                                                                                                   tGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCg                                                                                                                                                                                                                                                                 >  2:39963/1‑81 (MQ=255)
                                                                                                                                                                                                       aGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGga                                                                                                                                                                                                                                                             <  1:39963/81‑1 (MQ=255)
                                                                                                                                                                                                         ccgccAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAg                                                                                                                                                                                                                                    >  2:192273/1‑104 (MQ=255)
                                                                                                                                                                                                             cAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGa                                                                                                                                                                                                                                <  1:192276/104‑1 (MQ=255)
                                                                                                                                                                                                                                      ttCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACAAACGGTTCCAGATCGATAGCACCTTTCATCGCATCTTCAACCATGCCCGGTAAATGGGAACGACATTTAACGCAGCCAAAcgag                             >  1:165933/1‑272 (MQ=255)
                                                                                                                                                                                                                                          ggaggTCATAAACGCACGGTGAAGCCCCCGGCCGGCGTCCCGGGGGGAGTAGTAATAACCCTCGACTGAAGTGGACGGGTTAAACCTGGGGTAAATCAGGAACGAATTACGGTTCGAATGGATTTGCCGTCATGCATCAGGTAGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACAAACGGTTCCAGATCGATATCACCTTTCATCGCAGCTTCAACCATGCCCGGTAACCGGGAACGACCTTTCACGCCGCCAAACGCGGaa                          <  2:232294/269‑1 (MQ=255)
                                                                                                                                                                                                                                             aGTCATTATCTTAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAANTCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACAAACGGTTCCAGATCGATATCACCTTTCATCGCATCTTCAACCATGCCGGGTAACTGGGAACGACCTTTCACGCCGCCAAACGCGGAACCttt                     <  1:591949/275‑1 (MQ=255)
                                                                                                                                                                                                                                                   tATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGa                                                                                                                                                                                 >  1:561543/1‑113 (MQ=255)
                                                                                                                                                                                                                                                     tctcAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGt                                                                                                                     >  1:480333/1‑171 (MQ=255)
                                                                                                                                                                                                                                                     tctcAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACAAACGGTTCCAGATCGATATCACCTTTCATCGCATCTTCAACCATGCCCGGTAACTGGGAACGACCTTTCACGCCGCCAAACGCGGAACCTTTCCATAcg              <  1:287531/274‑1 (MQ=255)
                                                                                                                                                                                                                                                       tcAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTg                                                                                                                                                                             <  2:561538/113‑1 (MQ=255)
                                                                                                                                                                                                                                                         aGGGTTAATGCCCAGGACGGCTGACCGGGGAGAGTAGTAAGAACCCCCGCCTGAAGTGGACGGGGTAATCCTTCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACa                                                                                                                 <  2:480328/171‑1 (MQ=255)
                                                                                                                                                                                                                                                         aGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAg                                                                                                                                                                                                                                    >  2:673515/1‑56 (MQ=255)
                                                                                                                                                                                                                                                         aGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAAttt                                                                                                                                           >  1:134416/1‑145 (MQ=255)
                                                                                                                                                                                                                                                         aGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATc                                                                                                                                       >  1:500244/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         aGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACAAACGGTTCCAGATCGATATCACCTTTCATCGCATCTTCAACCATGCCCAGTAACCGGGAACGACCTTTCCCGCCGCCAAACTCGAACCCTTTcacacacctac          >  1:232298/1‑264 (MQ=255)
                                                                                                                                                                                                                                                             ttAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGa                                                                                                                                                                                                                                <  1:673522/56‑1 (MQ=255)
                                                                                                                                                                                                                                                             ttAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATc                                                                                                                                       <  2:134416/145‑1 (MQ=255)
                                                                                                                                                                                                                                                             ttAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAgg                                                                                                                                   <  2:500239/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                aTCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCa                                                                                                                                     >  1:364544/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                aTCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACAAACGGTTCCAGATCGATATCACCTTTCATCGCATCTTCAACCATGACCGTTAACTGGGAACGACCTTTAACGCCGCCAAACGCGNGACCGTTCCAGACGccaccaccgctc  >  2:550734/1‑269 (MQ=255)
                                                                                                                                                                                                                                                                    cccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGNGGGAAATCNGTAACGAATTACTGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACAAACGGTTCCAGATCGATATCACCTTTCATCGCATCTTCAACCATGCCCGGTAACTGGGAATGACCTTTGACGCCGCCAAACGCGGAACCTTTCCa                  <  1:550739/255‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cccGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCAAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCt                                                                                                                                 <  2:364539/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                      
GGCGCAACAACAATAATGTTGTGCTCAGCCGCGTAACGCTGCATCCCCGATTTAGTAATGAAATTCTGCTCGTTGCAGGTCAGGCCTGAAAGCCAGTACAGCACCGGCAATTTTTCATTCGCGGCTTTTGGTGGGAGATAGACGCCGACATTCATTTCACATTTCAGTGATTGGGAATAATGCCGATACATATTTTGCCAGCCGCCAAAGCTGACATGTTTTTCAATGAGTTCCATCAGTCATTATCTCAGGGTTAATCCCCCGGACGGCTGCCCGGGGAGAGTATTAAGAACCCCCGCCTGAAGTGGACGGGGTAAACCTGCGGGAAATCAGTAACGAATTACGGTTCGAATGGATTTGCCTTCATGCATCAGGTCGAAGGCGTCATTAATTTCATCCAGGCTCATGGTATGCGTGACAAACGGTTCCAGATCGATATCACCTTTCATCGCATCTTCAACCATGCCCGGTAACTGGGAACGACCTTTCACGCCGCCAAACGCGGAACCTTTCCATACGCGACCAGTGACC  >  CP009273/373589‑374119

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑