Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,481,974 G→T 88.0% S225* (TCA→TAA)  yfdE ← acetyl‑CoA:oxalate CoA‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,481,9740GT88.0% 568.1 / 79.8 233S225* (TCA→TAA) yfdEacetyl‑CoA:oxalate CoA‑transferase
Reads supporting (aligned to +/- strand):  ref base G (12/16);  new base T (110/95);  total (122/111)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 237 total aligned reads displayed.

TTAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAGACAC  >  CP009273/2481830‑2482118
                                                                                                                                                |                                                                                                                                                
ttAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACt                                                                                                                                               >  1:1448407/1‑148 (MQ=255)
 tAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACtt                                                                                                                                              >  2:1124599/1‑148 (MQ=255)
   aCCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTcc                                                                                                                                            <  2:679130/148‑1 (MQ=255)
    ccAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTccc                                                                                                                                           <  2:921933/148‑1 (MQ=255)
    ccAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTccc                                                                                                                                           <  1:1569413/148‑1 (MQ=255)
      aGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAg                                                                                                                                         >  1:1814073/1‑148 (MQ=255)
         tCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTg                                                                                                                                      <  2:2951710/148‑1 (MQ=255)
             tAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGat                                                                                                                                  <  2:1394849/148‑1 (MQ=255)
              aaGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGata                                                                                                                                 >  2:96561/1‑148 (MQ=255)
              aaGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGata                                                                                                                                 >  2:2736733/1‑148 (MQ=255)
              aaGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGata                                                                                                                                 >  1:2656233/1‑148 (MQ=255)
              aaGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGata                                                                                                                                 >  2:773171/1‑148 (MQ=255)
              aaGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGata                                                                                                                                 >  2:773225/1‑148 (MQ=255)
              aaGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGata                                                                                                                                 >  1:643375/1‑148 (MQ=255)
                gCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGatata                                                                                                                               >  1:805008/1‑148 (MQ=255)
                gCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGatata                                                                                                                               >  1:601255/1‑148 (MQ=255)
                  tCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATg                                                                                                                             >  2:1987560/1‑148 (MQ=255)
                  tCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATg                                                                                                                             >  2:2334703/1‑148 (MQ=255)
                   ccAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGc                                                                                                                            >  1:2297668/1‑148 (MQ=255)
                    cAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGcc                                                                                                                           <  1:3036453/148‑1 (MQ=255)
                      gTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCAt                                                                                                                         <  2:1281258/148‑1 (MQ=255)
                       tGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATc                                                                                                                        <  1:812791/148‑1 (MQ=255)
                          cTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGa                                                                                                                     <  1:1226297/148‑1 (MQ=255)
                           tGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGAc                                                                                                                    >  1:1881256/1‑148 (MQ=255)
                              cATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCAt                                                                                                                 >  1:293301/1‑148 (MQ=255)
                               aTAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATg                                                                                                                <  2:2442214/148‑1 (MQ=255)
                                  aCGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTc                                                                                                             >  1:1996612/1‑148 (MQ=255)
                                  aCGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTc                                                                                                             >  2:2562857/1‑148 (MQ=255)
                                  aCGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTc                                                                                                             <  2:907544/148‑1 (MQ=255)
                                   cGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTcc                                                                                                            <  1:2553307/148‑1 (MQ=255)
                                   cGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTcc                                                                                                            >  1:1144829/1‑148 (MQ=255)
                                     cAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAg                                                                                                          >  2:679931/1‑148 (MQ=255)
                                      aGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGa                                                                                                         >  2:2625567/1‑148 (MQ=255)
                                        aaaaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGaaa                                                                                                       >  2:330292/1‑148 (MQ=255)
                                        aaaaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGaaa                                                                                                       >  1:1045754/1‑148 (MQ=255)
                                        aaaaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGaaa                                                                                                       >  1:3020055/1‑148 (MQ=255)
                                         aaaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGaaaa                                                                                                      <  1:2583372/148‑1 (MQ=255)
                                         aaaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGaaaa                                                                                                      <  2:578315/148‑1 (MQ=255)
                                          aaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAAc                                                                                                     >  1:1969174/1‑148 (MQ=255)
                                          aaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAAc                                                                                                     >  1:709757/1‑148 (MQ=255)
                                            aaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTc                                                                                                   <  1:88321/148‑1 (MQ=255)
                                             aGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCa                                                                                                  <  2:2946471/148‑1 (MQ=255)
                                              gCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAg                                                                                                 >  2:165936/1‑148 (MQ=255)
                                               cTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGc                                                                                                >  1:2093099/1‑148 (MQ=255)
                                               cTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGc                                                                                                >  2:443640/1‑148 (MQ=255)
                                                   tCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTgg                                                                                            >  2:2670855/1‑148 (MQ=255)
                                                    cATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGc                                                                                           >  1:2673024/1‑148 (MQ=255)
                                                     aTTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCa                                                                                          <  2:2859016/148‑1 (MQ=255)
                                                      ttACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCAt                                                                                         >  2:2832812/1‑148 (MQ=255)
                                                        aCCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCa                                                                                       >  1:3243678/1‑148 (MQ=255)
                                                         ccACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCaa                                                                                      >  1:2499831/1‑148 (MQ=255)
                                                         ccACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCaa                                                                                      <  2:1592329/148‑1 (MQ=255)
                                                           accacaaATCGGAATCGGCTTAGCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAAc                                                                                    <  2:2261783/145‑1 (MQ=255)
                                                           acaacaAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAAc                                                                                    <  1:497630/148‑1 (MQ=255)
                                                            caacaaATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACa                                                                                   <  1:367401/148‑1 (MQ=255)
                                                            caacaaATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACa                                                                                   >  1:138199/1‑148 (MQ=255)
                                                                  aTCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCta                                                                             <  2:502631/148‑1 (MQ=255)
                                                                    cGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCtata                                                                           >  2:3230288/1‑148 (MQ=255)
                                                                     gTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTatat                                                                          <  1:2219609/148‑1 (MQ=255)
                                                                      tAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATc                                                                         >  2:1029040/1‑148 (MQ=255)
                                                                          cGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACa                                                                     >  1:2364318/1‑148 (MQ=255)
                                                                           ggCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACAt                                                                    >  2:1927547/1‑148 (MQ=255)
                                                                             cTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATgc                                                                  <  2:2604224/148‑1 (MQ=255)
                                                                                 tCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGccc                                                              <  1:1728356/148‑1 (MQ=255)
                                                                                  ccTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGcccc                                                             >  2:618582/1‑148 (MQ=255)
                                                                                  ccTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGcccc                                                             >  1:1032611/1‑148 (MQ=255)
                                                                                  ccTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGcccc                                                             <  1:40599/148‑1 (MQ=255)
                                                                                  ccTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGcccc                                                             <  2:389478/148‑1 (MQ=255)
                                                                                   cTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCct                                                            <  2:727221/148‑1 (MQ=255)
                                                                                    tGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCctc                                                           <  1:2395733/148‑1 (MQ=255)
                                                                                    tGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCctc                                                           <  2:1144829/148‑1 (MQ=255)
                                                                                      aGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCctctc                                                         >  1:3156492/1‑148 (MQ=255)
                                                                                      aGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCctctc                                                         >  1:12734/1‑148 (MQ=255)
                                                                                        tATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTg                                                       >  1:1323142/1‑148 (MQ=255)
                                                                                           tGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGct                                                    <  2:1448407/148‑1 (MQ=255)
                                                                                            gAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGctc                                                   <  1:1422885/148‑1 (MQ=255)
                                                                                               aaCATCAACAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCttt                                                <  2:2998021/148‑1 (MQ=255)
                                                                                               aaCATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCttt                                                >  1:1084518/1‑148 (MQ=255)
                                                                                                aCATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCtttt                                               >  1:625356/1‑148 (MQ=255)
                                                                                                   tCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTcgc                                            >  2:1887022/1‑148 (MQ=255)
                                                                                                   tCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTcgc                                            >  2:1904695/1‑148 (MQ=255)
                                                                                                       aaGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGcc                                        >  2:2207333/1‑148 (MQ=255)
                                                                                                            gCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATaaa                                   >  2:1554286/1‑148 (MQ=255)
                                                                                                                tGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTGCGCGGCCATACAGTGc                               >  2:1021611/1‑148 (MQ=255)
                                                                                                                 gTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGca                              <  1:2704438/148‑1 (MQ=255)
                                                                                                                   aGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGcact                            <  1:1072563/148‑1 (MQ=255)
                                                                                                                    gggATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGcactc                           <  2:276405/148‑1 (MQ=255)
                                                                                                                       aTGGTGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGcactcact                        <  2:2066666/148‑1 (MQ=255)
                                                                                                                        tGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTa                       <  1:2889359/148‑1 (MQ=255)
                                                                                                                        tGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTa                       <  1:2888956/148‑1 (MQ=255)
                                                                                                                         ggCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTAt                      <  2:183239/148‑1 (MQ=255)
                                                                                                                         ggCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTAt                      <  1:1602919/148‑1 (MQ=255)
                                                                                                                           cGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTc                    <  2:3154296/148‑1 (MQ=255)
                                                                                                                                 cccAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGa              <  2:3114437/148‑1 (MQ=255)
                                                                                                                                 cccAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGa              <  1:724581/148‑1 (MQ=255)
                                                                                                                                  ccAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGaa             <  1:750769/148‑1 (MQ=255)
                                                                                                                                        gTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGaataaata       >  1:1842089/1‑148 (MQ=255)
                                                                                                                                         ttGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAg      <  1:445492/148‑1 (MQ=255)
                                                                                                                                         ttGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAg      >  2:205358/1‑148 (MQ=255)
                                                                                                                                             ggtgACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAGacac  <  2:2100406/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TTAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAGACAC  >  CP009273/2481830‑2482118

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