Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,279,849 A→G intergenic (+21/‑131) ghxP → / → yjcE guanine/hypoxanthine permease, high affinity, guanine/hypoxanthine:H+ symporter/putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,279,8490AG98.7% 268.8 / ‑5.7 75intergenic (+21/‑131)ghxP/yjcEguanine/hypoxanthine permease, high affinity, guanine/hypoxanthine:H+ symporter/putative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (36/38);  minor base T (0/1);  total (36/39)
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

TCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCAT  >  NC_000913/4279713‑4279981
                                                                                                                                        |                                                                                                                                    
tCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCg                                                                                                                                    <  2:328596/139‑1 (MQ=255)
 cgcgACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                                   >  2:1786718/1‑139 (MQ=255)
  gcgACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCggg                                                                                                                                  >  2:638425/1‑139 (MQ=255)
       tctGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCa                                                                                                                             >  1:899572/1‑139 (MQ=255)
                  gtcgtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCg                                                                                                                  <  1:1626162/139‑1 (MQ=255)
                     gtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCa                                                                                                                             >  2:99621/1‑125 (MQ=255)
                     gtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCa                                                                                                                             <  1:99621/125‑1 (MQ=255)
                     gtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGttt                                                                                                               >  2:1000823/1‑139 (MQ=255)
                     gtcTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGttt                                                                                                               >  2:286966/1‑139 (MQ=255)
                            ttCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTAtttt                                                                                                        >  1:767071/1‑139 (MQ=255)
                                   gAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGca                                                                                                 <  2:1852659/139‑1 (MQ=255)
                                     aTGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGcaca                                                                                               <  2:899572/139‑1 (MQ=255)
                                     aTGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGcaca                                                                                               <  2:136429/139‑1 (MQ=255)
                                           cAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgt                                                                                         <  1:286966/139‑1 (MQ=255)
                                              gTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAg                                                                                      <  1:638425/139‑1 (MQ=255)
                                                      tctgTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGTAATCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCg                                                                              <  1:1359761/136‑1 (MQ=255)
                                                      tCGGTACGGTCGTTATCGCCGTTGCGCTGGATACCTTCTAGGCGTGTTGCTGGGCTATCTAATCTTTCGTCGGATTCTGAAAGCGGGTGGCAGGGGCTGCCCGCTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCg                                                                              <  1:1913950/139‑1 (MQ=255)
                                                       cGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGa                                                                             >  1:1623459/1‑139 (MQ=255)
                                                        ggTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                                   <  2:687776/84‑1 (MQ=255)
                                                        ggTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                                   >  1:687776/1‑84 (MQ=255)
                                                         gTACGGTCGTTATCGCCGTGGCGTTGGTGACCTACTGTGCGGGTGGCTGGGCTATCTAATCTTTCATCTGATTCTGAAAGCGGGCGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGata                                                                           >  2:1913950/1‑139 (MQ=255)
                                                         gTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGt                                                                                                                                 >  2:173757/1‑85 (MQ=255)
                                                         gTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGt                                                                                                                                 <  1:173757/85‑1 (MQ=255)
                                                               tCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTAt                                                                                                           >  1:1071674/1‑101 (MQ=255)
                                                               tCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTAt                                                                                                           <  2:1071674/101‑1 (MQ=255)
                                                                  ttATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTtctc                                                                                                     >  2:239606/1‑104 (MQ=255)
                                                                  ttATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTtctc                                                                                                     <  1:239606/104‑1 (MQ=255)
                                                                    aTCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTAt                                                                                                           <  2:1222428/96‑1 (MQ=255)
                                                                    aTCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTAt                                                                                                           >  1:1222428/1‑96 (MQ=255)
                                                                        ccGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgt                                                                                         >  2:1967931/1‑110 (MQ=255)
                                                                        ccGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgt                                                                                         <  1:1967931/110‑1 (MQ=255)
                                                                         cGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATgag                                                           >  1:1842890/1‑139 (MQ=255)
                                                                           tGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTAt                                                                                                           >  2:597937/1‑89 (MQ=255)
                                                                           tGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTAt                                                                                                           <  1:597937/89‑1 (MQ=255)
                                                                                   gTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCa                                                 <  1:1000823/139‑1 (MQ=255)
                                                                                     gacagtCTATGCGGGTTGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGTTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACAGTGTGTGAGCTTATGCGATATTCTGAAAAAAATGAGATTTCAGGCATa                                               <  2:591590/134‑1 (MQ=255)
                                                                                     gACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATa                                               <  1:1786718/139‑1 (MQ=255)
                                                                                     gACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATa                                               <  2:1871302/139‑1 (MQ=255)
                                                                                       ccTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATTGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAAc                                             >  2:722568/1‑139 (MQ=255)
                                                                                        cTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATa                                               <  1:1298294/136‑1 (MQ=255)
                                                                                        cTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATa                                               >  2:1298294/1‑136 (MQ=255)
                                                                                             aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                                   >  2:1927169/1‑47 (MQ=255)
                                                                                             aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                                   >  1:1858646/1‑47 (MQ=255)
                                                                                             aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                                   <  1:1927169/47‑1 (MQ=255)
                                                                                             aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCgg                                                                                                                                   <  2:1858646/47‑1 (MQ=255)
                                                                                             aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTc                                          <  2:980733/136‑1 (MQ=255)
                                                                                             aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTc                                          >  1:980733/1‑136 (MQ=255)
                                                                                             aTGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCttt                                       >  2:967191/1‑139 (MQ=255)
                                                                                                cGGGTGGCTGGGCTATCTAATCTTTCATTGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTtctc                                    <  1:722568/139‑1 (MQ=255)
                                                                                                 gggTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCa                                   <  1:122158/139‑1 (MQ=255)
                                                                                                  ggTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTc                                                      <  2:616722/119‑1 (MQ=255)
                                                                                                  ggTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTc                                                      >  1:616722/1‑119 (MQ=255)
                                                                                                   gTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                                 <  2:1842890/139‑1 (MQ=255)
                                                                                                   gTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                                 >  1:1185626/1‑139 (MQ=255)
                                                                                                     ggCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTtctcctgtc                               >  1:591590/1‑134 (MQ=255)
                                                                                                      gCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACg                                            >  2:1817487/1‑125 (MQ=255)
                                                                                                      gCTGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACg                                            <  1:1817487/125‑1 (MQ=255)
                                                                                                        tGGGCTATCTAATCTTTCATCGGATTCTGAAAGCGTGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCATATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATa                            <  2:1185626/139‑1 (MQ=255)
                                                                                                          ggCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCg                                                                              <  2:1119336/87‑1 (MQ=255)
                                                                                                          ggCTATCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCg                                                                              >  1:1119336/1‑87 (MQ=255)
                                                                                                               tCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGcc                                                                                                                    >  1:169376/1‑44 (MQ=255)
                                                                                                               tCTAATCTTTCATCGGATTCTGAAAGCGGGTGGCAATGGCTGcc                                                                                                                    <  2:169376/44‑1 (MQ=255)
                                                                                                                        tCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgtg                                                                                        >  2:1247886/1‑63 (MQ=255)
                                                                                                                        tCATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATtgtgtg                                                                                        <  1:1247886/63‑1 (MQ=255)
                                                                                                                         cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCg                                                                                                                  >  1:469446/1‑36 (MQ=255)
                                                                                                                         cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCg                                                                                                                  <  2:469446/36‑1 (MQ=38)
                                                                                                                         cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTtct                                     <  1:1004432/113‑1 (MQ=255)
                                                                                                                         cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTtct                                     >  2:1004432/1‑113 (MQ=255)
                                                                                                                         cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                                 <  2:640227/117‑1 (MQ=255)
                                                                                                                         cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGa                                 >  1:640227/1‑117 (MQ=255)
                                                                                                                         cATCGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAg           >  1:2037359/1‑139 (MQ=255)
                                                                                                                            cGGATTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGa        <  2:1706071/139‑1 (MQ=255)
                                                                                                                               aTTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATg                                                             <  1:1746359/83‑1 (MQ=255)
                                                                                                                               aTTCTGAAAGCGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATg                                                             >  2:1746359/1‑83 (MQ=255)
                                                                                                                                  cTGAAAGCGGGTGACAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCAt  >  2:1592168/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                    
TCGCGACTCTGGTGATTGGTCGTCTGGTTTCCGGCGAATGGCGCAAGTTGAACATCGGTACGGTCGTTATCGCCGTGGCGCTGGTGACCTTCTATGCGGGTGGCTGGGCTATCTAATCTTTCATCGGATTCTGAAAACGGGTGGCAATGGCTGCCCGTTTTTATTTTCTCCGCACATTGTGTGAGCTTTTGCGATATTCTGAAAAAAATGAGAATTCAGGCATAACGTCTTTCTCAGATCATAAAAAACATCGCAAACAGGGAACGCAT  >  NC_000913/4279713‑4279981

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