Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,889,129 C→A 100% A4S (GCA→TCA)  ygcQ ← putative flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,889,1290CA98.3% 580.9 / 0.8 181A4S (GCA→TCA) ygcQputative flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (94/83);  total (96/85)
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.97e-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 183 total aligned reads displayed.

CACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAT  >  CP009273/2888984‑2889269
                                                                                                                                                 |                                                                                                                                            
cACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTTCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTgaga                                                                                                                                            <  1:2001181/148‑1 (MQ=255)
 aCTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGat                                                                                                                                           <  2:2082683/148‑1 (MQ=255)
   tGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCACTCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGatat                                                                                                                                         <  2:2190907/148‑1 (MQ=255)
   tGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGatat                                                                                                                                         <  2:1878960/148‑1 (MQ=255)
      tCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCa                                                                                                                                      <  2:857150/148‑1 (MQ=255)
       ccaccaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCAt                                                                                                                                     <  2:295669/148‑1 (MQ=255)
       ccaccaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCAt                                                                                                                                     >  1:3411164/1‑148 (MQ=255)
       ccaccaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCAt                                                                                                                                     >  2:3025454/1‑148 (MQ=255)
         accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATgg                                                                                                                                   <  1:1314274/148‑1 (MQ=255)
          ccaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGc                                                                                                                                  <  2:2275337/148‑1 (MQ=255)
            aaCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTg                                                                                                                                <  2:368514/148‑1 (MQ=255)
             aCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGc                                                                                                                               >  1:254255/1‑148 (MQ=255)
               gCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCAt                                                                                                                             <  1:3322004/148‑1 (MQ=255)
                cATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATc                                                                                                                            <  1:2933334/148‑1 (MQ=255)
                cATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATc                                                                                                                            <  2:501508/148‑1 (MQ=255)
                cATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATc                                                                                                                            <  1:1644686/148‑1 (MQ=255)
                 aTCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATcc                                                                                                                           <  2:146341/148‑1 (MQ=255)
                 aTCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATcc                                                                                                                           <  2:146370/148‑1 (MQ=255)
                  tCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCt                                                                                                                          >  2:3342930/1‑148 (MQ=255)
                  tCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCt                                                                                                                          >  1:3254366/1‑148 (MQ=255)
                    cAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCttt                                                                                                                        >  1:2347517/1‑148 (MQ=255)
                     aGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTg                                                                                                                       >  1:2993494/1‑148 (MQ=255)
                     aGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTg                                                                                                                       >  2:2961184/1‑148 (MQ=255)
                      gCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTgg                                                                                                                      >  2:2973328/1‑148 (MQ=255)
                         ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGt                                                                                                                   >  1:2197230/1‑148 (MQ=255)
                           tGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTaa                                                                                                                 <  1:418842/148‑1 (MQ=255)
                            gTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTaaa                                                                                                                >  2:2202317/1‑148 (MQ=255)
                                  ccACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTcc                                                                                                          <  2:1019587/148‑1 (MQ=255)
                                     cTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGc                                                                                                       >  1:2108488/1‑148 (MQ=255)
                                     cTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGc                                                                                                       >  1:1222243/1‑148 (MQ=255)
                                         aGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATa                                                                                                   <  1:2153821/148‑1 (MQ=255)
                                          ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                  <  2:2308172/148‑1 (MQ=255)
                                             tGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCtt                                                                                               <  2:980674/148‑1 (MQ=255)
                                              gAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCttt                                                                                              >  2:418252/1‑148 (MQ=255)
                                               aGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCtttt                                                                                             >  1:1176547/1‑148 (MQ=255)
                                                ggTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTg                                                                                            <  1:2903055/148‑1 (MQ=255)
                                                 gTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGa                                                                                           <  2:1085112/148‑1 (MQ=255)
                                                    cAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCt                                                                                        <  2:1291821/148‑1 (MQ=255)
                                                    cAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCt                                                                                        >  1:2579757/1‑148 (MQ=255)
                                                      aTCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttt                                                                                      >  1:1332066/1‑148 (MQ=255)
                                                       tCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCtttt                                                                                     >  1:760348/1‑148 (MQ=255)
                                                        cTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttttt                                                                                    <  2:2812662/148‑1 (MQ=255)
                                                        cTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttttt                                                                                    >  1:1826580/1‑148 (MQ=255)
                                                          gCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCt                                                                                  >  2:2046221/1‑148 (MQ=255)
                                                                gcgcCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTAcc                                                                            <  2:1236065/148‑1 (MQ=255)
                                                                     aGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTg                                                                       >  2:523131/1‑148 (MQ=255)
                                                                      gggTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGc                                                                      >  1:706734/1‑148 (MQ=255)
                                                                       ggTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGcc                                                                     >  1:392479/1‑148 (MQ=255)
                                                                       ggTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGcc                                                                     >  1:392444/1‑148 (MQ=255)
                                                                          aCAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCat                                                                  >  2:1910349/1‑148 (MQ=255)
                                                                          aCAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCat                                                                  <  2:2791224/148‑1 (MQ=255)
                                                                          aCAACCGCTGAAATCCGGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCat                                                                  >  2:45054/1‑148 (MQ=255)
                                                                           cAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatc                                                                 >  2:1770339/1‑148 (MQ=255)
                                                                           cAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatc                                                                 >  1:2904391/1‑148 (MQ=255)
                                                                              ccGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaat                                                              <  1:2418758/148‑1 (MQ=255)
                                                                               cGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaatc                                                             <  2:1826580/148‑1 (MQ=255)
                                                                               cGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAAGTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaatc                                                             >  2:1520403/1‑148 (MQ=255)
                                                                                    aaaTCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGt                                                        <  1:1167018/148‑1 (MQ=255)
                                                                                    aaaTCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGt                                                        >  2:2975482/1‑148 (MQ=255)
                                                                                              ccgccAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTACATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGcgt                                              <  2:3498904/148‑1 (MQ=255)
                                                                                               cgccAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTAGTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGcgtc                                             <  2:3049585/148‑1 (MQ=255)
                                                                                                gccagccaGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGcgtcg                                            >  1:2291711/1‑148 (MQ=255)
                                                                                                gccagccaGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGcgtcg                                            <  2:2320415/148‑1 (MQ=255)
                                                                                                 ccagccagCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGcgtcgt                                           <  1:2501925/148‑1 (MQ=255)
                                                                                                  cagccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAACTAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGtt                                          >  2:3255041/1‑148 (MQ=255)
                                                                                                     ccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgct                                       >  2:2571205/1‑148 (MQ=255)
                                                                                                     ccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgct                                       >  2:2571100/1‑148 (MQ=255)
                                                                                                        gcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctctg                                    <  2:2331578/148‑1 (MQ=255)
                                                                                                        gcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctctg                                    >  1:2326450/1‑148 (MQ=255)
                                                                                                         cTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctctgc                                   >  1:366816/1‑148 (MQ=255)
                                                                                                          ttGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctctgct                                  <  2:3086024/148‑1 (MQ=255)
                                                                                                           tGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctctgctc                                 <  1:476554/148‑1 (MQ=255)
                                                                                                                ttGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGcc                            >  2:2957936/1‑148 (MQ=255)
                                                                                                                 tGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCa                           >  1:1880789/1‑148 (MQ=255)
                                                                                                                   cGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGt                         <  1:1079670/148‑1 (MQ=255)
                                                                                                                      cATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGc                      >  1:80586/1‑148 (MQ=255)
                                                                                                                       aTTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCa                     >  1:2010404/1‑148 (MQ=255)
                                                                                                                         tctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGa                   <  2:2993494/148‑1 (MQ=255)
                                                                                                                         tctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGa                   >  2:3165660/1‑148 (MQ=255)
                                                                                                                          ctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGac                  >  1:1299155/1‑148 (MQ=255)
                                                                                                                          ctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCAATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGac                  >  2:282262/1‑148 (MQ=255)
                                                                                                                          ctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCAATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGac                  >  2:282168/1‑148 (MQ=255)
                                                                                                                           tcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGaca                 <  1:3209915/148‑1 (MQ=255)
                                                                                                                              tGATTTATTGTCACAAGTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTg              <  2:1661935/148‑1 (MQ=255)
                                                                                                                                aTTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCa            <  2:2805608/148‑1 (MQ=255)
                                                                                                                                aTTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCa            >  1:2614880/1‑148 (MQ=255)
                                                                                                                                aTTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCa            <  1:3193412/148‑1 (MQ=255)
                                                                                                                                aTTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCa            >  2:3502266/1‑148 (MQ=255)
                                                                                                                                aTTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCa            <  2:3468099/148‑1 (MQ=255)
                                                                                                                                  ttATTTACACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCAAAGCTTATGAGCATTTTCTCATACCGTTTGCCCATCAATCCGCGTCGCAACGCGTCGTTGATCTGATCTTGCCAGATGCAGACAATCCATc          >  2:443363/1‑148 (MQ=255)
                                                                                                                                  ttATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc          >  2:94692/1‑148 (MQ=255)
                                                                                                                                  ttATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc          >  1:3081185/1‑148 (MQ=255)
                                                                                                                                  ttATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc          >  2:3253017/1‑148 (MQ=255)
                                                                                                                                  ttATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc          >  1:2238547/1‑148 (MQ=255)
                                                                                                                                     ttGTCACATTTGTTATATTGATGTCTTCATGCTTTTGCGTAAATAGTACTGACATAGCTTTTGAGATTTTTCTGCTAGCGTTTGCCCATCAATCAGCGTGGCTCCGCGTCGTTGCTCTGCTCTTGACAGTTGCAGTCACTTCATCGCt       <  2:3274018/148‑1 (MQ=255)
                                                                                                                                      tGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTg      <  1:102957/148‑1 (MQ=255)
                                                                                                                                      tGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTg      >  2:2689810/1‑148 (MQ=255)
                                                                                                                                         cacaATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCa   >  2:741267/1‑148 (MQ=255)
                                                                                                                                         cacaATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCa   >  2:3491405/1‑148 (MQ=255)
                                                                                                                                         cacaATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCa   >  2:1572422/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
CACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAT  >  CP009273/2888984‑2889269

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: