Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 777,644 A→G intergenic (+37/‑96) tolA → / → tolB membrane anchored protein in TolA‑TolQ‑TolR complex/periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913777,6440AG96.2% 272.5 / ‑0.9 79intergenic (+37/‑96)tolA/tolBmembrane anchored protein in TolA‑TolQ‑TolR complex/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (35/41);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTAGCATTTGGTTT  >  NC_000913/777509‑777771
                                                                                                                                       |                                                                                                                               
gCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCg                                                                                                                              >  2:1411047/1‑139 (MQ=255)
  tGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGt                                                                                                                            >  2:28553/1‑139 (MQ=255)
  tGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGt                                                                                                                            >  1:625762/1‑139 (MQ=255)
  tGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGt                                                                                                                            >  2:120024/1‑139 (MQ=255)
    cGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTg                                                                                                                          >  1:435828/1‑139 (MQ=255)
       tGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCa                                                                                                                       <  2:435828/139‑1 (MQ=255)
       tGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCa                                                                                                                       <  1:28553/139‑1 (MQ=255)
        ggCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCAc                                                                                                                      <  1:459202/139‑1 (MQ=255)
         gcagcaGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACg                                                                                                                     <  2:1702583/139‑1 (MQ=255)
         gcagcaGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACg                                                                                                                     <  1:120024/139‑1 (MQ=255)
          cagcagCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACgg                                                                                                                    <  2:1207318/139‑1 (MQ=255)
           agcagcTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTg                                                                                                                          >  2:1548981/1‑132 (MQ=255)
           agcagcTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTg                                                                                                                          <  1:1548981/132‑1 (MQ=255)
             cagcTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGt                                                                                                                 <  2:523645/139‑1 (MQ=255)
               gcTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTc                                                                                                               <  1:1411047/139‑1 (MQ=255)
               gcTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTc                                                                                                               <  1:1454026/139‑1 (MQ=255)
               gcTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTc                                                                                                               <  2:196692/139‑1 (MQ=255)
                    aCTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTa                                                                                                          <  2:785646/139‑1 (MQ=255)
                           aaGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTATTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTtgt                                                                                                   <  2:1003842/139‑1 (MQ=255)
                               tCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTAt                                                                                               >  1:262664/1‑139 (MQ=255)
                                cccGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGt                                                                                                                            >  2:491771/1‑109 (MQ=255)
                                cccGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGt                                                                                                                            <  1:491771/109‑1 (MQ=255)
                                  cGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTAtttt                                                                                            >  2:1341245/1‑139 (MQ=255)
                                        caccaAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCAc                                                                                                                      >  1:706029/1‑107 (MQ=255)
                                        caccaAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCAc                                                                                                                      <  2:706029/107‑1 (MQ=255)
                                        caccaAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGtt                                                                                                       <  2:1639383/122‑1 (MQ=255)
                                        caccaAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGCAGtt                                                                                                       >  1:1639383/1‑122 (MQ=255)
                                          ccaAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGtt                                                                                    >  1:652318/1‑139 (MQ=255)
                                           caAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTa                                                                                   >  1:343336/1‑139 (MQ=255)
                                                  ggCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCt                                                                            <  2:343336/139‑1 (MQ=255)
                                                        atatGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTaaa                                                                      <  2:262664/139‑1 (MQ=255)
                                                        atatGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTaaa                                                                      <  1:952038/139‑1 (MQ=255)
                                                                 gtTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCg                                                                                                                              >  2:953565/1‑74 (MQ=255)
                                                                 gtTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCg                                                                                                                              <  1:953565/74‑1 (MQ=255)
                                                                 gtTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACg                                                                                                                     <  1:1313620/83‑1 (MQ=255)
                                                                 gtTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACg                                                                                                                     >  2:1313620/1‑83 (MQ=255)
                                                                 gtTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTgg                                                                                                            <  1:727866/92‑1 (MQ=255)
                                                                 gtTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTgg                                                                                                            >  2:727866/1‑92 (MQ=255)
                                                                 gtTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTgg                                                             >  2:510577/1‑139 (MQ=255)
                                                                 gtTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTgg                                                             >  2:1310714/1‑139 (MQ=255)
                                                                 gtTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGt                                                                                                                            >  1:934767/1‑76 (MQ=255)
                                                                 gtTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGt                                                                                                                            <  2:934767/76‑1 (MQ=255)
                                                                  ttCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGt                                                                                                                 >  1:544545/1‑86 (MQ=255)
                                                                  ttCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGt                                                                                                                 <  2:544545/86‑1 (MQ=255)
                                                                  ttCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTggg                                                            >  1:1118734/1‑139 (MQ=255)
                                                                   tCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTg                                                                           <  2:1429312/123‑1 (MQ=255)
                                                                   tCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTg                                                                           >  1:1429312/1‑123 (MQ=255)
                                                                     aaaaaCGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCa                                                         <  2:906333/139‑1 (MQ=255)
                                                                          cgcgCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGt                                                    <  1:1310714/139‑1 (MQ=255)
                                                                          cgcgCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGt                                                    <  2:625762/139‑1 (MQ=255)
                                                                          cgcgCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGt                                                    <  1:805302/139‑1 (MQ=255)
                                                                          cgcgCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGt                                                    <  1:697986/139‑1 (MQ=255)
                                                                                   ggACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAg                                           <  2:652318/139‑1 (MQ=255)
                                                                                         cAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGt                                                    >  2:980463/1‑124 (MQ=255)
                                                                                         cAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGt                                                    >  2:1361907/1‑124 (MQ=255)
                                                                                         cAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGt                                                    <  1:980463/124‑1 (MQ=255)
                                                                                         cAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGt                                                    <  1:1361907/124‑1 (MQ=255)
                                                                                           aaCCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGatat                                   >  1:844645/1‑139 (MQ=255)
                                                                                              cGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCa                                                                                                                       <  1:452125/52‑1 (MQ=255)
                                                                                              cGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCa                                                                                                                       >  2:796065/1‑52 (MQ=255)
                                                                                              cGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCa                                                                                                                       >  2:452125/1‑52 (MQ=255)
                                                                                              cGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCa                                                                                                                       <  1:796065/52‑1 (MQ=255)
                                                                                              cGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGt                                                                                                                 <  1:729354/58‑1 (MQ=255)
                                                                                              cGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGt                                                                                                                 >  2:729354/1‑58 (MQ=255)
                                                                                                tAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATAtgatga                              >  2:773145/1‑139 (MQ=255)
                                                                                                 aaTCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTc                                                                             <  1:105414/91‑1 (MQ=255)
                                                                                                 aaTCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTc                                                                             >  2:105414/1‑91 (MQ=255)
                                                                                                 aaTCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGaa                             >  2:1442223/1‑139 (MQ=255)
                                                                                                          gTTGACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCAtt                    <  2:844645/139‑1 (MQ=255)
                                                                                                             gACTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACg                 <  2:1644420/139‑1 (MQ=255)
                                                                                                              aCTGTTCGGACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGa                >  2:1210132/1‑139 (MQ=255)
                                                                                                                 gTTCGGACGGTCAACATCAGGCGCCGGTTGCCACggg                                                                                                                   <  2:1316798/37‑1 (MQ=255)
                                                                                                                 gTTCGGACGGTCAACATCAGGCGCCGGTTGCCACggg                                                                                                                   >  1:1316798/1‑37 (MQ=255)
                                                                                                                 gTTCGGACGGTCAACATCAGGCGCCGGTTGCCACg                                                                                                                     >  1:1354242/1‑35 (MQ=255)
                                                                                                                 gTTCGGACGGTCAACATCAGGCGCCGGTTGCCACg                                                                                                                     <  2:1354242/35‑1 (MQ=255)
                                                                                                                     ggACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGtttt                                                                                                     >  2:223351/1‑47 (MQ=255)
                                                                                                                     ggACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGtttt                                                                                                     <  1:223351/47‑1 (MQ=255)
                                                                                                                     ggACGGTCAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTAGCAt         <  1:1442223/139‑1 (MQ=255)
                                                                                                                            cAACATCAGGCGCCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTAGCATTTGGttt  >  1:969625/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                               
GCTGCGTTGGCAGCAGCTAAACTTGCGAAGATCCCGAAACCACCAAGCCAGGCAGTATATGAAGTGTTCAAAAACGCGCCATTGGACTTCAAACCGTAATCGCGATGTTGACTGTTCGGACGGTCAACATCAGGCACCGGTTGCCACGGGGTTCTGGTAGTTTTGTGTATTTTAGTTTGTTAACATTCTGCTAAATTATCGTGGGCCATCGGTCCAGATAAGGGAGATATGATGAAGCAGGCATTACGAGTAGCATTTGGTTT  >  NC_000913/777509‑777771

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: