Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,481,974 G→T 100% 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,9740GT99.5% 601.7 / ‑4.3 183S225* (TCA→TAA) yfdEacetyl‑CoA:oxalate CoA‑transferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (100/82);  total (100/83)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 191 total aligned reads displayed.

TTAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAGACAC  >  CP009273/2481830‑2482118
                                                                                                                                                |                                                                                                                                                
ttAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACt                                                                                                                                               >  1:536852/1‑148 (MQ=255)
ttAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACt                                                                                                                                               >  2:2468993/1‑148 (MQ=255)
 tAACCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACtt                                                                                                                                              >  1:3395661/1‑148 (MQ=255)
  aaCCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTc                                                                                                                                             >  1:3119657/1‑148 (MQ=255)
   aCCAGTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTcc                                                                                                                                            >  2:33313/1‑148 (MQ=255)
       gTTCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGt                                                                                                                                        <  2:2154573/148‑1 (MQ=255)
        ttCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGtt                                                                                                                                       >  1:2378962/1‑148 (MQ=255)
        ttCCGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGtt                                                                                                                                       >  1:614378/1‑148 (MQ=255)
           cGTAAGCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCg                                                                                                                                    <  1:993816/148‑1 (MQ=255)
                gCTCCAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGatata                                                                                                                               <  1:2619124/148‑1 (MQ=255)
                   ccAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGc                                                                                                                            >  1:964874/1‑148 (MQ=255)
                    cAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGcc                                                                                                                           <  1:2454145/148‑1 (MQ=255)
                    cAGTGCCTGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGcc                                                                                                                           <  1:1811152/148‑1 (MQ=255)
                     aGTGCCGGGCATAACGCATAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCa                                                                                                                          <  2:1850149/148‑1 (MQ=255)
                       tGCCTGGCATAACGCAGAAAAAAGCTTTTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATc                                                                                                                        >  1:1429465/1‑148 (MQ=255)
                           tGGCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGAc                                                                                                                    >  1:523678/1‑148 (MQ=255)
                            ggCATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGAcc                                                                                                                   >  1:1028681/1‑148 (MQ=255)
                              cATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCAt                                                                                                                 >  2:1684254/1‑148 (MQ=255)
                              cATAACGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCAt                                                                                                                 >  2:3491970/1‑148 (MQ=255)
                                  aCGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTc                                                                                                             >  1:2715365/1‑148 (MQ=255)
                                   cGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTcc                                                                                                            <  2:292874/148‑1 (MQ=255)
                                   cGCAGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTcc                                                                                                            >  2:1742285/1‑148 (MQ=255)
                                      aGAAAAAAGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGa                                                                                                         >  2:786265/1‑148 (MQ=255)
                                        aaaaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGaaa                                                                                                       >  1:233917/1‑148 (MQ=255)
                                        aaaaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGaaa                                                                                                       >  1:36232/1‑148 (MQ=255)
                                         aaaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGaaaa                                                                                                      <  2:201961/148‑1 (MQ=255)
                                          aaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTTAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGAGATATATGCCATCAGACCATGATCCAGAAAAc                                                                                                     >  1:600666/1‑148 (MQ=255)
                                          aaaaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAAc                                                                                                     >  2:1177820/1‑148 (MQ=255)
                                            aaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTc                                                                                                   <  1:105328/148‑1 (MQ=255)
                                            aaGCTTGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTc                                                                                                   >  1:1543398/1‑148 (MQ=255)
                                              gattGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAg                                                                                                 <  2:2099395/146‑1 (MQ=255)
                                                 tGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGt                                                                                              >  1:2326665/1‑148 (MQ=255)
                                                 tGTCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGt                                                                                              <  1:1742285/148‑1 (MQ=255)
                                                   tCATTACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTgg                                                                                            >  2:3090896/1‑148 (MQ=255)
                                                      ttACCACAACACATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCAt                                                                                         <  2:2073004/148‑1 (MQ=255)
                                                      ttACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACCTCTAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCAt                                                                                         <  2:1156861/148‑1 (MQ=255)
                                                      ttACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCAt                                                                                         <  2:1249485/148‑1 (MQ=255)
                                                      ttACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCAt                                                                                         <  2:1663236/148‑1 (MQ=255)
                                                       tACCACAACAAATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATc                                                                                        >  1:3392399/1‑148 (MQ=255)
                                                             aacaaATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACAt                                                                                  <  1:2997860/148‑1 (MQ=255)
                                                               caaATCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCg                                                                                >  2:87478/1‑148 (MQ=255)
                                                                aaaTCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGc                                                                               <  1:1177820/148‑1 (MQ=255)
                                                                 aaTCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCt                                                                              >  2:2624727/1‑148 (MQ=255)
                                                                 aaTCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCt                                                                              >  1:772209/1‑148 (MQ=255)
                                                                  aTCGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCta                                                                             >  1:2022006/1‑148 (MQ=255)
                                                                    cGTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCtata                                                                           >  2:1125489/1‑148 (MQ=255)
                                                                     gTAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCCGCGTGGCATCAAACATCGCTatat                                                                          <  1:33313/148‑1 (MQ=255)
                                                                      tAATCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATc                                                                         >  2:505124/1‑148 (MQ=255)
                                                                       aaTCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCg                                                                        >  1:1703412/1‑148 (MQ=255)
                                                                        aTCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGa                                                                       >  1:2168647/1‑148 (MQ=255)
                                                                        aTCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGa                                                                       >  1:22706/1‑148 (MQ=255)
                                                                         tCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGAc                                                                      >  2:564378/1‑148 (MQ=255)
                                                                         tCGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGAc                                                                      <  1:2057579/148‑1 (MQ=255)
                                                                          cGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACa                                                                     <  1:842521/148‑1 (MQ=255)
                                                                          cGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACa                                                                     >  1:1996648/1‑148 (MQ=255)
                                                                          cGGCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACa                                                                     >  1:1996861/1‑148 (MQ=255)
                                                                            gCTTATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATg                                                                   <  1:2889444/148‑1 (MQ=255)
                                                                              ttATCCTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTGACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATgcg                                                                 <  2:2738959/148‑1 (MQ=255)
                                                                                  ccTGAGTATTGAAAACATCAAAAGGAGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGcccc                                                             >  2:2955843/1‑148 (MQ=255)
                                                                                   cTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCct                                                            <  2:935634/148‑1 (MQ=255)
                                                                                   cTGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCct                                                            <  1:752655/148‑1 (MQ=255)
                                                                                    tGAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCctc                                                           >  2:1911565/1‑148 (MQ=255)
                                                                                     gAGTATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCctct                                                          <  2:320845/148‑1 (MQ=255)
                                                                                        tATTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTg                                                       >  2:1499562/1‑148 (MQ=255)
                                                                                         aTTGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTgg                                                      >  1:2093309/1‑148 (MQ=255)
                                                                                          ttGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGc                                                     >  2:1823379/1‑148 (MQ=255)
                                                                                          ttGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGc                                                     <  1:2593506/148‑1 (MQ=255)
                                                                                           tGAAAACATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGct                                                    >  2:819438/1‑148 (MQ=255)
                                                                                              aaaCATCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCtt                                                 >  2:2469564/1‑148 (MQ=255)
                                                                                                  aTCAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTcg                                             >  1:231622/1‑148 (MQ=255)
                                                                                                    cAAAAGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTcgcg                                           >  1:2220709/1‑148 (MQ=255)
                                                                                                       aaGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGcc                                        <  1:3316716/148‑1 (MQ=255)
                                                                                                       aaGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGcc                                        <  1:139467/148‑1 (MQ=255)
                                                                                                        aGGTGCCATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCa                                       >  2:2960720/1‑148 (MQ=255)
                                                                                                             ccATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAg                                  <  1:1082587/148‑1 (MQ=255)
                                                                                                              cATGTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGt                                 >  1:2511042/1‑148 (MQ=255)
                                                                                                                 gTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGca                              <  1:536973/148‑1 (MQ=255)
                                                                                                                 gTAGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGca                              >  1:2128216/1‑148 (MQ=255)
                                                                                                                   aGGGATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGcact                            <  2:2597015/148‑1 (MQ=255)
                                                                                                                    gggATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCATAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGcactc                           >  1:328414/1‑148 (MQ=255)
                                                                                                                     ggATGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGcactca                          <  1:2929852/148‑1 (MQ=255)
                                                                                                                       aTGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGcactcact                        <  2:1332117/148‑1 (MQ=255)
                                                                                                                        tGGCGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCNTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTa                       <  1:277661/148‑1 (MQ=255)
                                                                                                                           cGATTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTc                    <  1:1277698/148‑1 (MQ=255)
                                                                                                                             aTTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCg                  <  2:466051/148‑1 (MQ=255)
                                                                                                                             aTTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCg                  <  2:628728/148‑1 (MQ=255)
                                                                                                                             aTTTCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCg                  >  2:1767980/1‑148 (MQ=255)
                                                                                                                              tttCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGc                 >  1:1147661/1‑148 (MQ=255)
                                                                                                                              tttCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGc                 >  2:1871160/1‑148 (MQ=255)
                                                                                                                               ttCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCt                >  1:924161/1‑148 (MQ=255)
                                                                                                                               ttCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCt                <  2:3494890/148‑1 (MQ=255)
                                                                                                                                tCCCAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTg               <  1:663086/148‑1 (MQ=255)
                                                                                                                                  ccAGACGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGaa             <  2:756820/148‑1 (MQ=255)
                                                                                                                                      aCGTTGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGaataaa         <  1:1810604/148‑1 (MQ=255)
                                                                                                                                         ttGTGGTTACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAg      <  2:273112/148‑1 (MQ=255)
                                                                                                                                           gtggttACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAGac    <  2:3402572/148‑1 (MQ=255)
                                                                                                                                           gtggttACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAGac    <  1:845101/148‑1 (MQ=255)
                                                                                                                                            tggttACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAGaca   <  2:360229/148‑1 (MQ=255)
                                                                                                                                             ggttACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAGacac  <  2:3468591/148‑1 (MQ=255)
                                                                                                                                             ggttACTTCCCAGTTGCGATATATGCCATCAGACCATGCTCCAGAAAACTCAGCGTGGCATCAAACATCGCTATATCGACATGCGCCCCTCTCTGGCTCTTTTCGCGGCCATAAAGTGCACTCACTATTCCGCTGAATAAATAGacac  >  1:3051242/1‑148 (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: