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,0550CA96.6% 626.1 / 11.5 204G20G (GGC→GGA) argGargininosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base C (1/6);  new base A (108/89);  total (109/95)
Fisher's exact test for biased strand distribution p-value = 5.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
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 204 total aligned reads displayed.

AATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGGAGAACGC  >  CP009273/3311910‑3312198
                                                                                                                                                 |                                                                                                                                               
  tGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTc                                                                                                                                             >  2:18010/1‑148 (MQ=255)
   gAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCt                                                                                                                                            >  2:2438898/1‑148 (MQ=255)
     aaaCTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTgg                                                                                                                                          >  1:2202015/1‑148 (MQ=255)
      aaCTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGa                                                                                                                                         <  1:3163576/148‑1 (MQ=255)
       aCTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGac                                                                                                                                        >  2:2481974/1‑148 (MQ=255)
        cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGaca                                                                                                                                       <  1:2276037/148‑1 (MQ=255)
        cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGaca                                                                                                                                       >  1:1159231/1‑148 (MQ=255)
        cTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGaca                                                                                                                                       >  2:2063637/1‑148 (MQ=255)
          cATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACAcc                                                                                                                                     >  1:822450/1‑148 (MQ=255)
           atttatttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCa                                                                                                                                    <  1:363129/148‑1 (MQ=255)
            tttatttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAg                                                                                                                                   >  2:3118064/1‑148 (MQ=255)
             ttatttTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGt                                                                                                                                  >  2:519572/1‑148 (MQ=255)
                      aTAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCActg                                                                                                                         >  1:2788780/1‑148 (MQ=255)
                      aTAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCActg                                                                                                                         >  1:2789229/1‑148 (MQ=255)
                        aaaaaTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCActgct                                                                                                                       >  2:836349/1‑148 (MQ=255)
                         aaaaTTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCActgctg                                                                                                                      >  2:432100/1‑148 (MQ=255)
                             ttCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGa                                                                                                                  >  2:2897348/1‑148 (MQ=255)
                               cAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATg                                                                                                                <  2:261466/148‑1 (MQ=255)
                                  tGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGa                                                                                                             <  1:2469470/148‑1 (MQ=255)
                                   gagaGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGAc                                                                                                            >  1:264572/1‑148 (MQ=255)
                                     gagCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACaa                                                                                                          >  2:795950/1‑148 (MQ=255)
                                     gagCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACaa                                                                                                          >  1:1117264/1‑148 (MQ=255)
                                      agCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACaaa                                                                                                         >  1:2091632/1‑148 (MQ=255)
                                       gCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACaaaa                                                                                                        <  1:2468277/148‑1 (MQ=255)
                                       gCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACaaaa                                                                                                        <  1:2471539/148‑1 (MQ=255)
                                          gAAATACAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGTATGCGACAAAAggg                                                                                                     >  2:188977/1‑148 (MQ=255)
                                           aaaTCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGa                                                                                                    <  2:1601720/148‑1 (MQ=255)
                                            aaTCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAg                                                                                                   <  1:434675/148‑1 (MQ=255)
                                            aaTCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAg                                                                                                   <  2:344453/148‑1 (MQ=255)
                                              tCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCg                                                                                                 <  2:1564541/148‑1 (MQ=255)
                                                 aGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGtt                                                                                              >  2:3178129/1‑148 (MQ=255)
                                                  ggCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTc                                                                                             >  1:865015/1‑148 (MQ=255)
                                                   gCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTcc                                                                                            <  1:2288754/148‑1 (MQ=255)
                                                   gCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTcc                                                                                            <  2:1413787/148‑1 (MQ=255)
                                                    cTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCt                                                                                           <  1:2063637/148‑1 (MQ=255)
                                                     tcatcatcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCtt                                                                                          <  2:480461/148‑1 (MQ=255)
                                                       atcatcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTAt                                                                                        >  1:1305670/1‑148 (MQ=255)
                                                       atcatcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTAt                                                                                        >  1:1850347/1‑148 (MQ=255)
                                                        tcatcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATg                                                                                       >  2:2596224/1‑148 (MQ=255)
                                                        tcatcaGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATg                                                                                       >  1:3510675/1‑148 (MQ=255)
                                                           tcaGTTAATTAAGCAGGGTGTTATTTAATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCat                                                                                    >  2:174627/1‑148 (MQ=255)
                                                             aGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCatat                                                                                  <  2:2940858/148‑1 (MQ=255)
                                                              gTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCAtata                                                                                 >  2:3118202/1‑148 (MQ=255)
                                                               ttaattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATAc                                                                                <  1:2406886/148‑1 (MQ=255)
                                                               ttaattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATAc                                                                                <  1:519281/148‑1 (MQ=255)
                                                                taattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACt                                                                               >  1:109551/1‑148 (MQ=255)
                                                                 aattaaGCAGTGTGTTATTTTATGACTAAGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCTCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGATTTTGATGCGTCAAAAGGGATCGGTTCCTTATGCATATACTg                                                                              <  1:742158/148‑1 (MQ=255)
                                                                 aattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTg                                                                              <  1:742103/148‑1 (MQ=255)
                                                                 aattaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTGTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTg                                                                              <  2:921509/148‑1 (MQ=255)
                                                                   ttaaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCa                                                                            >  2:1933593/1‑148 (MQ=255)
                                                                     aaGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCaaa                                                                          <  1:2932873/148‑1 (MQ=255)
                                                                        cAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCt                                                                       <  1:2741389/148‑1 (MQ=255)
                                                                         aGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTg                                                                      <  1:317806/148‑1 (MQ=255)
                                                                          gggTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTgg                                                                     <  2:3367921/148‑1 (MQ=255)
                                                                             tgtTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgcc                                                                  <  2:240640/148‑1 (MQ=255)
                                                                               ttATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccag                                                                <  2:1194924/148‑1 (MQ=255)
                                                                               ttATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccag                                                                >  1:1533750/1‑148 (MQ=255)
                                                                               ttATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccag                                                                >  2:2195340/1‑148 (MQ=255)
                                                                                tATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccagc                                                               <  2:1209819/148‑1 (MQ=255)
                                                                                tATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccagc                                                               >  1:3340586/1‑148 (MQ=255)
                                                                                  ttttATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGgccagcca                                                             >  2:2077590/1‑148 (MQ=255)
                                                                                   tttATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGccagccag                                                            <  1:1695336/148‑1 (MQ=255)
                                                                                   tttATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGccagccag                                                            <  1:335315/148‑1 (MQ=255)
                                                                                        gacgacGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGaa                                                       >  1:361879/1‑148 (MQ=255)
                                                                                         acgacgATTCTCAAGCATCTCCCGTTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAag                                                      <  2:1517899/148‑1 (MQ=255)
                                                                                         acgacgATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAag                                                      <  1:836349/148‑1 (MQ=255)
                                                                                          cgacgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAaga                                                     <  1:588770/148‑1 (MQ=255)
                                                                                            acgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAgg                                                   <  1:1456157/148‑1 (MQ=255)
                                                                                            acgaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCGTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAgg                                                   >  2:2981327/1‑148 (MQ=255)
                                                                                              gaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGAt                                                 <  1:1019152/148‑1 (MQ=255)
                                                                                              gaTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGAt                                                 <  1:2104464/148‑1 (MQ=255)
                                                                                               aTTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGAtt                                                >  2:705913/1‑148 (MQ=255)
                                                                                                ttCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTa                                               >  1:1367015/1‑148 (MQ=255)
                                                                                                ttCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTa                                               >  2:1736529/1‑148 (MQ=255)
                                                                                                  ctcAAGCATCGCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTatg                                             <  1:2569457/148‑1 (MQ=255)
                                                                                                   tcAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTatga                                            <  2:2684219/148‑1 (MQ=255)
                                                                                                       gCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCg                                        >  2:2218336/1‑148 (MQ=255)
                                                                                                        cATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGa                                       <  1:416998/148‑1 (MQ=255)
                                                                                                        cATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGa                                       >  1:2973670/1‑148 (MQ=255)
                                                                                                          tctcCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATc                                     >  1:2531536/1‑148 (MQ=255)
                                                                                                           ctcCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATcc                                    <  2:822450/148‑1 (MQ=255)
                                                                                                             cccGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCt                                  >  1:2752959/1‑148 (MQ=255)
                                                                                                                  tAGGTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCtcgtcg                             >  1:3273713/1‑148 (MQ=255)
                                                                                                                    ggTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTg                           >  1:2330251/1‑148 (MQ=255)
                                                                                                                    ggTCAACGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTg                           >  1:544396/1‑148 (MQ=255)
                                                                                                                         aCGTATTGGTATCGCTTTTTCTGGAGGTCTTGACACCATTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGACAGCCAGACGAAGTGGATTATGATGCGATCCCTCGTCGTGCCATg                      >  2:442875/1‑148 (MQ=255)
                                                                                                                         aCGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATg                      >  2:2048016/1‑148 (MQ=255)
                                                                                                                          cGTATTTGTATCGCTTTTTCTGGATGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATAGACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATgg                     <  1:174627/148‑1 (MQ=255)
                                                                                                                          cGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATgg                     >  1:73861/1‑148 (MQ=255)
                                                                                                                          cGTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATgg                     <  1:1812791/148‑1 (MQ=255)
                                                                                                                           gTATTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGa                    <  1:1645438/148‑1 (MQ=255)
                                                                                                                             aTTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAAt                  <  2:664072/148‑1 (MQ=255)
                                                                                                                             aTTGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAAt                  >  1:2887942/1‑148 (MQ=255)
                                                                                                                              ttGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATa                 <  1:3263500/148‑1 (MQ=255)
                                                                                                                              ttGGTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATa                 >  1:3458519/1‑148 (MQ=255)
                                                                                                                                 gTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACgg              <  2:1710247/148‑1 (MQ=255)
                                                                                                                                 gTATCGCTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACgg              >  1:803111/1‑148 (MQ=255)
                                                                                                                                       cTTTTTCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGga        >  1:63681/1‑148 (MQ=255)
                                                                                                                                         ttttCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGgaga      >  1:2683111/1‑148 (MQ=255)
                                                                                                                                            tCTGGAGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGGAGAACg   >  2:25652/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AATGAAAACTCATTTATTTTGCATAAAAATTCAGTGAGAGCGGAAATCCAGGCTCATCATCAGTTAATTAAGCAGGGTGTTATTTTATGACGACGATTCTCAAGCATCTCCCGGTAGGTCAACGTATTGGTATCGCTTTTTCTGGCGGTCTGGACACCAGTGCCGCACTGCTGTGGATGCGACAAAAGGGAGCGGTTCCTTATGCATATACTGCAAACCTGGGCCAGCCAGACGAAGAGGATTATGATGCGATCCCTCGTCGTGCCATGGAATACGGCGCGGAGAACGC  >  CP009273/3311910‑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: ‑