Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,426,592 T→A 17.5% intergenic (+35/+36) rplI → / ← yjfZ 50S ribosomal subunit protein L9/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,426,5920TA17.5% 376.3 / 70.1 149intergenic (+35/+36)rplI/yjfZ50S ribosomal subunit protein L9/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (64/59);  new base A (10/16);  total (74/75)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 149 total aligned reads displayed.

TGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGG  >  NC_000913/4426496‑4426690
                                                                                                |                                                                                                  
tGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTgg                                                                                                 <  2:1063759/100‑1 (MQ=255)
 gAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGt                                                                                                >  2:3525057/1‑100 (MQ=255)
  aGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTc                                                                                               <  2:2422668/100‑1 (MQ=255)
    cTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCgg                                                                                             >  2:2899753/1‑100 (MQ=255)
     ttCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGc                                                                                            >  2:2780532/1‑100 (MQ=255)
      tCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCg                                                                                           <  2:1455601/100‑1 (MQ=255)
      tCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCg                                                                                           >  1:617027/1‑100 (MQ=255)
        cAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGtt                                                                                         >  1:3367149/1‑100 (MQ=255)
        cAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGtt                                                                                         >  2:413473/1‑100 (MQ=255)
        cAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtt                                                                                         <  1:936798/100‑1 (MQ=255)
        cAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtt                                                                                         <  1:2899733/100‑1 (MQ=255)
        cAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtt                                                                                         <  2:2915739/100‑1 (MQ=255)
         aGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGttt                                                                                        <  1:3492601/100‑1 (MQ=255)
         aGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGttt                                                                                        <  1:951887/100‑1 (MQ=255)
          ggTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtttt                                                                                       <  1:3215045/100‑1 (MQ=255)
          ggTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGttt                                                                                        <  2:3403154/99‑1 (MQ=255)
          ggTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAAAGCCGACCATTGCTCGGCGtttt                                                                                       >  1:2728099/1‑100 (MQ=255)
           gTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTg                                                                                      <  2:424097/100‑1 (MQ=255)
           gTTCACAGCGAAGTATTCGCGAAAGTGAACGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGttttt                                                                                      <  1:3379187/100‑2 (MQ=255)
             tCACAGCGAAGTATTCGCGAAAGTGGTCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtttttac                                                                                    <  1:611973/100‑4 (MQ=255)
              cacaGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCtt                                                                                   >  1:3147247/1‑100 (MQ=255)
              cacaGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCtt                                                                                   <  2:2392337/100‑1 (MQ=255)
              cacaGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATGGGTCGGCGTTTTGCtt                                                                                   >  2:1050295/1‑100 (MQ=255)
              cacaGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtttttacg                                                                                   <  2:1617854/100‑5 (MQ=255)
              cacaGCGAAGTATTCGCGAAAGCGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtttttacg                                                                                   <  2:1881869/100‑5 (MQ=38)
              atctgCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtttttacg                                                                                   <  2:699985/96‑5 (MQ=38)
               acaGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCttt                                                                                  <  2:3450884/100‑1 (MQ=255)
                caGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTc                                                                                 >  1:2239212/1‑100 (MQ=255)
                caGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTc                                                                                 <  2:792817/100‑1 (MQ=255)
                caGCGAAGTATTCGCGAAAGTGAGCGTAAACGGAGAAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtttttacg                                                                                   <  1:2837449/98‑5 (MQ=37)
                  gCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTa                                                                               >  2:1225434/1‑100 (MQ=255)
                  gCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTa                                                                               >  1:1088191/1‑100 (MQ=255)
                  gCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGCTTTGCTTTCTa                                                                               >  2:447855/1‑100 (MQ=255)
                  gCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCAATGGTCGGCGtttttacgagat                                                                               >  1:699981/1‑92 (MQ=38)
                      aGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTAtttt                                                                           <  1:3032618/100‑1 (MQ=255)
                      aGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTAtttt                                                                           <  2:950374/100‑1 (MQ=255)
                        tATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTc                                                                         >  2:2756013/1‑100 (MQ=255)
                         aTTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCg                                                                        <  2:1027819/100‑1 (MQ=255)
                         aTTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCg                                                                        <  1:1448822/100‑1 (MQ=255)
                         aTTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTCCTTTCTATTTTTCg                                                                        >  1:2421634/1‑100 (MQ=255)
                           tCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTc                                                                      <  2:194013/100‑1 (MQ=255)
                            cgcgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCa                                                                     >  1:569382/1‑100 (MQ=255)
                             gcgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAg                                                                    <  2:326454/100‑1 (MQ=255)
                             gcgAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAg                                                                    >  2:208824/1‑100 (MQ=255)
                               gAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGt                                                                  >  2:1653401/1‑100 (MQ=255)
                                aaaGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTa                                                                 <  2:2515295/100‑1 (MQ=255)
                                 aaGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTAt                                                                >  1:2659577/1‑100 (MQ=255)
                                 aaGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTAt                                                                >  1:3107596/1‑100 (MQ=255)
                                  aGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTTGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTAtt                                                               >  1:998887/1‑100 (MQ=255)
                                  aGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTAtt                                                               <  1:2825581/100‑1 (MQ=255)
                                   gTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTa                                                              >  2:2029994/1‑100 (MQ=255)
                                   gTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTa                                                              >  2:265977/1‑100 (MQ=255)
                                    tGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAg                                                             >  1:3427658/1‑100 (MQ=255)
                                     gATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGt                                                            >  1:3179550/1‑100 (MQ=255)
                                       tCGTAAACGTAGTAGCTGAATAATTCGTTATTCAAGGAGGCGTAAAAACGCCGACCATTGGGCGGAGTTTTGCGTTCTATTTTTCGTGAGGTATTAGttt                                                          >  1:2832184/1‑100 (MQ=255)
                                       tCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGttt                                                          >  2:876377/1‑100 (MQ=255)
                                        cGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTc                                                         >  1:829925/1‑100 (MQ=255)
                                        cGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTc                                                         <  1:2645880/100‑1 (MQ=255)
                                              cGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGt                                                   >  2:1553506/1‑100 (MQ=255)
                                               gtagtaGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTa                                                  <  2:3323349/100‑1 (MQ=255)
                                               gtagtaGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTCCTATTTTTCGTCAGGTATTAGTTTCGCAAGTa                                                  <  1:86831/100‑1 (MQ=255)
                                                   tagCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATc                                              <  2:3147269/100‑1 (MQ=255)
                                                      cTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCAt                                           >  1:2758574/1‑100 (MQ=255)
                                                       tGAGTAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGTCGCTTTGCTTTATATTTTTCGTCACGTACTAGTTTGGCAAGTTCATCCAt                                           >  1:3060374/1‑99 (MQ=255)
                                                        gAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTc                                         >  1:2094085/1‑100 (MQ=255)
                                                          ataatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCa                                        <  1:554865/99‑1 (MQ=255)
                                                          ataatTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCAc                                       <  1:1177774/100‑1 (MQ=255)
                                                             atTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCAt                                    <  2:617030/100‑1 (MQ=255)
                                                             atTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCAt                                    >  2:2531855/1‑100 (MQ=255)
                                                              ttCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCAtt                                   >  1:1730541/1‑100 (MQ=255)
                                                                cGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATtct                                 <  1:1842775/100‑1 (MQ=255)
                                                                  ttATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTctct                               <  1:767148/100‑1 (MQ=255)
                                                                    aTTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtg                             >  2:3307692/1‑100 (MQ=255)
                                                                     tgcAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGTTCCATTCACCATTCTCTtgt                            <  1:2163638/98‑1 (MQ=255)
                                                                      tCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtgtg                           <  1:1087646/100‑1 (MQ=255)
                                                                      tCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTtgtg                           <  2:2039437/100‑1 (MQ=255)
                                                                       cAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTgtgt                          <  2:886368/100‑1 (MQ=255)
                                                                       cAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTgtgt                          >  1:2538151/1‑100 (MQ=255)
                                                                        aaCGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTCGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTc                         >  1:1898539/1‑100 (MQ=255)
                                                                             gaaagcaaaaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATaa                    >  2:2561640/7‑100 (MQ=255)
                                                                               cGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGt                  >  2:1871453/1‑100 (MQ=255)
                                                                               cGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGt                  >  2:2087243/1‑100 (MQ=255)
                                                                                  aaaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAg               >  1:1024465/1‑100 (MQ=255)
                                                                                  aaaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAg               <  1:1321421/100‑1 (MQ=255)
                                                                                  aaaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAg               <  1:2599257/100‑1 (MQ=255)
                                                                                   aaaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAgg              >  2:2175750/1‑100 (MQ=255)
                                                                                   aaaaCGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAgg              >  1:2253542/1‑100 (MQ=255)
                                                                                    aaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa             >  1:1390873/1‑100 (MQ=255)
                                                                                    aaaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGa             >  2:1082877/1‑100 (MQ=255)
                                                                                     aaCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGAt            <  2:2418052/100‑1 (MQ=255)
                                                                                      aCGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATg           >  2:1907484/1‑100 (MQ=255)
                                                                                       cGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTTTTTTTTTATTTTTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATgg          >  1:1605009/1‑100 (MQ=255)
                                                                                       cGCCGACCAATGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATgg          >  2:2447866/1‑100 (MQ=255)
                                                                                         ccGACCATTGGTCGGCGTTTTGTTTTTTTTTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGtt        >  1:1851697/1‑100 (MQ=255)
                                                                                          cGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGttt       >  2:1307293/1‑100 (MQ=255)
                                                                                            aCCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATCAGTCAGGATGGTTTGc     <  2:825971/100‑1 (MQ=255)
                                                                                             ccATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTATTGTGTCCATAAGTCAGGATGGTTTGc     >  1:2105590/1‑99 (MQ=255)
                                                                                              cATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTg   <  1:2448403/100‑1 (MQ=255)
                                                                                              cATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTg   >  2:2127028/1‑100 (MQ=255)
                                                                                              cATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTg   >  2:858996/1‑100 (MQ=255)
                                                                                                |                                                                                                  
TGAGCTTCCAGGTTCACAGCGAAGTATTCGCGAAAGTGATCGTAAACGTAGTAGCTGAATAATTCGTTATTCAACGAGACGTAAAAACGCCGACCATTGGTCGGCGTTTTGCTTTCTATTTTTCGTCAGGTATTAGTTTCGCAAGTAGATCCATTCACCATTCTCTTGTGTCCATAAGTCAGGATGGTTTGCTGG  >  NC_000913/4426496‑4426690

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