Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,884,978 C→A 5.3% D183Y (GAT→TAT)  tsaB ← tRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; binding partner and protease for TsaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,884,9780CA 5.3% 795.5 / 13.0 206D183Y (GAT→TAT) tsaBtRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; binding partner and protease for TsaD
Reads supporting (aligned to +/- strand):  ref base C (115/80);  new base A (5/6);  total (120/86)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 207 total aligned reads displayed.

TTCTTTGCCCGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGA  >  CP009273/1884832‑1885124
                                                                                                                                                  |                                                                                                                                                  
   tttGCCCGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCgcgc                                                                                                                                                >  1:2555372/1‑148 (MQ=255)
    ttGCCCGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCa                                                                                                                                               >  2:761670/1‑148 (MQ=255)
         cGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAAcc                                                                                                                                          >  2:1687945/1‑148 (MQ=255)
         cGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAAcc                                                                                                                                          >  1:622885/1‑148 (MQ=255)
          ggAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCa                                                                                                                                         <  1:35535/148‑1 (MQ=255)
          ggAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATAGCGCAAAACCa                                                                                                                                         <  1:830120/148‑1 (MQ=255)
           gAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAg                                                                                                                                        <  1:2191177/148‑1 (MQ=255)
           gAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAg                                                                                                                                        <  1:2188377/148‑1 (MQ=255)
           gAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAg                                                                                                                                        >  1:2159517/1‑148 (MQ=255)
            aaGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGc                                                                                                                                       >  2:18984/1‑148 (MQ=255)
            aaGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGc                                                                                                                                       >  2:2486820/1‑148 (MQ=255)
            aaGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGCTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGc                                                                                                                                       >  2:1539888/1‑148 (MQ=255)
            aaGTTTCTTCCATGCGACGTTGTTAAGTAAATAAACCGTTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGc                                                                                                                                       <  2:3504493/148‑1 (MQ=255)
               tttCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCg                                                                                                                                    >  1:1380215/1‑148 (MQ=255)
                ttcttcCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGc                                                                                                                                   >  2:1926333/1‑148 (MQ=255)
                 tcttcCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCAAATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGct                                                                                                                                  >  2:3256257/1‑148 (MQ=255)
                  cttcCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGctc                                                                                                                                 >  2:2710929/1‑148 (MQ=255)
                   ttcCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGctct                                                                                                                                >  1:317436/1‑148 (MQ=255)
                     ccATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGctctct                                                                                                                              <  2:2216713/148‑1 (MQ=255)
                      cATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCtt                                                                                                                             <  2:680061/148‑1 (MQ=255)
                         gCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTAc                                                                                                                          >  1:3012222/1‑148 (MQ=255)
                            aCGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCga                                                                                                                       >  2:2745000/1‑148 (MQ=255)
                                tgttACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATc                                                                                                                   >  2:45655/1‑148 (MQ=255)
                                 gttACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATcc                                                                                                                  <  1:298862/148‑1 (MQ=255)
                                    aCGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGc                                                                                                               >  1:2244580/1‑148 (MQ=255)
                                      gTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCa                                                                                                             <  2:2585421/148‑1 (MQ=255)
                                      gTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATAGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCa                                                                                                             <  2:692044/148‑1 (MQ=255)
                                        aaataaaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCTAGATCCGGCCAAg                                                                                                           >  2:1434523/1‑148 (MQ=255)
                                         aataaaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGc                                                                                                          >  1:723710/1‑148 (MQ=255)
                                          ataaaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCt                                                                                                         <  2:1328935/148‑1 (MQ=255)
                                             aaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGc                                                                                                      <  2:1603594/148‑1 (MQ=255)
                                             aaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGc                                                                                                      >  2:1724581/1‑148 (MQ=255)
                                              aCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGcc                                                                                                     <  2:473567/148‑1 (MQ=255)
                                                  gTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGcc                                                                                                 >  1:2745116/1‑148 (MQ=255)
                                                    tCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCg                                                                                               >  1:3313098/1‑148 (MQ=255)
                                                       gCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTcc                                                                                            <  1:2856694/148‑1 (MQ=255)
                                                       gCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTcc                                                                                            <  1:2466444/148‑1 (MQ=255)
                                                       gCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTcc                                                                                            <  1:390366/148‑1 (MQ=255)
                                                         aTGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCt                                                                                          >  1:1044310/1‑148 (MQ=255)
                                                             tCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCg                                                                                      <  1:784351/148‑1 (MQ=255)
                                                               caccgccaccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTc                                                                                    >  1:1039819/1‑148 (MQ=255)
                                                                     caccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCAt                                                                              <  1:906921/148‑1 (MQ=255)
                                                                       ccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTc                                                                            >  2:859508/1‑148 (MQ=255)
                                                                       ccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTc                                                                            >  2:2764014/1‑148 (MQ=255)
                                                                       ccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCCTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTc                                                                            >  2:625528/1‑148 (MQ=255)
                                                                           tttACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCg                                                                        >  1:1554791/1‑148 (MQ=255)
                                                                           tttAACCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGAAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCg                                                                        >  1:3655250/1‑148 (MQ=255)
                                                                            ttACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCgg                                                                       >  2:3220112/1‑148 (MQ=255)
                                                                            ttACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCgg                                                                       >  2:3220288/1‑148 (MQ=255)
                                                                                ccTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAg                                                                   <  1:1962765/148‑1 (MQ=255)
                                                                                  tCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGACAGCt                                                                 >  2:1495984/1‑148 (MQ=255)
                                                                                   cAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTg                                                                >  2:2322816/1‑148 (MQ=255)
                                                                                   cAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTg                                                                <  1:18984/148‑1 (MQ=255)
                                                                                   cAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTg                                                                >  2:397152/1‑148 (MQ=255)
                                                                                    aGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGt                                                               <  1:2927159/148‑1 (MQ=255)
                                                                                    aGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGt                                                               <  2:1105123/148‑1 (MQ=255)
                                                                                      cAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTg                                                             >  1:1163637/1‑148 (MQ=255)
                                                                                            tCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCg                                                       >  2:3690544/1‑148 (MQ=255)
                                                                                                acacGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCa                                                   >  2:3436779/1‑148 (MQ=255)
                                                                                                 cacGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCAt                                                  <  2:3777323/148‑1 (MQ=255)
                                                                                                      aaTCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATAGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACg                                             <  2:2922218/148‑1 (MQ=255)
                                                                                                        tCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGAt                                           <  2:3067359/148‑1 (MQ=255)
                                                                                                         cGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGAtt                                          >  2:3737995/1‑148 (MQ=255)
                                                                                                          ggCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGAttt                                         <  2:1378423/148‑1 (MQ=255)
                                                                                                          ggCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGAttt                                         <  2:2526757/148‑1 (MQ=255)
                                                                                                            cagcaTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCg                                       <  1:1864938/148‑1 (MQ=255)
                                                                                                            cagcaTATCTTCAGCAGCAGGCAGTAACACTTCGCCATAGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCg                                       >  2:78496/1‑148 (MQ=255)
                                                                                                            cagcaTATCTTCAGCAGCAGGCAGTAACACTTCGCCATAGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATGTCg                                       >  2:79242/1‑148 (MQ=255)
                                                                                                              gcaTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCggg                                     <  2:1930733/148‑1 (MQ=255)
                                                                                                               caTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGCTCCGGCCNAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGt                                    <  1:1938330/148‑1 (MQ=255)
                                                                                                               caTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGt                                    >  1:3710905/1‑148 (MQ=255)
                                                                                                               caTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGt                                    >  1:26376/1‑148 (MQ=255)
                                                                                                                 tatCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAACGCTGTTGCATTCGTTCATGGACGATTTCGGGttt                                  >  2:2618860/1‑148 (MQ=255)
                                                                                                                 tatCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGttt                                  >  2:3396158/1‑148 (MQ=255)
                                                                                                                  atCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACAAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTg                                 >  1:663211/1‑148 (MQ=255)
                                                                                                                    cTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAg                               >  1:1036247/1‑148 (MQ=255)
                                                                                                                      tCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTa                             >  2:2270678/1‑148 (MQ=255)
                                                                                                                      tCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTa                             >  2:1999858/1‑148 (MQ=255)
                                                                                                                        agcagcagGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACg                           <  2:9508/148‑1 (MQ=255)
                                                                                                                         gcagcagGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACgg                          <  1:2463125/148‑1 (MQ=255)
                                                                                                                          cagcagGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGc                         <  1:2745000/148‑1 (MQ=255)
                                                                                                                            gcaggcagTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCtt                       <  2:1240416/148‑1 (MQ=255)
                                                                                                                              aggcagTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCg                     <  1:2900562/148‑1 (MQ=255)
                                                                                                                               ggcagTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCgg                    >  1:2112950/1‑148 (MQ=255)
                                                                                                                                gcagTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGt                   >  1:1357323/1‑148 (MQ=255)
                                                                                                                                gcagTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGt                   <  1:1518827/148‑1 (MQ=255)
                                                                                                                                 cagTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGtt                  >  2:3115501/1‑148 (MQ=255)
                                                                                                                                  agTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGttt                 <  1:1444635/148‑1 (MQ=255)
                                                                                                                                   gTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttc                >  1:33316/1‑148 (MQ=255)
                                                                                                                                    tAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttct               >  2:1751853/1‑148 (MQ=255)
                                                                                                                                    tAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttct               >  2:2448696/1‑148 (MQ=255)
                                                                                                                                     aaCACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttctt              >  1:1563004/1‑148 (MQ=255)
                                                                                                                                      acacTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttcttc             <  2:2159517/148‑1 (MQ=255)
                                                                                                                                       cacTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCa            <  1:1667884/148‑1 (MQ=255)
                                                                                                                                         cTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCAcc          >  2:1408749/1‑148 (MQ=255)
                                                                                                                                           tCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGt        >  2:1250559/1‑148 (MQ=255)
                                                                                                                                              ccATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGcc     <  2:26376/148‑1 (MQ=255)
                                                                                                                                               cATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCa    >  2:1377935/1‑148 (MQ=255)
                                                                                                                                               cATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCa    >  2:2426673/1‑148 (MQ=255)
                                                                                                                                                 tCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGa  <  2:2184588/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTCTTTGCCCGGAAGTTTCTTCCATGCGACGTTGTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGA  >  CP009273/1884832‑1885124

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: