Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,152,655 G→T 5.0% R182S (CGT→AGT)  trmA ← tRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,152,6550GT 5.0% 759.6 / 9.9 198R182S (CGT→AGT) trmAtRNA m(5)U54 methyltransferase, SAM‑dependent; tmRNA m(5)U341 methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (96/91);  new base T (8/3);  total (104/94)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 202 total aligned reads displayed.

ACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACG  >  CP009273/4152509‑4152802
                                                                                                                                                  |                                                                                                                                                   
aCGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTGTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACgg                                                                                                                                                    >  1:3040191/1‑148 (MQ=255)
           cTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCAtctc                                                                                                                                         >  2:1763551/1‑148 (MQ=255)
            tAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACTGTAGATCATctct                                                                                                                                        <  2:2984476/148‑1 (MQ=255)
             aaTGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCtt                                                                                                                                       <  2:1893220/148‑1 (MQ=255)
              aTGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCttt                                                                                                                                      <  2:1195404/148‑1 (MQ=255)
              aTGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCttt                                                                                                                                      <  1:1271663/148‑1 (MQ=255)
                gAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTcc                                                                                                                                    >  2:1733362/1‑148 (MQ=255)
                gAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTcc                                                                                                                                    >  1:3315034/1‑148 (MQ=255)
                gAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTcc                                                                                                                                    >  1:3314847/1‑148 (MQ=255)
                 aaaaGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTccc                                                                                                                                   >  1:1079445/1‑148 (MQ=255)
                 aaaaGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTccc                                                                                                                                   <  1:25946/148‑1 (MQ=255)
                 aaaaGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTccc                                                                                                                                   <  2:3505770/148‑1 (MQ=255)
                   aaGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTg                                                                                                                                 >  1:1200225/1‑148 (MQ=255)
                        aCCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGAcc                                                                                                                            >  1:1840578/1‑148 (MQ=255)
                         ccGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCg                                                                                                                           <  1:1451381/148‑1 (MQ=255)
                         ccGGTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCg                                                                                                                           <  2:2576591/148‑1 (MQ=255)
                            ttGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTTCGACCGGCa                                                                                                                        >  2:150767/1‑148 (MQ=255)
                               ccGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCGCTGCGACCGGCAGAc                                                                                                                     >  2:2550837/1‑148 (MQ=255)
                                 gCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGt                                                                                                                   <  1:3127973/148‑1 (MQ=255)
                                 gCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGt                                                                                                                   <  2:972438/148‑1 (MQ=255)
                                  cAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGtt                                                                                                                  >  1:2025725/1‑148 (MQ=255)
                                  cAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGtt                                                                                                                  >  2:2666580/1‑148 (MQ=255)
                                     tACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCAt                                                                                                               <  2:1989796/148‑1 (MQ=255)
                                     tACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCAt                                                                                                               >  2:3025855/1‑148 (MQ=255)
                                     tACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCAt                                                                                                               >  2:1251880/1‑148 (MQ=255)
                                       cAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCg                                                                                                             <  1:26985/148‑1 (MQ=255)
                                       cAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCg                                                                                                             <  1:2853401/148‑1 (MQ=255)
                                        aGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGa                                                                                                            <  1:2888696/148‑1 (MQ=255)
                                        aGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGa                                                                                                            <  1:1886091/148‑1 (MQ=255)
                                        aGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGa                                                                                                            <  2:2347431/148‑1 (MQ=255)
                                        aGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGa                                                                                                            <  1:1159682/148‑1 (MQ=255)
                                        aGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGa                                                                                                            <  1:1267282/148‑1 (MQ=255)
                                        aGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGa                                                                                                            <  1:139315/148‑1 (MQ=255)
                                         gCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGAt                                                                                                           <  2:731249/148‑1 (MQ=255)
                                         gCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGAt                                                                                                           <  2:2919249/148‑1 (MQ=255)
                                          cTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATg                                                                                                          >  1:1671533/1‑148 (MQ=255)
                                           tCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGt                                                                                                         >  1:1130311/1‑148 (MQ=255)
                                            ccAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTa                                                                                                        <  1:246843/148‑1 (MQ=255)
                                            ccAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTa                                                                                                        <  1:172546/148‑1 (MQ=255)
                                             cAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTaa                                                                                                       >  2:3076456/1‑148 (MQ=255)
                                               gTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATc                                                                                                     <  2:665112/148‑1 (MQ=255)
                                                tAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATcc                                                                                                    >  2:3452538/1‑148 (MQ=255)
                                                   aTCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGa                                                                                                 >  2:115092/1‑148 (MQ=255)
                                                     cGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATc                                                                                               <  1:2197522/148‑1 (MQ=255)
                                                     cGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATc                                                                                               <  1:440479/148‑1 (MQ=255)
                                                        cTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAg                                                                                            >  1:1311489/1‑148 (MQ=255)
                                                         ttttGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGc                                                                                           <  1:729115/148‑1 (MQ=255)
                                                         ttttGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGc                                                                                           <  1:2803445/148‑1 (MQ=255)
                                                          tttGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCt                                                                                          >  2:1875054/1‑148 (MQ=255)
                                                             gAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGa                                                                                       >  1:2280526/1‑148 (MQ=255)
                                                             gAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGa                                                                                       >  2:1478073/1‑148 (MQ=255)
                                                              acccTTTGTTAACGACCAGAGCCCATACCAGCATCTGAATATTCATAGCAGCGTTCGGCTGGGTAAAGATGATTTATACCTGACGGTAGATCATCTCTTTACCTCAGACCGGCAGACGATCATCGATGAAATAAAGAACCAGCTCTAt                                                                                      >  2:226313/3‑148 (MQ=255)
                                                              aGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGAt                                                                                      >  2:226327/1‑148 (MQ=255)
                                                                ccTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGAttt                                                                                    >  2:3374728/1‑148 (MQ=255)
                                                                ccTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGAttt                                                                                    <  1:1728184/148‑1 (MQ=255)
                                                                 cTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGAtttt                                                                                   >  2:3473937/1‑148 (MQ=255)
                                                                  tttGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACTGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTg                                                                                  >  2:2377727/1‑148 (MQ=255)
                                                                  tttGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTg                                                                                  >  2:2246592/1‑148 (MQ=255)
                                                                     ggTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGtt                                                                               <  1:1997088/148‑1 (MQ=255)
                                                                            tCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTg                                                                        <  2:2645898/148‑1 (MQ=255)
                                                                                gcgcCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCg                                                                    >  1:1981710/1‑148 (MQ=255)
                                                                                gcgcCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCg                                                                    >  1:174617/1‑148 (MQ=255)
                                                                                gcgcCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCg                                                                    >  1:1116450/1‑148 (MQ=255)
                                                                                   cccATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGAcc                                                                 >  2:2628679/1‑148 (MQ=255)
                                                                                      aTTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAAt                                                              <  1:526605/148‑1 (MQ=255)
                                                                                       ttCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATc                                                             >  1:2631689/1‑148 (MQ=255)
                                                                                          cAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGa                                                          <  1:2961539/148‑1 (MQ=255)
                                                                                            gCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATg                                                        >  1:153542/1‑148 (MQ=255)
                                                                                            gCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATg                                                        >  2:1156623/1‑148 (MQ=255)
                                                                                               tCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGcac                                                     >  2:3553004/1‑148 (MQ=255)
                                                                                               tCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGcac                                                     >  2:2740677/1‑148 (MQ=255)
                                                                                                   aaTATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTc                                                 <  2:2025725/148‑1 (MQ=255)
                                                                                                   aaTATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTc                                                 >  1:2557165/1‑148 (MQ=255)
                                                                                                    atatTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCa                                                <  1:1472695/148‑1 (MQ=255)
                                                                                                      atTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGa                                              <  1:1017095/148‑1 (MQ=255)
                                                                                                      atTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGa                                              >  2:2553786/1‑148 (MQ=255)
                                                                                                      atTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGa                                              <  1:3025039/148‑1 (MQ=255)
                                                                                                        tCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGAtt                                            <  2:594548/148‑1 (MQ=255)
                                                                                                             gccgcGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTgcg                                       <  1:3452538/148‑1 (MQ=255)
                                                                                                              ccgcGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTgcgc                                      >  1:572225/1‑148 (MQ=255)
                                                                                                               cgcgTTCGGCTGGGTAAAGCTGTTTTCTACCTGACTGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTAAGATTCTgcgcg                                     >  2:3523337/1‑148 (MQ=255)
                                                                                                                   ttCGGCTGGGTAAAGCTGTTTTCTACCTGACTGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCa                                 <  2:3397697/148‑1 (MQ=255)
                                                                                                                   ttCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCa                                 >  2:1567677/1‑148 (MQ=255)
                                                                                                                    tCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAg                                <  1:1547302/148‑1 (MQ=255)
                                                                                                                      ggCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTg                              >  2:805010/1‑148 (MQ=255)
                                                                                                                             tAAAGCTGTTTTCTACCTGACTGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATcgcg                       >  1:3155696/1‑148 (MQ=255)
                                                                                                                                aGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAg                    >  1:385457/1‑148 (MQ=255)
                                                                                                                                aGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAg                    >  2:3353946/1‑148 (MQ=255)
                                                                                                                                aGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAg                    <  2:1338982/148‑1 (MQ=255)
                                                                                                                                   tGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGc                 >  2:2182410/1‑148 (MQ=255)
                                                                                                                                     ttttCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCt               >  1:1639126/1‑148 (MQ=255)
                                                                                                                                      tttCTACCTGACGGTAGATCATCTCTTCCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACAGTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTc              >  2:227487/1‑148 (MQ=255)
                                                                                                                                        tCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCg            >  1:180692/1‑148 (MQ=255)
                                                                                                                                         cTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGc           <  2:517966/148‑1 (MQ=255)
                                                                                                                                            ccTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTc        >  1:1043759/1‑148 (MQ=255)
                                                                                                                                             cTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTcc       <  1:958546/148‑1 (MQ=255)
                                                                                                                                             cTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTcc       <  2:1936664/148‑1 (MQ=255)
                                                                                                                                               gACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTg     >  2:2587951/1‑148 (MQ=255)
                                                                                                                                                aCGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGa    >  1:2145429/1‑148 (MQ=255)
                                                                                                                                                  ggTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACg  <  1:1112591/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACGCGCCAGCGCTAATGAAAAGTTACCGTTGCCGCAGTACAGCTCCAGTAAATCGCCTTTTGAGCCTTTGGTTACGTCCAGCGCCCATTCCAGCATCTGAATATTCATCGCCGCGTTCGGCTGGGTAAAGCTGTTTTCTACCTGACGGTAGATCATCTCTTTCCCTGCGACCGGCAGACGTTCATCGATGTAATCCTGATCCAGCTCGATTTTGGTTTTCGTTGCCCGACCAATCAGATGCACATTCAGATTCTGCGCGCGCAGTGCATCGCGCAGGGCCTCCGCTTCCTGACG  >  CP009273/4152509‑4152802

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