Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,312,055 C→A 100% G20G (GGC→GGA)  argG → argininosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,312,0550CA95.2% 613.7 / 19.2 206G20G (GGC→GGA) argGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base A (104/92);  total (108/98)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 206 total aligned reads displayed.

ATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGGAGAACGC  >  CP009273/3311911‑3312198
                                                                                                                                                |                                                                                                                                               
aTGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGt                                                                                                                                              >  1:1563411/1‑148 (MQ=255)
 tGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTc                                                                                                                                             >  2:2794318/1‑148 (MQ=255)
  gAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCt                                                                                                                                            <  2:821590/148‑1 (MQ=255)
   aaaaCTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTg                                                                                                                                           >  2:3301817/1‑148 (MQ=255)
    aaaCTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTgg                                                                                                                                          <  1:1796010/148‑1 (MQ=255)
       cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGaca                                                                                                                                       >  1:2318724/1‑148 (MQ=255)
        tCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGacac                                                                                                                                      >  2:41897/1‑148 (MQ=255)
         cATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACAcc                                                                                                                                     >  2:3029860/1‑148 (MQ=255)
          atttatttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCa                                                                                                                                    >  2:2270859/1‑148 (MQ=255)
            ttatttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTCGGTCAACGGATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGt                                                                                                                                  >  2:2583909/1‑148 (MQ=255)
             tatttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTg                                                                                                                                 <  1:2930742/148‑1 (MQ=255)
               ttttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGcc                                                                                                                               >  2:611406/1‑148 (MQ=255)
                tttGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCg                                                                                                                              >  2:745704/1‑148 (MQ=255)
                     aTAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCActg                                                                                                                         <  1:2390486/148‑1 (MQ=255)
                     aTAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCActg                                                                                                                         <  1:602271/148‑1 (MQ=255)
                      tAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCActgc                                                                                                                        <  1:39232/148‑1 (MQ=255)
                       aaaaaTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCActgct                                                                                                                       <  1:776053/148‑1 (MQ=255)
                           aTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTgg                                                                                                                   <  2:2111710/148‑1 (MQ=255)
                            ttCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGa                                                                                                                  >  1:3004251/1‑148 (MQ=255)
                              cAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATg                                                                                                                <  1:1439203/148‑1 (MQ=255)
                                 tGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGa                                                                                                             <  2:1502368/148‑1 (MQ=255)
                                   agagCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACa                                                                                                           >  1:1051497/1‑148 (MQ=255)
                                   agagCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACa                                                                                                           >  1:61642/1‑148 (MQ=255)
                                     agCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACaaa                                                                                                         >  1:147088/1‑148 (MQ=255)
                                     agCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACaaa                                                                                                         <  1:2114106/148‑1 (MQ=255)
                                       cGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAg                                                                                                       >  1:1958163/1‑148 (MQ=255)
                                       cGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAg                                                                                                       >  1:1193338/1‑148 (MQ=255)
                                        ggAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAgg                                                                                                      >  2:2905543/1‑148 (MQ=255)
                                         gAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAggg                                                                                                     <  1:456637/148‑1 (MQ=255)
                                          aaaTCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGa                                                                                                    >  2:1212001/1‑148 (MQ=255)
                                          aaaTCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGa                                                                                                    <  2:2925723/148‑1 (MQ=255)
                                           aaTCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAg                                                                                                   <  2:378773/148‑1 (MQ=255)
                                                 ggCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTc                                                                                             <  2:2412994/148‑1 (MQ=255)
                                                 ggCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTc                                                                                             <  2:692739/148‑1 (MQ=255)
                                                 ggCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTc                                                                                             <  2:888533/148‑1 (MQ=255)
                                                 ggCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTc                                                                                             <  1:1146310/148‑1 (MQ=255)
                                                  gCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTcc                                                                                            >  1:1109585/1‑148 (MQ=255)
                                                  gCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTcc                                                                                            >  1:1490668/1‑148 (MQ=255)
                                                  gCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTcc                                                                                            >  1:3084161/1‑148 (MQ=255)
                                                     catcatcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTa                                                                                         >  1:3001836/1‑148 (MQ=255)
                                                       tcatcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATg                                                                                       >  2:3026498/1‑148 (MQ=255)
                                                       tcatcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATg                                                                                       >  1:3201352/1‑148 (MQ=255)
                                                       tcatcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATg                                                                                       <  2:1261996/148‑1 (MQ=255)
                                                          tcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCat                                                                                    >  1:1453432/1‑148 (MQ=255)
                                                          tcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCat                                                                                    >  1:3095596/1‑148 (MQ=255)
                                                             gTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCAtata                                                                                 >  2:182028/1‑148 (MQ=255)
                                                             gTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCAtata                                                                                 <  1:1699665/148‑1 (MQ=255)
                                                             gTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCAtata                                                                                 >  2:2655716/1‑148 (MQ=255)
                                                              ttaattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATAc                                                                                >  1:2374115/1‑148 (MQ=255)
                                                              ttaattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATAc                                                                                >  2:1205992/1‑148 (MQ=255)
                                                                aattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTg                                                                              <  2:3080962/148‑1 (MQ=255)
                                                                  ttaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCa                                                                            >  2:2401384/1‑148 (MQ=255)
                                                                   taaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCaa                                                                           <  2:2472323/148‑1 (MQ=255)
                                                                      gCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAAcc                                                                        <  2:2371704/148‑1 (MQ=255)
                                                                      gCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAAcc                                                                        <  2:1377772/148‑1 (MQ=255)
                                                                      gCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAAcc                                                                        <  2:1537327/148‑1 (MQ=255)
                                                                        aGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTg                                                                      <  1:1844061/148‑1 (MQ=255)
                                                                          ggTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTggg                                                                    >  1:2132397/1‑148 (MQ=255)
                                                                           gtgtTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgc                                                                   <  1:517079/148‑1 (MQ=255)
                                                                            tgtTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgcc                                                                  <  2:2083533/148‑1 (MQ=255)
                                                                               tATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccagc                                                               <  1:2798770/148‑1 (MQ=255)
                                                                               tATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccagc                                                               <  1:2699061/148‑1 (MQ=255)
                                                                               tATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccagc                                                               >  2:1904235/1‑148 (MQ=255)
                                                                                   ttATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGa                                                           <  1:1212001/148‑1 (MQ=255)
                                                                                   ttATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGa                                                           <  1:2794318/148‑1 (MQ=255)
                                                                                      tGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGa                                                        <  2:2545208/148‑1 (MQ=255)
                                                                                        acgacgATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAag                                                      <  2:3023257/148‑1 (MQ=255)
                                                                                          gacgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAagag                                                    <  2:201909/148‑1 (MQ=255)
                                                                                           acgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAgg                                                   <  2:1563411/148‑1 (MQ=255)
                                                                                           acgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAgg                                                   >  1:2567906/1‑148 (MQ=255)
                                                                                           acgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAgg                                                   <  1:1902440/148‑1 (MQ=255)
                                                                                           acgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTAACGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAgg                                                   <  1:1282325/148‑1 (MQ=255)
                                                                                            cgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGa                                                  <  1:1991695/148‑1 (MQ=255)
                                                                                            cgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGa                                                  <  1:191642/148‑1 (MQ=255)
                                                                                              aTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGAtt                                                <  2:147088/148‑1 (MQ=255)
                                                                                                    aaGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTatgatg                                          <  1:3287895/148‑1 (MQ=255)
                                                                                                      gCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCg                                        >  1:2577405/1‑148 (MQ=255)
                                                                                                       cATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGa                                       >  1:913582/1‑148 (MQ=255)
                                                                                                         tctcCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATc                                     <  1:1341184/148‑1 (MQ=255)
                                                                                                         tctcCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATc                                     <  1:1649310/148‑1 (MQ=255)
                                                                                                          ctcCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATTCATATACTGCAAAACTGGGCCAGCCAGACGAAGAGGATTATGATGCGATcc                                    >  1:3062394/1‑148 (MQ=255)
                                                                                                          ctcCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATcc                                    >  2:1937613/1‑148 (MQ=255)
                                                                                                            cccGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCt                                  >  2:2305736/1‑148 (MQ=255)
                                                                                                             ccGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtc                                 >  2:1298691/1‑148 (MQ=255)
                                                                                                              cGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtcg                                <  1:2477859/148‑1 (MQ=255)
                                                                                                              cGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtcg                                <  1:819337/148‑1 (MQ=255)
                                                                                                                 tAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtcgtcg                             <  1:2273633/148‑1 (MQ=255)
                                                                                                                   ggTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTg                           >  1:961840/1‑148 (MQ=255)
                                                                                                                   ggTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTg                           >  1:2293054/1‑148 (MQ=255)
                                                                                                                         cGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATgg                     >  2:825802/1‑148 (MQ=255)
                                                                                                                            aTTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAAt                  <  1:3036501/148‑1 (MQ=255)
                                                                                                                             ttGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATa                 >  2:3251971/1‑148 (MQ=255)
                                                                                                                             ttGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATa                 >  2:1614704/1‑148 (MQ=255)
                                                                                                                             ttGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATa                 <  1:1190280/148‑1 (MQ=255)
                                                                                                                              tGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATAc                <  1:1859671/148‑1 (MQ=255)
                                                                                                                                gTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACgg              <  1:2036957/148‑1 (MQ=255)
                                                                                                                                   tCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGgcgc           >  1:374082/1‑148 (MQ=255)
                                                                                                                                    cGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGcgcg          >  2:2444342/1‑148 (MQ=255)
                                                                                                                                       tttttCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGgag       <  1:648042/148‑1 (MQ=255)
                                                                                                                                        ttttCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGgaga      <  1:1723570/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
ATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGGAGAACGC  >  CP009273/3311911‑3312198

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