Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,878,098 T→A 45.7% 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,0980TA45.7% 78.4 / inf 225A827A (GCA→GCT) ygcBR‑loop helicase‑annealase Cas3 needed for Cascade anti‑viral activity
Reads supporting (aligned to +/- strand):  ref base T (56/66);  new base A (46/57);  total (102/123)
Fisher's exact test for biased strand distribution p-value = 8.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 228 total aligned reads displayed.

ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTA  >  CP009273/2877951‑2878245
                                                                                                                                                   |                                                                                                                                                   
  ccccTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGc                                                                                                                                                   >  1:1770759/1‑148 (MQ=255)
    ccTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCaa                                                                                                                                                 >  1:308076/1‑148 (MQ=255)
    ccTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCaa                                                                                                                                                 >  1:3369189/1‑148 (MQ=255)
     cTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAg                                                                                                                                                <  2:942909/148‑1 (MQ=255)
           aGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTc                                                                                                                                          >  1:47753/1‑148 (MQ=255)
           aGGTAATAAAAATACTGTTACCCTGCCAGACCCATCCATCCATATTCTGTTTCCCTTCAAGCCAAAGTAACCCATCTTCATATACTACTTCAGAAAAACTACTTTTCAAGGTGAAGGTTACTTTGACGCGATTAAGTGCAAGCGCCTc                                                                                                                                          <  1:1974581/148‑1 (MQ=255)
              taataaCAATACTGTTACCCTGCCAGACCCATCAATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTAATTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATa                                                                                                                                       <  1:2213902/148‑1 (MQ=255)
                ataaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATAct                                                                                                                                     >  2:2148281/1‑148 (MQ=255)
                 taaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATActg                                                                                                                                    >  2:274640/1‑148 (MQ=255)
                  aaCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATActgc                                                                                                                                   <  1:3153042/148‑1 (MQ=255)
                   aaaaTACTGTTACCCTGCCATACCCATCCATCCATTTTCTGTTTCCCTTCAAGACAAAGTAACCCATCTTCATCTACTAATTCAGAAAAACTACTTTTCCTGGTGAAGGTTACATTTACGAGATTAAGAGCAAGCGCCTCATActgct                                                                                                                                  <  1:728357/146‑1 (MQ=255)
                   aCAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATActgct                                                                                                                                  <  2:2339018/148‑1 (MQ=255)
                    cAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTGTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATActgctg                                                                                                                                 <  1:1924822/148‑1 (MQ=255)
                     aaTACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGt                                                                                                                                <  2:982496/148‑1 (MQ=255)
                      aTACTTTTACCCTGCCACACCCATCCATCCAGTTTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTAATAATACATAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGtt                                                                                                                               <  1:1516728/148‑1 (MQ=255)
                       tACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCACGCGCCTCATACTGCTGTTc                                                                                                                              >  1:2855148/1‑148 (MQ=255)
                       tACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTc                                                                                                                              <  2:1141736/148‑1 (MQ=255)
                            ttACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGAc                                                                                                                         >  1:434203/1‑148 (MQ=255)
                            ttACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGAc                                                                                                                         <  2:3112391/148‑1 (MQ=255)
                             tACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACt                                                                                                                        >  1:2919512/1‑148 (MQ=255)
                              aCCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACtt                                                                                                                       >  2:2588598/1‑148 (MQ=255)
                              aCCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCATCCAGCCAAAGTAACCCATCTACATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCACGCGCCTCATACTGCTGTTCATGACtt                                                                                                                       >  2:1536632/1‑148 (MQ=255)
                                   gCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTc                                                                                                                  <  1:831829/148‑1 (MQ=255)
                                    ccAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTcc                                                                                                                 <  1:1445410/148‑1 (MQ=255)
                                    ccAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCACGCGCCTCATACTGCTGTTCATGACTTAGGTcc                                                                                                                 >  2:397585/1‑148 (MQ=255)
                                    ccAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTcc                                                                                                                 <  2:2341794/148‑1 (MQ=255)
                                     cAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCt                                                                                                                <  1:1289448/148‑1 (MQ=255)
                                      aGACCCATCCATCCAGAGTCTGTGTCCCTTCCATCCAAAGTAACACATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTc                                                                                                               <  1:3040581/148‑1 (MQ=255)
                                       gACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCg                                                                                                              >  1:506985/1‑148 (MQ=255)
                                       gACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCg                                                                                                              <  1:2290700/148‑1 (MQ=255)
                                         cccATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTa                                                                                                            >  1:2375797/1‑148 (MQ=255)
                                            atccatccAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGAc                                                                                                         <  2:3217297/148‑1 (MQ=255)
                                             tccatccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGAcc                                                                                                        <  1:289086/148‑1 (MQ=255)
                                               catccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTg                                                                                                      >  2:427972/1‑148 (MQ=255)
                                                atccaGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTgg                                                                                                     <  2:902402/148‑1 (MQ=255)
                                                   caGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCa                                                                                                  <  1:1380816/148‑1 (MQ=255)
                                                     gATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATc                                                                                                >  2:2772820/1‑148 (MQ=255)
                                                     gATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATc                                                                                                >  2:498225/1‑148 (MQ=255)
                                                      aTTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTTAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCg                                                                                               <  1:3122097/148‑1 (MQ=255)
                                                       ttCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGa                                                                                              >  2:2336823/1‑148 (MQ=255)
                                                           gTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAg                                                                                          <  2:2534655/148‑1 (MQ=255)
                                                            tttCACTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGGACATTGACGCGATTAAGAGCAAGCGCCGCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCCTCGAGCAgt                                                                                         >  2:1897294/1‑148 (MQ=255)
                                                                 cTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGttt                                                                                    <  1:1892278/148‑1 (MQ=255)
                                                                       gCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGa                                                                              <  1:2832536/148‑1 (MQ=255)
                                                                       gCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGa                                                                              >  2:1419088/1‑148 (MQ=255)
                                                                         cAATGTAACCCTTCTTCATCTAGTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGTCGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCTTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAg                                                                            <  1:642262/148‑1 (MQ=255)
                                                                            aGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACg                                                                         <  1:2148281/148‑1 (MQ=255)
                                                                            aGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACg                                                                         >  1:3165280/1‑148 (MQ=255)
                                                                             gTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGt                                                                        >  1:1068039/1‑148 (MQ=255)
                                                                               aaCCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGttt                                                                      >  2:1994756/1‑148 (MQ=255)
                                                                                 cccATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGt                                                                    <  2:414892/148‑1 (MQ=255)
                                                                                  ccATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTa                                                                   <  1:2332984/148‑1 (MQ=255)
                                                                                    aTCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACa                                                                 >  2:2397261/1‑148 (MQ=255)
                                                                                    aTCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACa                                                                 >  2:564528/1‑148 (MQ=255)
                                                                                     tCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACAt                                                                <  1:2598547/148‑1 (MQ=255)
                                                                                      cTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATa                                                               >  1:1619318/1‑148 (MQ=255)
                                                                                        tCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAg                                                             <  1:1833628/148‑1 (MQ=255)
                                                                                        tCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAg                                                             >  2:335778/1‑148 (MQ=255)
                                                                                        tCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGGTTACCTGAAGACGTTTGTACATAAg                                                             >  2:1850589/1‑148 (MQ=255)
                                                                                         cATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAgg                                                            <  2:2358066/148‑1 (MQ=255)
                                                                                         cATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAgg                                                            <  2:2228239/148‑1 (MQ=255)
                                                                                           tCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCa                                                          <  2:536323/148‑1 (MQ=255)
                                                                                            ctactaCTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaa                                                         <  2:2938740/148‑1 (MQ=255)
                                                                                             tactacTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaat                                                        <  1:765913/148‑1 (MQ=255)
                                                                                              actactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaata                                                       <  2:2845062/148‑1 (MQ=255)
                                                                                               ctactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaataa                                                      <  2:1205381/148‑1 (MQ=255)
                                                                                               ctactTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCaataa                                                      <  1:272441/148‑1 (MQ=255)
                                                                                                      aGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAgg                                               <  2:1724633/148‑1 (MQ=255)
                                                                                                       gAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGc                                              <  2:1705774/148‑1 (MQ=255)
                                                                                                        aaaaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCt                                             >  2:245859/1‑148 (MQ=255)
                                                                                                         aaaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTc                                            <  1:2572561/148‑1 (MQ=255)
                                                                                                          aaaCTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCa                                           >  1:802697/1‑148 (MQ=255)
                                                                                                             cTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCAttt                                        >  1:3491632/1‑148 (MQ=255)
                                                                                                             cTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCAttt                                        >  1:1712334/1‑148 (MQ=255)
                                                                                                              tACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGAAATAATGGCAGGCTCATTTc                                       >  2:2689847/1‑148 (MQ=255)
                                                                                                               aCGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTcc                                      <  2:1612425/148‑1 (MQ=255)
                                                                                                               aCGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTcc                                      >  2:1175976/1‑148 (MQ=255)
                                                                                                                cGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTccc                                     >  2:1268750/1‑148 (MQ=255)
                                                                                                                  tttCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCa                                   >  2:2852623/1‑148 (MQ=255)
                                                                                                                    tCCAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGGGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATc                                 <  2:2735377/148‑1 (MQ=255)
                                                                                                                     ccAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATcc                                <  2:3407073/148‑1 (MQ=255)
                                                                                                                      cAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATccc                               <  2:21223/148‑1 (MQ=255)
                                                                                                                      cAGGTGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATccc                               >  1:90660/1‑148 (MQ=255)
                                                                                                                        ggTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTc                             >  1:2323214/1‑148 (MQ=255)
                                                                                                                          tGAAGGGTACATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGt                           >  2:1841851/1‑148 (MQ=255)
                                                                                                                                 tACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCa                    <  2:1668677/148‑1 (MQ=255)
                                                                                                                                  aCATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCaa                   >  1:1756589/1‑148 (MQ=255)
                                                                                                                                   cATTGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAg                  <  1:2345218/148‑1 (MQ=255)
                                                                                                                                    aTTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGa                 <  2:43715/148‑1 (MQ=255)
                                                                                                                                     ttGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGaa                >  2:3241317/1‑148 (MQ=255)
                                                                                                                                     ttGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGaa                <  1:536773/148‑1 (MQ=255)
                                                                                                                                     ttGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGaa                >  2:2416620/1‑148 (MQ=255)
                                                                                                                                      tGACGCGATTAAGAGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAAt               <  2:3091947/148‑1 (MQ=255)
                                                                                                                                        aCGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATgg             <  1:997099/148‑1 (MQ=255)
                                                                                                                                         cgcgATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGt            <  1:3463363/148‑1 (MQ=255)
                                                                                                                                         cgcgATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGt            <  1:1994133/148‑1 (MQ=255)
                                                                                                                                          gcgATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGtt           <  1:1237105/148‑1 (MQ=255)
                                                                                                                                             aTTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCa        <  1:806098/148‑1 (MQ=255)
                                                                                                                                             aTTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCa        <  1:3385573/148‑1 (MQ=255)
                                                                                                                                                   tGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTa  <  2:1770788/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGACGCGATTAAGTGCAAGCGCCTCATACTGCTGTTCATGACTTAGGTCCTCGTAGACCTGGCCATCGAGCAGTTGTTTACCTGAAGACGTTTGTACATAAGGCAATAATGGCAGGCTCATTTCCCCATCCCTCGTTACCGCAAGAATGGTTTCATCGTTA  >  CP009273/2877951‑2878245

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑