Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,878,098 T→A 78.2% A827A (GCA→GCT)  ygcB ← R‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,878,0980TA78.2% 497.3 / 137.7 229A827A (GCA→GCT) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base A (98/81);  total (124/105)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

CATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTAT  >  CP009273/2877950‑2878246
                                                                                                                                                    |                                                                                                                                                    
 aTCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGc                                                                                                                                                    >  1:2564213/1‑150 (MQ=255)
 aTCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGc                                                                                                                                                    >  1:96867/1‑150 (MQ=255)
 aTCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGc                                                                                                                                                    >  2:1919321/1‑150 (MQ=255)
   ccccTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGAGTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCaa                                                                                                                                                  <  2:3735585/150‑1 (MQ=255)
    cccTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAg                                                                                                                                                 <  2:879848/150‑1 (MQ=255)
            aGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCAt                                                                                                                                         >  1:2359026/1‑150 (MQ=255)
               taataaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATAct                                                                                                                                      <  1:3786828/150‑1 (MQ=255)
               taataaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATAct                                                                                                                                      <  1:3783183/150‑1 (MQ=255)
                aataaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATActg                                                                                                                                     <  2:3541298/150‑1 (MQ=255)
                    aCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCACAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGt                                                                                                                                 <  1:1262093/150‑1 (MQ=255)
                     cAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGtt                                                                                                                                <  2:2237183/150‑1 (MQ=255)
                     cAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGtt                                                                                                                                <  1:4176194/150‑1 (MQ=255)
                      aaTACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTc                                                                                                                               >  1:1809085/1‑150 (MQ=255)
                        tACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCAt                                                                                                                             <  2:2312163/150‑1 (MQ=255)
                         aCTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATg                                                                                                                            >  2:4447578/1‑150 (MQ=255)
                           tGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGAc                                                                                                                          >  1:4383330/1‑150 (MQ=255)
                               aCCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAg                                                                                                                      >  1:1067108/1‑150 (MQ=255)
                                 ccTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGt                                                                                                                    <  2:5014659/150‑1 (MQ=255)
                                 ccTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGt                                                                                                                    >  1:1557475/1‑150 (MQ=255)
                                    gCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCt                                                                                                                 <  2:1849471/150‑1 (MQ=255)
                                       aGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGt                                                                                                              >  2:807955/1‑150 (MQ=255)
                                       aGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGAGTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGt                                                                                                              >  2:259768/1‑150 (MQ=255)
                                       aGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACACGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGt                                                                                                              <  2:2576593/150‑1 (MQ=255)
                                        gACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTa                                                                                                             >  1:3044298/1‑150 (MQ=255)
                                        gACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTa                                                                                                             <  2:4637906/150‑1 (MQ=255)
                                        gACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTa                                                                                                             >  2:2289853/1‑150 (MQ=255)
                                         aCCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAg                                                                                                            >  1:4762122/1‑150 (MQ=255)
                                         aCCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAg                                                                                                            >  2:671418/1‑150 (MQ=255)
                                             atccatccAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCt                                                                                                        <  2:202783/150‑1 (MQ=255)
                                              tccatccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTg                                                                                                       >  1:1242923/1‑150 (MQ=255)
                                              tccatccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTg                                                                                                       >  1:407510/1‑150 (MQ=255)
                                              tccatccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTg                                                                                                       >  2:1201837/1‑150 (MQ=255)
                                              tccatccaGATTCTGTTTCCCTTCCAGCCAAAGGAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTg                                                                                                       <  2:2548244/150‑1 (MQ=255)
                                               ccatccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCCTGACTTAGGTCCTCGTAGACCTgg                                                                                                      <  1:4054803/150‑1 (MQ=255)
                                                     aGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCg                                                                                                <  2:51484/150‑1 (MQ=255)
                                                     aGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCg                                                                                                <  1:552787/150‑1 (MQ=255)
                                                       aTTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAg                                                                                              <  2:1723937/150‑1 (MQ=255)
                                                       aTTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAg                                                                                              >  2:4735412/1‑150 (MQ=255)
                                                        ttCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGc                                                                                             >  1:2248078/1‑150 (MQ=255)
                                                            gTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAgtt                                                                                         <  1:764351/150‑1 (MQ=255)
                                                                 ccTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTa                                                                                    >  1:2716642/1‑150 (MQ=255)
                                                                  cTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTAc                                                                                   >  1:2071505/1‑150 (MQ=255)
                                                                  cTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTAc                                                                                   <  1:4454684/150‑1 (MQ=255)
                                                                   ttCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTAcc                                                                                  >  2:4173144/1‑150 (MQ=255)
                                                                       aGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGaa                                                                              <  2:4273716/150‑1 (MQ=255)
                                                                       aGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGaa                                                                              <  2:1696186/150‑1 (MQ=255)
                                                                         ccAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGa                                                                            >  1:4341080/1‑150 (MQ=255)
                                                                           aaaGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACg                                                                          <  1:1373456/150‑1 (MQ=255)
                                                                            aaGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGt                                                                         >  2:3650899/1‑150 (MQ=255)
                                                                               tAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTg                                                                      >  1:4685838/1‑150 (MQ=255)
                                                                                  cccATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTAc                                                                   >  1:2810500/1‑150 (MQ=255)
                                                                                  cccATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTAc                                                                   <  2:3279871/150‑1 (MQ=255)
                                                                                  cccATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTAc                                                                   >  1:1603721/1‑150 (MQ=255)
                                                                                   ccATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACa                                                                  <  2:211149/150‑1 (MQ=255)
                                                                                    cATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACAt                                                                 <  2:1769911/150‑1 (MQ=255)
                                                                                       cTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAg                                                              <  1:2888171/150‑1 (MQ=255)
                                                                                       cTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAg                                                              <  2:4149761/150‑1 (MQ=255)
                                                                                       cTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAg                                                              <  2:2842053/150‑1 (MQ=255)
                                                                                        ttCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAgg                                                             >  1:4518970/1‑150 (MQ=255)
                                                                                         tCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGc                                                            >  1:67068/1‑150 (MQ=255)
                                                                                           aTCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaa                                                          >  2:2333026/1‑150 (MQ=255)
                                                                                            tCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaat                                                         <  1:1546168/150‑1 (MQ=255)
                                                                                            tCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaat                                                         >  1:4868984/1‑150 (MQ=255)
                                                                                               actactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaataat                                                      <  2:4015015/150‑1 (MQ=255)
                                                                                                 tactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATgg                                                    >  1:2562356/1‑150 (MQ=255)
                                                                                                  actTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGc                                                   >  1:1697905/1‑150 (MQ=255)
                                                                                                   ctTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCa                                                  >  1:3877453/1‑150 (MQ=255)
                                                                                                   ctTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCa                                                  <  1:1489572/150‑1 (MQ=255)
                                                                                                         aaaaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCa                                            <  2:1441123/150‑1 (MQ=255)
                                                                                                           aaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCAtt                                          <  1:3312337/150‑1 (MQ=255)
                                                                                                               tACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTccc                                      <  1:1922936/150‑1 (MQ=255)
                                                                                                               tACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTccc                                      >  1:4974491/1‑150 (MQ=255)
                                                                                                                     tCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATccc                                <  2:3043863/150‑1 (MQ=255)
                                                                                                                       cAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTc                              <  2:4304372/150‑1 (MQ=255)
                                                                                                                           tgcaGGGTACATTGACGCGTTTAAGAGCAAGCGCCTCATACTGCTGTTCTAGACTTAGGTCCTAGTAGTCCTGGCCATCGAGCAGTTGATTACCTGAAGACGTTTGTACATAAGGCAATATTGGCAGGCTCTTTTCCCCATCCCTCGTTa                          <  2:4804947/147‑1 (MQ=255)
                                                                                                                           tGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTa                          >  2:3121511/1‑150 (MQ=255)
                                                                                                                             aaGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTAcc                        <  2:4003685/150‑1 (MQ=255)
                                                                                                                             aaGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTAcc                        >  1:1152828/1‑150 (MQ=255)
                                                                                                                             aaGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTAcc                        >  1:3733065/1‑150 (MQ=255)
                                                                                                                             aaGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTAcc                        <  1:3543689/150‑1 (MQ=255)
                                                                                                                             aaGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATGAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTAcc                        <  2:1881047/150‑1 (MQ=255)
                                                                                                                               gggTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGc                      >  1:120958/1‑150 (MQ=255)
                                                                                                                               gggTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGc                      >  2:4953123/1‑150 (MQ=255)
                                                                                                                                 gTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCaa                    <  2:5034483/150‑1 (MQ=255)
                                                                                                                                 gTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCaa                    <  2:2021115/150‑1 (MQ=255)
                                                                                                                                  tACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAg                   >  1:138241/1‑150 (MQ=255)
                                                                                                                                     aTTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAAt                <  2:1927259/150‑1 (MQ=255)
                                                                                                                                     aTTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAAt                >  1:5046447/1‑150 (MQ=255)
                                                                                                                                     aTTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAAt                <  1:376472/150‑1 (MQ=255)
                                                                                                                                      ttGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAAta               <  2:2424212/150‑2 (MQ=255)
                                                                                                                                       tGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATgg              >  2:1400612/1‑150 (MQ=255)
                                                                                                                                           gcgATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTc          <  1:262130/150‑1 (MQ=255)
                                                                                                                                            cgATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCa         <  2:3044298/150‑1 (MQ=255)
                                                                                                                                             gATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCAt        >  1:337103/1‑150 (MQ=255)
                                                                                                                                             gATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCAt        >  1:5142492/1‑150 (MQ=255)
                                                                                                                                              aTTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTATACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACAGCAAGAATGGTTTCATc       >  1:2301862/1‑150 (MQ=255)
                                                                                                                                              aTTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTTCCCATCCCTCGTTACCGCAAGAATGGTTTCATc       >  2:2657862/1‑150 (MQ=255)
                                                                                                                                              aTTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATc       >  1:2751955/1‑150 (MQ=255)
                                                                                                                                                 aaGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGtt    >  2:1525226/1‑150 (MQ=255)
                                                                                                                                                 aaGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGtt    >  1:3156564/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTAT  >  CP009273/2877950‑2878246

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