Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,688,879 C→A 6.1% pseudogene (467/753 nt) bcsQ ← involved in cellulose production, minD superfamily (pseudogene)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,688,8790CA 6.1% 732.3 / 16.8 195pseudogene (467/753 nt)bcsQinvolved in cellulose production, minD superfamily (pseudogene)
Reads supporting (aligned to +/- strand):  ref base C (91/92);  new base A (7/5);  total (98/97)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 198 total aligned reads displayed.

ATCGCTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTAAGCCGG  >  CP009273/3688732‑3689026
                                                                                                                                                   |                                                                                                                                                   
aTCGCTTCATCACGATGAATGAGCATCGTCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTc                                                                                                                                                     >  1:1437427/1‑148 (MQ=255)
 tCGCTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCg                                                                                                                                                    >  1:390665/1‑148 (MQ=255)
  cGCTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCgg                                                                                                                                                   >  2:3236063/1‑148 (MQ=255)
    cTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGat                                                                                                                                                 <  1:586608/148‑1 (MQ=255)
    cTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGat                                                                                                                                                 <  1:2741383/148‑1 (MQ=255)
    cTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGat                                                                                                                                                 >  1:2027430/1‑148 (MQ=255)
     ttCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGata                                                                                                                                                <  2:2954273/148‑1 (MQ=255)
     ttCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGata                                                                                                                                                <  2:1178746/148‑1 (MQ=255)
      tcatcaCGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGatat                                                                                                                                               <  1:699418/148‑1 (MQ=255)
         tcaCGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGc                                                                                                                                            >  2:1938668/1‑148 (MQ=255)
         tcaCGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGc                                                                                                                                            <  1:3493588/148‑1 (MQ=255)
             gATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGtt                                                                                                                                        >  1:1458562/1‑148 (MQ=255)
               tgaatgaGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTgg                                                                                                                                      <  1:1588121/148‑1 (MQ=255)
                gaatgaGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGc                                                                                                                                     >  2:2527283/1‑148 (MQ=255)
                  atgaGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCAt                                                                                                                                   >  1:3391824/1‑148 (MQ=255)
                    gaGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATcc                                                                                                                                 >  1:2040746/1‑148 (MQ=255)
                     aGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCca                                                                                                                                >  2:1973703/1‑148 (MQ=255)
                     aGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCca                                                                                                                                >  2:2602789/1‑148 (MQ=255)
                     aGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCca                                                                                                                                >  2:1154070/1‑148 (MQ=255)
                      gCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCcac                                                                                                                               >  1:2740770/1‑148 (MQ=255)
                       cATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCcaca                                                                                                                              <  1:1842782/148‑1 (MQ=255)
                       cATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCcaca                                                                                                                              <  2:1687577/148‑1 (MQ=255)
                       cATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTAGGATATGGCAGTTGGCATCcaca                                                                                                                              <  2:1694899/148‑1 (MQ=255)
                         tCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACAtt                                                                                                                            >  2:3356117/1‑148 (MQ=255)
                             cAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACATGACTGCCAATACGGAAGTCATTAATCAAAATAAGTGCGCCAAACGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTTGCATCCACATAGACg                                                                                                                        >  1:1966258/1‑148 (MQ=255)
                                tAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGAtt                                                                                                                     >  1:3030686/1‑148 (MQ=255)
                                tAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGAtt                                                                                                                     >  1:2657504/1‑148 (MQ=255)
                                 aaTCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTg                                                                                                                    >  2:1275909/1‑148 (MQ=255)
                                      gcgcTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTAGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGc                                                                                                               >  1:1069705/1‑148 (MQ=255)
                                       cgcTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCg                                                                                                              <  2:705448/148‑1 (MQ=255)
                                       cgcTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTAGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCg                                                                                                              <  2:1730640/148‑1 (MQ=255)
                                         cTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAg                                                                                                            >  1:261424/1‑148 (MQ=255)
                                         cTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTAGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAg                                                                                                            >  1:776393/1‑148 (MQ=255)
                                           ggCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTg                                                                                                          >  1:2084937/1‑148 (MQ=255)
                                            gCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGa                                                                                                         <  2:1883293/148‑1 (MQ=255)
                                             cTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCACCCATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGAt                                                                                                        <  2:297107/148‑1 (MQ=255)
                                                tGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGc                                                                                                     >  1:663080/1‑148 (MQ=255)
                                                      cAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAAc                                                                                               <  1:2020849/148‑1 (MQ=255)
                                                       aaaGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACt                                                                                              <  1:2210132/148‑1 (MQ=255)
                                                        aaGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTc                                                                                             <  1:480705/148‑1 (MQ=255)
                                                        aaGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTc                                                                                             <  1:3525708/148‑1 (MQ=255)
                                                                  atatCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTcagcag                                                                                        >  1:1637154/1‑143 (MQ=255)
                                                                  atatCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTcagcag                                                                                        <  2:1637154/143‑1 (MQ=255)
                                                                   tatCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCtggtg                                                                                  <  2:945661/148‑1 (MQ=255)
                                                                    atCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCtggtgg                                                                                 <  1:2604404/148‑1 (MQ=255)
                                                                    atCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCtggtgg                                                                                 >  2:1567359/1‑148 (MQ=255)
                                                                    atCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCtggtgg                                                                                 <  1:2985462/148‑1 (MQ=255)
                                                                     tCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTggg                                                                                >  2:2964862/1‑148 (MQ=255)
                                                                      cGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGt                                                                               >  1:871513/1‑148 (MQ=255)
                                                                        tCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTa                                                                             >  2:57298/1‑148 (MQ=255)
                                                                        tCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTa                                                                             >  2:3131466/1‑148 (MQ=255)
                                                                         ccTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTAt                                                                            <  1:1278172/148‑1 (MQ=255)
                                                                         ccTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTAt                                                                            <  1:236256/148‑1 (MQ=255)
                                                                         ccTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTAt                                                                            >  1:1006496/1‑148 (MQ=255)
                                                                         ccTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTAt                                                                            <  2:1324979/148‑1 (MQ=255)
                                                                         ccTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTAt                                                                            >  1:1733000/1‑148 (MQ=255)
                                                                            gAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTg                                                                         >  1:2159727/1‑148 (MQ=255)
                                                                                  gACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGc                                                                   <  1:3198738/148‑1 (MQ=255)
                                                                                   aCTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCa                                                                  >  1:2522104/1‑148 (MQ=255)
                                                                                     tGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATc                                                                <  1:2688870/148‑1 (MQ=255)
                                                                                      gCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCa                                                               >  1:1708104/1‑148 (MQ=255)
                                                                                      gCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCa                                                               >  1:3472222/1‑148 (MQ=255)
                                                                                        cAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCAcg                                                             <  2:2167929/148‑1 (MQ=255)
                                                                                         aaTACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCAcgc                                                            >  2:784066/1‑148 (MQ=255)
                                                                                          aTACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCTcgcg                                                           >  1:1227875/1‑148 (MQ=255)
                                                                                           tACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCgg                                                          >  1:1236516/1‑148 (MQ=255)
                                                                                           tACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCgg                                                          >  1:2251304/1‑148 (MQ=255)
                                                                                             cGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTa                                                        <  2:2522077/148‑1 (MQ=255)
                                                                                               gAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAg                                                      >  2:826007/1‑148 (MQ=255)
                                                                                                  gTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTAGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGa                                                                      <  2:2667911/129‑1 (MQ=255)
                                                                                                        aaTCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTaaa                                             >  2:2158459/1‑148 (MQ=255)
                                                                                                        aaTCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTaaa                                             <  1:1938668/148‑1 (MQ=255)
                                                                                                          tCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAAt                                           <  1:3431952/148‑1 (MQ=255)
                                                                                                          tCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTAGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAAt                                           >  1:1191218/1‑148 (MQ=255)
                                                                                                            aaaaTATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATcc                                         <  2:3310570/148‑1 (MQ=255)
                                                                                                             aaaTATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCa                                        <  2:2811175/148‑1 (MQ=255)
                                                                                                              aaTATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCAc                                       <  1:1401554/148‑1 (MQ=255)
                                                                                                              aaTATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCAc                                       <  2:871513/148‑1 (MQ=255)
                                                                                                                tatGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTg                                     >  2:1734852/1‑148 (MQ=255)
                                                                                                                tatGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTg                                     >  2:3497795/1‑148 (MQ=255)
                                                                                                                 atGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTgg                                    >  1:370156/1‑148 (MQ=255)
                                                                                                                 atGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTgg                                    <  2:2629869/148‑1 (MQ=255)
                                                                                                                      cgcCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACg                               <  2:2933021/148‑1 (MQ=255)
                                                                                                                       gcCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGc                              <  1:812718/148‑1 (MQ=255)
                                                                                                                        ccATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGcc                             >  2:3120069/1‑148 (MQ=255)
                                                                                                                        ccATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGcc                             <  2:3436332/148‑1 (MQ=255)
                                                                                                                          aTCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCg                           >  2:1257186/1‑148 (MQ=255)
                                                                                                                           tCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGc                          >  2:1348800/1‑148 (MQ=255)
                                                                                                                               gCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGc                      <  2:1238760/148‑1 (MQ=255)
                                                                                                                                cAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCt                     <  2:2027430/148‑1 (MQ=255)
                                                                                                                                  gcgcTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCttt                   <  1:1318290/148‑1 (MQ=255)
                                                                                                                                   cgcTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCtttt                  >  1:3457030/1‑148 (MQ=255)
                                                                                                                                     cTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAg                <  1:2054890/148‑1 (MQ=255)
                                                                                                                                     cTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAg                >  2:42328/1‑148 (MQ=255)
                                                                                                                                     cTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAg                >  1:3220412/1‑148 (MQ=255)
                                                                                                                                         cTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGc            >  2:275769/1‑148 (MQ=255)
                                                                                                                                          tGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCt           <  1:2327876/148‑1 (MQ=255)
                                                                                                                                             tGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTa        <  2:2188684/148‑1 (MQ=255)
                                                                                                                                              gCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTaa       <  2:609310/148‑1 (MQ=255)
                                                                                                                                                aGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTAAGc     >  2:429196/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCGCTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTAAGCCGG  >  CP009273/3688732‑3689026

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: