Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,055,132 G→T 100% E141* (GAA→TAA)  torA → trimethylamine N‑oxide (TMAO) reductase I, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,055,1320GT99.5% 662.7 / ‑4.1 206E141* (GAA→TAA) torAtrimethylamine N‑oxide (TMAO) reductase I, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (101/104);  total (101/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-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 213 total aligned reads displayed.

TCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAG  >  CP009273/1054987‑1055276
                                                                                                                                                 |                                                                                                                                                
tCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATgaataa                                                                                                                                                <  1:642666/148‑1 (MQ=255)
tCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATgaataa                                                                                                                                                >  1:1661264/1‑148 (MQ=255)
  aCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACt                                                                                                                                              >  1:1072469/1‑148 (MQ=255)
  aCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACt                                                                                                                                              <  2:472685/148‑1 (MQ=255)
  aCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACt                                                                                                                                              >  1:1861286/1‑148 (MQ=255)
  aCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACt                                                                                                                                              <  2:1607819/148‑1 (MQ=255)
   cGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTg                                                                                                                                             >  1:2231773/1‑148 (MQ=255)
    gTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTgg                                                                                                                                            >  1:2924126/1‑148 (MQ=255)
    gTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTgg                                                                                                                                            <  1:832532/148‑1 (MQ=255)
     tACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGa                                                                                                                                           >  1:3370494/1‑148 (MQ=255)
          aaCGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAAcgcg                                                                                                                                      >  1:3234913/1‑148 (MQ=255)
          aaCGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAAcgcg                                                                                                                                      >  1:2531524/1‑148 (MQ=255)
                gcgcGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGa                                                                                                                                >  2:479952/1‑148 (MQ=255)
                  gcgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGaaa                                                                                                                              <  2:1395361/148‑1 (MQ=255)
                     tATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACt                                                                                                                           <  1:497623/148‑1 (MQ=255)
                      aTTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTc                                                                                                                          >  2:2097302/1‑148 (MQ=255)
                      aTTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTc                                                                                                                          >  1:790183/1‑148 (MQ=255)
                       ttCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCa                                                                                                                         >  2:2668402/1‑148 (MQ=255)
                        tCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCAc                                                                                                                        >  2:2644468/1‑148 (MQ=255)
                               ccGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGa                                                                                                                 >  1:1471280/1‑148 (MQ=255)
                                cGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAg                                                                                                                >  2:2793106/1‑148 (MQ=255)
                                 gATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGt                                                                                                               >  2:170247/1‑148 (MQ=255)
                                  aTGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTg                                                                                                              >  2:2557058/1‑148 (MQ=255)
                                    ggTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGcc                                                                                                            <  2:1376172/148‑1 (MQ=255)
                                    ggTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGcc                                                                                                            <  2:2118543/148‑1 (MQ=255)
                                        cgcgTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGc                                                                                                        >  1:2800755/1‑148 (MQ=255)
                                          cgTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTg                                                                                                      >  1:2992019/1‑148 (MQ=255)
                                           gTGGACTGGCTGCGTTAGCGCCATCTCAGCGATACCTCCCAGCGGGGTGATAACGGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTGGACTTGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGa                                                                                                     <  1:35863/148‑1 (MQ=255)
                                           gTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGa                                                                                                     <  2:1007476/148‑1 (MQ=255)
                                               aCTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGCCATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGc                                                                                                 <  2:1823814/148‑1 (MQ=255)
                                                         tAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTgg                                                                                       >  1:660851/1‑148 (MQ=255)
                                                           aGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCa                                                                                     <  1:755208/148‑1 (MQ=255)
                                                           aGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCa                                                                                     <  1:3495746/148‑1 (MQ=255)
                                                            gcgcCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCac                                                                                    >  1:2079259/1‑147 (MQ=255)
                                                               ccATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCg                                                                                 >  2:427418/1‑148 (MQ=255)
                                                                cATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGa                                                                                >  2:2733308/1‑148 (MQ=255)
                                                                   ctcAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgg                                                                             >  2:457569/1‑148 (MQ=255)
                                                                     cAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgggg                                                                           >  2:1413556/1‑148 (MQ=255)
                                                                     cAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgggg                                                                           >  2:1952313/1‑148 (MQ=255)
                                                                     cAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgggg                                                                           >  1:2891362/1‑148 (MQ=255)
                                                                     cAGCGATACCTCCCAGCGCGGTGATAACCGGTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgggg                                                                           >  1:1139833/1‑148 (MQ=255)
                                                                      aGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGa                                                                          <  2:3024447/148‑1 (MQ=255)
                                                                        cGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATg                                                                        <  2:1104044/148‑1 (MQ=255)
                                                                        cGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATg                                                                        <  1:183034/148‑1 (MQ=255)
                                                                           tACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTc                                                                     <  2:3309621/148‑1 (MQ=255)
                                                                           tACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTc                                                                     <  1:1237557/148‑1 (MQ=255)
                                                                           tACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTc                                                                     <  2:1352447/148‑1 (MQ=255)
                                                                             ccTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCa                                                                   >  1:2410091/1‑148 (MQ=255)
                                                                               tCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATa                                                                 <  2:337493/148‑1 (MQ=255)
                                                                                cccAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATaa                                                                >  1:2577049/1‑148 (MQ=255)
                                                                                   aGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGc                                                             <  2:850535/148‑1 (MQ=255)
                                                                                     cgcgGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCtt                                                           <  1:2473087/148‑1 (MQ=255)
                                                                                       cgGTGATACCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCg                                                         <  2:1015702/148‑1 (MQ=255)
                                                                                        ggTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCgg                                                        >  2:1211774/1‑148 (MQ=255)
                                                                                          tGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCgggg                                                      <  2:3213046/148‑1 (MQ=255)
                                                                                           gATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGa                                                     >  2:2966644/1‑148 (MQ=255)
                                                                                           gATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGa                                                     >  1:434424/1‑148 (MQ=255)
                                                                                            aTAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGAt                                                    >  2:616072/1‑148 (MQ=255)
                                                                                            aTAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGAt                                                    >  1:1090585/1‑148 (MQ=255)
                                                                                            aTAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGAt                                                    >  1:2548234/1‑148 (MQ=255)
                                                                                             tAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCTGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATg                                                   <  2:1261316/148‑1 (MQ=255)
                                                                                              aaCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATACCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGc                                                  <  2:789689/148‑1 (MQ=255)
                                                                                              aaCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGc                                                  <  1:884673/148‑1 (MQ=255)
                                                                                               aCCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCt                                                 >  1:692559/1‑148 (MQ=255)
                                                                                                 cGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTgg                                               <  2:3397581/148‑1 (MQ=255)
                                                                                                       gtgCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAg                                         <  1:3018430/148‑1 (MQ=255)
                                                                                                       gtgCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAg                                         <  1:3077513/148‑1 (MQ=255)
                                                                                                         gcgcGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCt                                       <  2:2217674/148‑1 (MQ=255)
                                                                                                           gcgTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTAt                                     <  1:526496/148‑1 (MQ=255)
                                                                                                           gcgTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTAt                                     <  1:3320924/148‑1 (MQ=255)
                                                                                                              tGAGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGc                                  <  2:2867770/148‑1 (MQ=255)
                                                                                                                aGCTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCt                                <  2:1567836/148‑1 (MQ=255)
                                                                                                                  cTGGGATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTg                              >  2:954847/1‑148 (MQ=255)
                                                                                                                   tGGGAGGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGGGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGc                             <  1:3231931/148‑1 (MQ=255)
                                                                                                                    gggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCa                            <  2:2522745/148‑1 (MQ=255)
                                                                                                                     ggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCAt                           <  2:1510079/148‑1 (MQ=255)
                                                                                                                     ggATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCAt                           >  1:1958211/1‑148 (MQ=255)
                                                                                                                      gATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATg                          <  2:2382248/148‑1 (MQ=255)
                                                                                                                      gATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATg                          >  1:326203/1‑148 (MQ=255)
                                                                                                                      gATGAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATg                          >  2:3215119/1‑148 (MQ=255)
                                                                                                                      gATGAAGCCCACGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATg                          <  2:3112130/148‑1 (MQ=255)
                                                                                                                         gAAGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTa                       >  2:2360153/1‑148 (MQ=255)
                                                                                                                          aaGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTaa                      <  1:1594232/148‑1 (MQ=255)
                                                                                                                          aaGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTaa                      <  1:634728/148‑1 (MQ=255)
                                                                                                                          aaGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTaa                      <  1:2397745/148‑1 (MQ=255)
                                                                                                                          aaGCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTaa                      <  2:308428/148‑1 (MQ=255)
                                                                                                                            gCCCTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATa                    >  1:1387735/1‑148 (MQ=255)
                                                                                                                             cccTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAg                   >  2:2479675/1‑148 (MQ=255)
                                                                                                                              ccTCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGc                  >  1:2131392/1‑148 (MQ=255)
                                                                                                                                tCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGt                <  1:566880/148‑1 (MQ=255)
                                                                                                                                tCGACATGTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGt                <  1:861413/148‑1 (MQ=255)
                                                                                                                                       gTTCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACg         >  2:536254/1‑148 (MQ=255)
                                                                                                                                        ttCTATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACgg        >  1:634889/1‑148 (MQ=255)
                                                                                                                                           tATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCg     >  1:2328541/1‑148 (MQ=255)
                                                                                                                                           tATGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCg     <  2:688701/148‑1 (MQ=255)
                                                                                                                                            aTGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCgg    >  1:3016015/1‑148 (MQ=255)
                                                                                                                                            aTGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCgg    >  1:1268369/1‑148 (MQ=255)
                                                                                                                                             tGAATAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGga   >  1:2684009/1‑148 (MQ=255)
                                                                                                                                              gaataaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGgag  >  1:2729382/1‑148 (MQ=255)
                                                                                                                                              gaataaCTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGgag  <  2:2924126/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTGCCTTGCATGGTAATAGCGTTGGTACGGGCGGAG  >  CP009273/1054987‑1055276

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: