Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,197,853 T→C 27.1% N129D (AAT→GAT)  ymfE ← e14 prophage; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,8530TC27.1% 142.6 / 59.8 92N129D (AAT→GAT) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (30/37);  new base C (11/14);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGA  >  NC_000913/1197756‑1197947
                                                                                                 |                                                                                              
gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGATGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                              <  2:1429648/100‑1 (MQ=255)
gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                              <  1:652749/100‑1 (MQ=255)
  aaaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTc                                                                                            <  1:2199802/100‑1 (MQ=255)
    aTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCAt                                                                                          >  2:2071824/1‑100 (MQ=255)
    aTTCAAAAAATATACAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCAt                                                                                          <  2:1563309/100‑1 (MQ=255)
      tCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTa                                                                                        >  1:386527/1‑100 (MQ=255)
       cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTAc                                                                                       >  2:292773/1‑100 (MQ=255)
       cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTAc                                                                                       >  2:247854/1‑100 (MQ=255)
       cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTAc                                                                                       >  1:1906989/1‑100 (MQ=255)
        aaaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACa                                                                                      >  2:123259/1‑100 (MQ=255)
               atCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGa                                                                                <  2:199464/99‑1 (MQ=255)
                tCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAg                                                                               <  1:560950/99‑1 (MQ=255)
                tCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGc                                                                              >  1:1569378/1‑100 (MQ=255)
                 caaaaaaTAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCt                                                                             <  2:1168676/98‑1 (MQ=255)
                  cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCtt                                                                            >  1:501267/1‑100 (MQ=255)
                    aaaaTAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCa                                                                          <  1:723492/100‑1 (MQ=255)
                    aaaaTAGTATACTTTATTCCAGAAGAATTCAATATAATGTTTATCTTCAATTTTTCTTACTTCAGATTAATATAGATTACTCATTACATTGTGAGCTTCa                                                                          >  2:1025086/1‑100 (MQ=255)
                        tAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCtt                                                                      <  1:497823/100‑1 (MQ=255)
                         aGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCttt                                                                     >  1:480673/1‑100 (MQ=255)
                             taCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCtttattt                                                                 >  1:1784610/1‑100 (MQ=255)
                             taCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCtttattt                                                                 >  2:1039742/1‑100 (MQ=255)
                               cTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTaa                                                               >  2:1231338/1‑100 (MQ=255)
                               cTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTaa                                                               >  2:933000/1‑100 (MQ=255)
                                 ttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAtt                                                             <  2:466958/100‑1 (MQ=255)
                                 ttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAtt                                                             <  2:481976/100‑1 (MQ=255)
                                   aTTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTTTTTTTTTTTTTTAATATAGATTGCTCATTACATTGTAAGCTTGATTTTTATTTAAtttt                                                           >  2:554777/1‑100 (MQ=255)
                                    ttCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTc                                                          <  2:1472223/100‑1 (MQ=255)
                                    ttCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTc                                                          >  1:761677/1‑100 (MQ=255)
                                    ttCCAGAAGAGTTAAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTc                                                          <  1:1073255/100‑1 (MQ=255)
                                      ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                        >  1:2078899/1‑100 (MQ=255)
                                      ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                        >  1:1970635/1‑100 (MQ=255)
                                      ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                        <  2:160320/100‑1 (MQ=255)
                                      ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                        <  1:1844249/100‑1 (MQ=255)
                                       cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                        <  1:925229/99‑1 (MQ=255)
                                       cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                       >  1:75706/1‑100 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                     <  2:162280/100‑1 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                     <  2:278229/100‑1 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                     <  2:2133597/100‑1 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                     <  1:1010833/100‑1 (MQ=255)
                                           agagTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGAc                                                   <  1:193296/100‑1 (MQ=255)
                                            gagTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATATTCTGTTGAc                                                   <  2:987403/99‑1 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTATTTAAGGGTAATATAGATTGCTCATTATATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                >  2:1553170/1‑100 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                >  1:628595/1‑100 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                <  2:1395143/100‑1 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAg                                              >  2:1475684/1‑99 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                             <  2:1232753/100‑1 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                             <  1:590483/100‑1 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTACTGTTGACTCCAGc                                             >  2:1526832/1‑100 (MQ=255)
                                                   atatAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACTTTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGct                                            >  2:1545605/1‑99 (MQ=255)
                                                   atatAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGctc                                           <  2:1167811/100‑1 (MQ=255)
                                                     ataATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCtctc                                         >  1:258233/1‑100 (MQ=255)
                                                       aaTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCg                                       <  1:1200510/100‑1 (MQ=255)
                                                       aaTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCg                                       <  1:1501263/100‑1 (MQ=255)
                                                        aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGctct                                          >  1:56607/1‑96 (MQ=255)
                                                        aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGt                                      <  1:167671/100‑1 (MQ=255)
                                                         tGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTg                                     >  1:906499/1‑100 (MQ=255)
                                                          gTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGa                                    <  1:1145995/100‑1 (MQ=255)
                                                          gTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGa                                    <  1:401463/100‑1 (MQ=255)
                                                           tttGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTGCATCTTTATTTCATTTTATATTGACTCAAGCTCTCCGTGAt                                   >  2:1125836/1‑100 (MQ=255)
                                                              gTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAAc                                <  1:2200659/100‑1 (MQ=255)
                                                              ctCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAAc                                <  2:29645/99‑1 (MQ=255)
                                                                cTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAAcag                              >  1:29645/1‑98 (MQ=255)
                                                                cTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                              >  1:1400399/1‑100 (MQ=255)
                                                                cTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACg                               >  1:858736/1‑99 (MQ=255)
                                                                cTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACCCCAGCTCTCCGTGATAACgg                              <  1:1854605/100‑1 (MQ=255)
                                                                cTTCAATTTTTCTAACTTCAGGGTAATATAGCTCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                              >  1:1514018/1‑100 (MQ=255)
                                                                 ttCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGt                             <  2:1414378/100‑1 (MQ=255)
                                                                 ttCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGt                             <  2:1474311/100‑1 (MQ=255)
                                                                  tCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGtt                            <  1:1951521/100‑1 (MQ=255)
                                                                  tCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGtt                            >  2:1742932/1‑100 (MQ=255)
                                                                    aaTTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGtttt                          <  1:79203/100‑1 (MQ=255)
                                                                    aaTTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGtttt                          <  1:509057/100‑1 (MQ=255)
                                                                      tttttCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTtat                        <  1:1683204/100‑1 (MQ=255)
                                                                       ttttCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTGTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTtata                       <  1:1657015/100‑1 (MQ=255)
                                                                         ttCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAAt                     <  2:1740449/100‑1 (MQ=255)
                                                                              aCTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGAt                <  1:644944/100‑1 (MQ=255)
                                                                              aCTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCAGCTTTATTTAATTTTCTGTTGTTTCCAGCTCTCCTTGATAACGGTTTTATAATTAGAt                <  2:756507/100‑1 (MQ=255)
                                                                               nttCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg               >  1:1157314/2‑100 (MQ=255)
                                                                                ttCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATgn              <  2:477743/100‑2 (MQ=255)
                                                                                 tCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt             >  2:1256327/1‑100 (MQ=255)
                                                                                   aGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTa           <  2:2078885/100‑1 (MQ=255)
                                                                                   aGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTa           >  2:304655/1‑100 (MQ=255)
                                                                                     ggTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATc         <  2:953602/100‑1 (MQ=255)
                                                                                       tAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGAGTCCAGCTCTCCGTTATAACGGTTTTATAATTAGATGCTTATccc       >  2:2112450/1‑100 (MQ=255)
                                                                                       tAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATccc       >  1:558283/1‑100 (MQ=255)
                                                                                         atatAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAATGGTTTTATTATTAGATGCTTATcc        <  2:910030/97‑1 (MQ=255)
                                                                                          tataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCaaa    >  1:375144/1‑100 (MQ=255)
                                                                                          tataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCaaa    >  2:2086983/1‑100 (MQ=255)
                                                                                           ataGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg   <  1:114586/100‑1 (MQ=255)
                                                                                            taGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTTATAACGGTTTTATAATTAGATGCTTATCCCAAAg   <  2:469932/99‑1 (MQ=255)
                                                                                            taGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg   <  1:1273226/99‑1 (MQ=255)
                                                                                            taGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGa  >  1:1238165/1‑100 (MQ=255)
                                                                                                 |                                                                                              
GAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGA  >  NC_000913/1197756‑1197947

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