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,1290CA97.7% 740.0 / 6.2 234A4S (GCA→TCA) ygcQputative flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base A (111/117);  total (113/121)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 236 total aligned reads displayed.

CACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCG  >  CP009273/2888984‑2889273
                                                                                                                                                 |                                                                                                                                                
 aCTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGat                                                                                                                                               <  1:1067338/148‑1 (MQ=255)
 aCTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGat                                                                                                                                               <  2:1129280/148‑1 (MQ=255)
  cTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGata                                                                                                                                              >  2:1015634/1‑148 (MQ=255)
    gTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATAtt                                                                                                                                            >  2:1654924/1‑148 (MQ=255)
     ttCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTc                                                                                                                                           >  1:1326096/1‑148 (MQ=255)
      tCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCa                                                                                                                                          >  2:3566595/1‑148 (MQ=255)
         accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATgg                                                                                                                                       <  1:829377/148‑1 (MQ=255)
         accaACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATgg                                                                                                                                       <  2:481481/148‑1 (MQ=255)
              cGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCa                                                                                                                                  >  2:2183635/1‑148 (MQ=255)
               gCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCAt                                                                                                                                 >  1:3475402/1‑148 (MQ=255)
                  tCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCt                                                                                                                              <  2:407793/148‑1 (MQ=255)
                  tCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCt                                                                                                                              >  2:230276/1‑148 (MQ=255)
                  tCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCt                                                                                                                              <  1:950602/148‑1 (MQ=255)
                  tCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCt                                                                                                                              >  2:1918562/1‑148 (MQ=255)
                   ccAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCtt                                                                                                                             >  1:2044894/1‑148 (MQ=255)
                        aCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCg                                                                                                                        <  2:1206057/148‑1 (MQ=255)
                         ccTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGt                                                                                                                       >  1:3785515/1‑148 (MQ=255)
                          cTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTa                                                                                                                      >  1:3430989/1‑148 (MQ=255)
                           tGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTaa                                                                                                                     <  1:225129/148‑1 (MQ=255)
                           tGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTaa                                                                                                                     >  1:1948587/1‑148 (MQ=255)
                           tGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTaa                                                                                                                     >  1:2763262/1‑148 (MQ=255)
                            gTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTaaa                                                                                                                    <  2:1077282/148‑1 (MQ=255)
                            gTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTaaa                                                                                                                    <  2:1077242/148‑1 (MQ=255)
                            gTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTaaa                                                                                                                    >  1:61189/1‑148 (MQ=255)
                              tCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATa                                                                                                                  <  2:1519958/148‑1 (MQ=255)
                              tCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATa                                                                                                                  <  1:1808221/148‑1 (MQ=255)
                                     cTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGc                                                                                                           >  2:2494252/1‑148 (MQ=255)
                                          ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                      <  1:923978/148‑1 (MQ=255)
                                          ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                      >  2:418287/1‑148 (MQ=255)
                                          ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                      >  2:2779300/1‑148 (MQ=255)
                                          ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                      >  2:2782081/1‑148 (MQ=255)
                                          ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                      >  2:3205993/1‑148 (MQ=255)
                                          ggCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAg                                                                                                      <  1:2150795/148‑1 (MQ=255)
                                             tGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCtt                                                                                                   <  2:2963849/148‑1 (MQ=255)
                                               aGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCtttt                                                                                                 <  2:244050/148‑1 (MQ=255)
                                                ggTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTg                                                                                                <  1:2781570/148‑1 (MQ=255)
                                                   tCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGc                                                                                             <  1:607197/148‑1 (MQ=255)
                                                     aaTCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCtt                                                                                           <  2:581626/148‑1 (MQ=255)
                                                        cTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCttttt                                                                                        >  2:1482305/1‑148 (MQ=255)
                                                         tGCCAGTGCGCCAGGTTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATTTTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTc                                                                                       <  1:3096388/148‑1 (MQ=255)
                                                         tGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTc                                                                                       <  1:1192876/148‑1 (MQ=255)
                                                            cAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGc                                                                                    <  1:2855495/148‑1 (MQ=255)
                                                              gTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTa                                                                                  <  1:1601757/148‑1 (MQ=255)
                                                              gTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCCGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTa                                                                                  <  2:2526486/148‑1 (MQ=255)
                                                                  gcCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGt                                                                              <  1:374154/148‑1 (MQ=255)
                                                                    cAGGGTACAACCGCTGAAATCCTGCGCCGCCAGACAGCTTGCTATTGCGGCATTCTCCTGATTTTTTGTCACAATTGAGATATTCATGGCTGCATCCTTTTGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGttt                                                                            <  2:1905917/148‑1 (MQ=255)
                                                                      gggTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGc                                                                          >  1:259033/1‑148 (MQ=255)
                                                                            aaCCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatca                                                                    <  1:1593820/148‑1 (MQ=255)
                                                                             aCCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaa                                                                   <  2:1729475/148‑1 (MQ=255)
                                                                              ccGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCatcaat                                                                  <  2:647592/148‑1 (MQ=255)
                                                                                  tGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGc                                                              >  1:1475228/1‑148 (MQ=255)
                                                                                   gAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCg                                                             <  1:546754/148‑1 (MQ=255)
                                                                                    aaaTCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGt                                                            <  1:3566595/148‑1 (MQ=255)
                                                                                       tCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGc                                                         >  2:215618/1‑148 (MQ=255)
                                                                                       tCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGc                                                         >  2:549845/1‑148 (MQ=255)
                                                                                        ccTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCt                                                        >  1:2233007/1‑148 (MQ=255)
                                                                                         cTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTc                                                       <  2:2590306/148‑1 (MQ=255)
                                                                                          tGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTAGCGTTTGCCCATCAATCAGCGTCGCTcc                                                      >  1:812342/1‑148 (MQ=255)
                                                                                             gccgccAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCcgcg                                                   <  2:1321044/148‑1 (MQ=255)
                                                                                              ccgccAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGcgt                                                  <  2:1103877/148‑1 (MQ=255)
                                                                                                gccagccaGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGcgtcg                                                <  2:644478/148‑1 (MQ=255)
                                                                                                   agccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtg                                             <  1:2907281/148‑1 (MQ=255)
                                                                                                    gccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgc                                            >  2:1520787/1‑148 (MQ=255)
                                                                                                    gccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgc                                            <  1:2021943/148‑1 (MQ=255)
                                                                                                     ccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgct                                           >  1:727415/1‑148 (MQ=255)
                                                                                                     ccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgct                                           >  1:2054867/1‑148 (MQ=255)
                                                                                                     ccagcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgct                                           <  1:2582622/148‑1 (MQ=255)
                                                                                                       agcTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctct                                         >  1:1002921/1‑148 (MQ=255)
                                                                                                          ttGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTtgctctgct                                      >  1:181933/1‑148 (MQ=255)
                                                                                                            gCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTgctctgctct                                    <  2:3770873/148‑1 (MQ=255)
                                                                                                                ttGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGcc                                >  1:674925/1‑148 (MQ=255)
                                                                                                                ttGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGcc                                >  2:111282/1‑148 (MQ=255)
                                                                                                                ttGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGcc                                >  2:3114011/1‑148 (MQ=255)
                                                                                                                ttGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGcc                                <  1:2627968/148‑1 (MQ=255)
                                                                                                                 tGCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTACTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCa                               >  2:1988295/1‑148 (MQ=255)
                                                                                                                  gCGGCATTCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAg                              <  1:3708853/148‑1 (MQ=255)
                                                                                                                   cGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGt                             <  1:2381102/148‑1 (MQ=255)
                                                                                                                        ttCTCCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAg                        <  2:1581920/148‑1 (MQ=255)
                                                                                                                         tctcCTGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGa                       >  1:1141039/1‑148 (MQ=255)
                                                                                                                              tGATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTg                  >  2:1797568/1‑148 (MQ=255)
                                                                                                                               gATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGc                 >  1:1641574/1‑148 (MQ=255)
                                                                                                                               gATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGc                 >  1:1872595/1‑148 (MQ=255)
                                                                                                                               gATTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGc                 >  2:1765992/1‑148 (MQ=255)
                                                                                                                                aTTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCa                <  2:3475402/148‑1 (MQ=255)
                                                                                                                                aTTTATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTg                  >  2:753567/1‑146 (MQ=255)
                                                                                                                                  ttATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc              >  1:1632584/1‑148 (MQ=255)
                                                                                                                                  ttATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc              >  1:3699281/1‑148 (MQ=255)
                                                                                                                                  ttATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATc              >  2:3138704/1‑148 (MQ=255)
                                                                                                                                   tATTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCg             <  1:2064633/148‑1 (MQ=255)
                                                                                                                                    aTTGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGc            <  2:1414291/148‑1 (MQ=255)
                                                                                                                                     ttGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCt           <  1:3161151/148‑1 (MQ=255)
                                                                                                                                     ttGTCACAATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCt           >  2:1773068/1‑148 (MQ=255)
                                                                                                                                         cacaATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCa       <  2:1948587/148‑1 (MQ=255)
                                                                                                                                         cacaATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCa       >  2:466238/1‑148 (MQ=255)
                                                                                                                                         cacaATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCa       <  1:549845/148‑1 (MQ=255)
                                                                                                                                          acaATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAt      >  1:2522828/1‑148 (MQ=255)
                                                                                                                                           caATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAtt     <  1:230276/148‑1 (MQ=255)
                                                                                                                                           caATTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAtt     >  2:709475/1‑148 (MQ=255)
                                                                                                                                            aaTTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAttt    <  1:1697303/148‑1 (MQ=255)
                                                                                                                                            aaTTGAGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCAttt    >  1:1379523/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTTGCTATTGCGGCATTCTCCTGATTTATTGTCACAATTGCGATATTCATGGCTGCATCCTTTGGCGTAAATAGTCCTGCCATAGCTTTTGAGCTTTTTCTGCTACCGTTTGCCCATCAATCAGCGTCGCTCCGCGTCGTTGCTCTGCTCTTGCCAGTTGCAGACACTGCATCGCTGGCATTTCG  >  CP009273/2888984‑2889273

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