Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,589 G→C 15.8% A97P (GCA→CCA)  yafO → mRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913252,5890GC15.8% 265.5 / 41.6 101A97P (GCA→CCA) yafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (42/43);  new base C (7/9);  total (49/52)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

Only 100 of 112 total aligned reads displayed.

GCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCG  >  NC_000913/252490‑252686
                                                                                                   |                                                                                                 
gCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                                                                                   <  2:151562/100‑1 (MQ=255)
gCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                                                                                   >  2:1551393/1‑100 (MQ=255)
gCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                                                                                   <  1:2120668/100‑1 (MQ=255)
 cAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaagc                                                                                                  >  1:1014455/1‑100 (MQ=255)
 cAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGagcaag                                                                                                   <  1:1366772/99‑1 (MQ=255)
 cAAATGAGAATAATACATTTAAGAGAAAGTTGAGAAAATTAAGAAGAACGAATGAAGAAGCGTATTTGGTATATTGTGAGGGGGCGTTTGATGagcaag                                                                                                   >  2:1902035/1‑99 (MQ=255)
      gagaATAATCCATTTCCGCCACAGTTGTGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGt                                                                                             <  2:664941/100‑1 (MQ=255)
      gagaATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGt                                                                                             <  1:2796190/100‑1 (MQ=255)
        gaATAATTCATGGCCGCCACATGTGCTTCAAGGCAGCAGAACGAATGACGAATCGCATTTGTTATATGGTCAGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                           <  1:2538140/100‑1 (MQ=255)
          ataatCCATTTCCGCCACAGTTGCGCCAATTCAGCATAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCt                                                                                         <  1:423189/100‑1 (MQ=255)
          ataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCt                                                                                         >  2:658451/1‑100 (MQ=255)
           ttatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTc                                                                                        >  1:2057311/3‑100 (MQ=255)
           taatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTc                                                                                        >  2:472484/1‑100 (MQ=255)
            aatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCa                                                                                       <  2:571014/100‑1 (MQ=255)
            aatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCa                                                                                       <  2:325537/100‑1 (MQ=255)
            aatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGAAGAAGCGCATTTGGTATATTATCAGGGAGCGTTTGATGAGCAAGAATGGTTGCTCa                                                                                       >  2:1156471/1‑100 (MQ=255)
             atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                      >  1:2390851/1‑100 (MQ=255)
             atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAt                                                                                      >  1:2858988/1‑100 (MQ=255)
             atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAt                                                                                      <  2:1976367/100‑1 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAtt                                                                                     >  1:102440/1‑100 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCAtt                                                                                     >  1:102276/1‑100 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                     <  2:2901334/100‑2 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                     <  1:1430103/100‑2 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                     <  1:1066501/100‑2 (MQ=255)
              tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                     <  1:2481304/100‑2 (MQ=255)
               ccATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtca                                                                                    <  1:1122385/100‑3 (MQ=255)
               ccATTTCCGCCAAAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                     <  2:1072258/99‑2 (MQ=255)
                cATTTCCGCCACAGTTGCGTCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGc                                                                                   <  1:2612136/100‑1 (MQ=255)
                  tttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtcaaac                                                                                 <  1:1284252/100‑6 (MQ=39)
                    tCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCAtt                                                                               >  1:2358521/1‑100 (MQ=255)
                     ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTc                                                                              <  2:2999502/100‑1 (MQ=255)
                     ccGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTc                                                                              >  1:147942/1‑100 (MQ=255)
                      cGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                             <  1:2238258/100‑1 (MQ=255)
                       gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTg                                                                            >  1:2377876/1‑100 (MQ=255)
                       gCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGAACGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTg                                                                            >  1:757706/1‑100 (MQ=255)
                         cacaGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGaa                                                                          <  1:1039234/100‑1 (MQ=255)
                           caGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAc                                                                        <  2:876782/100‑1 (MQ=255)
                               tGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGa                                                                    <  2:1167439/100‑1 (MQ=255)
                               tGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGa                                                                    >  2:1630867/1‑100 (MQ=255)
                                 cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGaa                                                                   <  2:256801/99‑1 (MQ=255)
                                 cgcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGaa                                                                   <  1:554372/99‑1 (MQ=255)
                                  gcCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAcc                                                                 >  2:1574454/1‑100 (MQ=255)
                                   ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGGCATTCTGAAACCTGAAc                                                                  >  2:2059295/1‑98 (MQ=255)
                                    cAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTc                                                               >  2:1283225/1‑100 (MQ=255)
                                     aaTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                              <  1:1017521/100‑1 (MQ=255)
                                     aaTTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCa                                                              >  1:517253/1‑100 (MQ=255)
                                       ttCAACAGAAAGAATCAAGAAGAGGGGTTTGTATTTTGTCATAGGACAGTTGATGAGCTAGCATGGTTGCTCATTGTCATGCTGTAACCTGAACCtggta                                                            >  2:1764182/1‑96 (MQ=255)
                                         cagcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaa                                                          >  1:2174333/1‑100 (MQ=255)
                                         cagcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaa                                                          <  2:1520324/100‑1 (MQ=255)
                                           gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                        <  1:1041286/100‑1 (MQ=255)
                                           gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                        <  1:2210710/100‑1 (MQ=255)
                                           gcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                        >  1:2539553/1‑100 (MQ=255)
                                           cttccgatcttGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                        <  1:1450817/90‑1 (MQ=255)
                                            cagAACGAATCACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTg                                                       <  2:1999374/100‑1 (MQ=255)
                                             agAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTgg                                                      <  1:2554936/100‑1 (MQ=255)
                                              gAACGAATGACTAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                     <  1:1252427/100‑1 (MQ=255)
                                                aCGAATGACGAAGCGCATTTGGTATATTGTCTGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTc                                                   <  1:757881/100‑1 (MQ=255)
                                                  gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCga                                                 <  2:466140/100‑1 (MQ=255)
                                                  gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCga                                                 <  1:2148745/100‑1 (MQ=255)
                                                   aaTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                <  2:1044960/100‑1 (MQ=255)
                                                   aaTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                >  1:1707547/1‑100 (MQ=255)
                                                    aTGACGAAGCGCATTTGGTATATTGTCAGGGGGTGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgaga                                               <  1:1048514/100‑1 (MQ=255)
                                                    aTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgaga                                               <  1:53341/100‑1 (MQ=255)
                                                    aTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgaga                                               <  2:1033365/100‑1 (MQ=255)
                                                     tggcGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTTCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGAt                                              <  2:1744394/97‑1 (MQ=255)
                                                     tGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACtagatcggaag                                              >  2:1450803/1‑90 (MQ=255)
                                                     tGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGATTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGAt                                              >  1:225501/1‑100 (MQ=255)
                                                       atgttgcgcATTTGTTATATGGACAGGGGTCGTTAGATGGGCAAGCATGGTTGCTCGTTGCAATTGTGAAACGCGAACTTCATAAACTGGATCGAGAGaa                                            <  2:808918/95‑1 (MQ=255)
                                                        cGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaac                                           >  2:1389817/1‑100 (MQ=255)
                                                         gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaaca                                          >  2:2281110/1‑100 (MQ=255)
                                                         gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaaca                                          <  1:1493507/100‑1 (MQ=255)
                                                             cgcATTTGGTATATTGTCAGGTGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCa                                      <  2:2941050/100‑1 (MQ=255)
                                                             cgcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCa                                      <  2:2390824/100‑1 (MQ=255)
                                                                 tttGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAAtt                                  <  2:2358494/100‑2 (MQ=255)
                                                                    ggTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCAt                               >  2:52179/1‑100 (MQ=255)
                                                                     gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATa                              >  1:1639997/1‑100 (MQ=255)
                                                                     gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATa                              >  2:1133755/1‑100 (MQ=255)
                                                                     gTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAAACTGAAACTCATAAACTGGCTCGAGATAAAAACCAAATGCATa                              >  1:1639775/1‑100 (MQ=255)
                                                                      tataTTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATaa                             >  2:2562555/1‑100 (MQ=255)
                                                                           tGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTg                        <  2:1013712/100‑1 (MQ=255)
                                                                            gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTgg                       <  2:603293/100‑1 (MQ=255)
                                                                             tCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                      >  2:1817713/1‑100 (MQ=255)
                                                                              cAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGa                     >  2:1850858/1‑100 (MQ=255)
                                                                               aGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGaa                    >  2:194810/1‑100 (MQ=255)
                                                                                gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATAGGGaaa                   >  1:2582845/1‑100 (MQ=255)
                                                                                  gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAAt                 >  1:2919644/1‑100 (MQ=255)
                                                                                     cGTTTGGTGAGCAAGCATGGTTGCTAATTGCTGTTTTGAAAGCTGAGCCTCATAAACTGGCTCGAGATAACAAGGACATGCATAAAATTGGGAACATGGc              >  2:2423091/1‑100 (MQ=255)
                                                                                     cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc              <  2:700709/100‑1 (MQ=255)
                                                                                      gTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCa             <  2:285658/100‑1 (MQ=255)
                                                                                        gcaatgaGCAACCATGCTTGCTCATTTCCATTCTGAAACCTGAACCTCATAAACTAGCACGAAATAACAACCAAATGCATAAAATTGGGAAAATGGCAGa           >  1:804114/4‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa          >  1:377765/3‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa          >  2:742828/3‑100 (MQ=255)
                                                                                         caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGACTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa          >  1:944993/3‑100 (MQ=255)
                                                                                          gatgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAg         <  2:2071350/100‑1 (MQ=255)
                                                                                          gatgaGCAAGCATGGTGGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAAGTGGGAAAATGGCGGAAg         >  2:2275359/1‑100 (MQ=255)
                                                                                           atgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGACTCGAAATAACAACCAAATGCATAAAATTGGGAAAATAGCAGAAg         >  1:2228284/1‑99 (MQ=255)
                                                                                            tgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTAGGAAAATGGCAGAAGCg       >  2:1122424/1‑100 (MQ=255)
                                                                                              agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGtt     >  1:1718951/1‑100 (MQ=255)
                                                                                              agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGtt     <  1:1260801/100‑1 (MQ=255)
                                                                                                 aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCg  >  1:2475792/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
GCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCG  >  NC_000913/252490‑252686

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