Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,688,879 C→A 6.8% 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.8% 758.3 / 19.9 204pseudogene (467/753 nt)bcsQinvolved in cellulose production, minD superfamily (pseudogene)
Reads supporting (aligned to +/- strand):  ref base C (112/78);  new base A (11/3);  total (123/81)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 208 total aligned reads displayed.

TCGCTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTAAGCCGG  >  CP009273/3688733‑3689026
                                                                                                                                                  |                                                                                                                                                   
tCGCTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCg                                                                                                                                                    >  2:1668728/1‑148 (MQ=255)
  gCTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGa                                                                                                                                                  >  1:295205/1‑148 (MQ=255)
   cTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGat                                                                                                                                                 <  2:2504022/148‑1 (MQ=255)
   cTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGat                                                                                                                                                 >  2:736699/1‑148 (MQ=255)
   cTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGat                                                                                                                                                 <  2:614005/148‑1 (MQ=255)
   cTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGat                                                                                                                                                 >  2:599476/1‑148 (MQ=255)
    ttCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGata                                                                                                                                                <  2:3166547/148‑1 (MQ=255)
        tcaCGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGc                                                                                                                                            >  2:508154/1‑148 (MQ=255)
           cGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGt                                                                                                                                         <  1:1210174/148‑1 (MQ=255)
           cGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTAGGATATGGCAGt                                                                                                                                         >  1:2004204/1‑148 (MQ=255)
             atgaatgaGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTg                                                                                                                                       >  2:3077231/1‑148 (MQ=255)
              tgaatgaGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTAGGATATGGCAGTTgg                                                                                                                                      >  2:659791/1‑148 (MQ=255)
               gaatgaGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGc                                                                                                                                     <  1:3290105/148‑1 (MQ=255)
               gaatgaGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGc                                                                                                                                     >  1:2057820/1‑148 (MQ=255)
                aatgaGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCa                                                                                                                                    >  2:2373744/1‑148 (MQ=255)
                         cGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTg                                                                                                                           >  2:2120044/1‑148 (MQ=255)
                         cGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTg                                                                                                                           <  1:2923363/148‑1 (MQ=255)
                         cGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTg                                                                                                                           >  2:1445700/1‑148 (MQ=255)
                         cGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTAGGATATGGCAGTTGGCATCCACATTg                                                                                                                           >  1:3182466/1‑148 (MQ=255)
                          ggCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGa                                                                                                                          <  2:1611314/148‑1 (MQ=255)
                           gCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGAc                                                                                                                         >  2:2796681/1‑148 (MQ=255)
                           gCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGAc                                                                                                                         >  2:2747007/1‑148 (MQ=255)
                              gTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGAt                                                                                                                      <  2:1720224/148‑1 (MQ=255)
                                aaTCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTg                                                                                                                    <  1:829356/148‑1 (MQ=255)
                                 aTCGGCGCTGGCTTTGCAACCAAAGCTGGTAAGTATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGc                                                                                                                   >  2:1780929/1‑148 (MQ=255)
                                   cGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCa                                                                                                                 <  1:2507303/148‑1 (MQ=255)
                                    ggCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAg                                                                                                                >  1:469413/1‑148 (MQ=255)
                                      cgcTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCg                                                                                                              <  1:1540808/148‑1 (MQ=255)
                                      cgcTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCg                                                                                                              <  1:3183260/148‑1 (MQ=255)
                                       gcTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGa                                                                                                             <  1:3278716/148‑1 (MQ=255)
                                         tGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGt                                                                                                           <  2:2913333/148‑1 (MQ=255)
                                         tGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGt                                                                                                           <  2:1264438/148‑1 (MQ=255)
                                          ggCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTg                                                                                                          <  1:1771234/148‑1 (MQ=255)
                                          ggCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTg                                                                                                          <  1:2925439/148‑1 (MQ=255)
                                            cTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGAt                                                                                                        >  2:2478316/1‑148 (MQ=255)
                                            cTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGAt                                                                                                        >  1:823010/1‑148 (MQ=255)
                                            cTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGAt                                                                                                        <  2:788543/148‑1 (MQ=255)
                                                   aCCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACaa                                                                                                 <  1:2373744/148‑1 (MQ=255)
                                                   aCCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAAAATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACaa                                                                                                 >  2:887960/1‑148 (MQ=255)
                                                    ccAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACaaa                                                                                                >  2:1356665/1‑148 (MQ=255)
                                                    ccAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACaaa                                                                                                >  2:1356676/1‑148 (MQ=255)
                                                     cAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAAc                                                                                               >  1:239857/1‑148 (MQ=255)
                                                       aaGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTc                                                                                             >  2:2022873/1‑148 (MQ=255)
                                                        aGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTca                                                                                            >  2:2795862/1‑148 (MQ=255)
                                                         gCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTcag                                                                                           >  2:756332/1‑148 (MQ=255)
                                                          cTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTcagc                                                                                          >  2:1570530/1‑148 (MQ=255)
                                                            ggTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTcagcag                                                                                        >  2:380222/1‑148 (MQ=255)
                                                            ggTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTcagcag                                                                                        >  1:1454554/1‑148 (MQ=255)
                                                              tAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCt                                                                                      >  1:433358/1‑148 (MQ=255)
                                                                aaTATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCtgg                                                                                    <  1:332865/148‑1 (MQ=255)
                                                                aaTATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCtgg                                                                                    >  2:1936279/1‑148 (MQ=255)
                                                                 atatCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCtggt                                                                                   >  2:1001022/1‑148 (MQ=255)
                                                                  tatCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCtggtg                                                                                  <  1:2478165/148‑1 (MQ=255)
                                                                   atCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCtggtgg                                                                                 >  1:1071948/1‑148 (MQ=255)
                                                                   atCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCtggtgg                                                                                 >  2:2261192/1‑148 (MQ=255)
                                                                    tCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTggg                                                                                <  1:1051079/148‑1 (MQ=255)
                                                                     cGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGt                                                                               >  1:2624462/1‑148 (MQ=255)
                                                                       tCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTa                                                                             >  1:2722036/1‑148 (MQ=255)
                                                                           gAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTAGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTg                                                                         >  1:1900630/1‑148 (MQ=255)
                                                                             aCCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCg                                                                       <  1:1288562/148‑1 (MQ=255)
                                                                             aCCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCg                                                                       <  1:2604066/148‑1 (MQ=255)
                                                                              ccTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGa                                                                      >  1:86545/1‑148 (MQ=255)
                                                                              ccTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGa                                                                      >  2:2613013/1‑148 (MQ=255)
                                                                                  aCTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCa                                                                  <  1:1226482/148‑1 (MQ=255)
                                                                                     gCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCa                                                               <  1:2668495/148‑1 (MQ=255)
                                                                                             ggAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTaa                                                       <  1:394370/148‑1 (MQ=255)
                                                                                              gAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAg                                                      <  1:956781/148‑1 (MQ=255)
                                                                                              gAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAg                                                      <  1:203930/148‑1 (MQ=255)
                                                                                              gAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTAGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAg                                                      <  1:1104968/148‑1 (MQ=255)
                                                                                               aaGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGt                                                     <  1:1187364/148‑1 (MQ=255)
                                                                                               aaGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGt                                                     <  1:1363420/148‑1 (MQ=255)
                                                                                               aaGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGt                                                     <  1:3017563/148‑1 (MQ=255)
                                                                                                    aTTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGAtt                                                >  2:2732958/1‑148 (MQ=255)
                                                                                                            aaaTATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCa                                        >  1:2089412/1‑148 (MQ=255)
                                                                                                            aaaTATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCa                                        >  1:2089367/1‑148 (MQ=255)
                                                                                                              atatGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACt                                      >  1:2773133/1‑148 (MQ=255)
                                                                                                              atatGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACt                                      >  1:380365/1‑148 (MQ=255)
                                                                                                                  gtgCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTa                                  >  2:1433111/1‑148 (MQ=255)
                                                                                                                     cgcCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACg                               <  2:818685/148‑1 (MQ=255)
                                                                                                                     cgcCATCAGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACTTTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACg                               <  2:3033167/148‑1 (MQ=255)
                                                                                                                       ccATCCGGCAGCTCTTGCTGATGCAGTCGGATATGGCAGTTGGAATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGcc                             <  1:1234757/148‑1 (MQ=255)
                                                                                                                       ccATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGcc                             >  2:44935/1‑148 (MQ=255)
                                                                                                                         aTCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCg                           >  2:3279795/1‑148 (MQ=255)
                                                                                                                          tCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGc                          >  1:548244/1‑148 (MQ=255)
                                                                                                                           ccGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCt                         >  2:741645/1‑148 (MQ=255)
                                                                                                                           ccGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCt                         >  2:1273092/1‑148 (MQ=255)
                                                                                                                               cAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCt                     >  1:836011/1‑148 (MQ=255)
                                                                                                                                  cgcTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCtttt                  >  2:2258183/1‑148 (MQ=255)
                                                                                                                                    cTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAg                <  2:2409208/148‑1 (MQ=255)
                                                                                                                                     ttGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGt               <  2:1020974/148‑1 (MQ=255)
                                                                                                                                          gATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTg          <  1:2709434/148‑1 (MQ=255)
                                                                                                                                           aTGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGt         >  2:403146/1‑148 (MQ=255)
                                                                                                                                           aTGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGt         >  2:403101/1‑148 (MQ=255)
                                                                                                                                           aTGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGt         >  2:679947/1‑148 (MQ=255)
                                                                                                                                             gCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTaa       >  2:1174308/1‑148 (MQ=255)
                                                                                                                                                gTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTAAGcc    >  1:169712/1‑148 (MQ=255)
                                                                                                                                                 tCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTAAGCCg   >  1:1120320/1‑148 (MQ=255)
                                                                                                                                                 tCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTAAGCCg   >  2:1666949/1‑148 (MQ=255)
                                                                                                                                                  cGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTAAGCCgg  <  1:2261164/148‑1 (MQ=255)
                                                                                                                                                  cGGATAGGGCAGTGGGCCACCACAGTGACGATTGCCAGCGAGGGATCGCACAAACACCGCAGCGGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTAAGCCgg  <  1:2261354/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TCGCTTCATCACGATGAATGAGCATCGGCAGTAATCGGCGCTGGCTTTGCAACCAAAGCTGGTAAATATCGTCCTGAACCTGACTGCCAATACGGAAGTCATTAATCAAAATATGTGCGCCATCCGGCAGCGCTTGCTGATGCAGTCGGATATGGCAGTTGGCATCCACATTGACGATTGCCAGCGAGTGATCGCACAAACTCAGCAGCTGGTGGGTTATCTGCGAGGCATCACGCGGTAAGTCGATTAAAATCCACTGGTAACGCCCGCTGGCTTTTAGTTGCTGTAAGCCGG  >  CP009273/3688733‑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: