Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,736,157 G→A *393* (TAG→TAA xylR → xylose divergent operon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,736,1570GA99.0% 309.5 / NA 95*393* (TAG→TAAxylRxylose divergent operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (45/49);  minor base C (0/1);  total (45/50)
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

AACCACCTTGTCGATCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAGCGATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAA  >  NC_000913/3736022‑3736287
                                                                                                                                       |                                                                                                                                   
aaCCACCTTGTCGATCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACga                                                                                                                                 <  2:71169/139‑1 (MQ=255)
aaCCACCTTGTCGATCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACga                                                                                                                                 <  1:181190/139‑1 (MQ=255)
  ccaccTTGTCGATCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACgaat                                                                                                                               >  2:58795/1‑137 (MQ=255)
     ccTTGTCGATCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACgaatgat                                                                                                                            <  1:176863/139‑6 (MQ=255)
         gTCGATCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATgaga                                                                                                                         >  2:699796/1‑139 (MQ=255)
               cAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTc                                                                                                                   <  1:58795/139‑1 (MQ=255)
                  tGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTACCGAATGATGAGAATTGTCGGc                                                                                                                <  2:100366/139‑1 (MQ=255)
                    agaTATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAAtgatga                                                                                                                           >  2:998664/1‑126 (MQ=255)
                    agaTATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTCTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAAtgatga                                                                                                                           <  1:998664/126‑1 (MQ=255)
                    agaTATCGCAAATGGGCGGTTATCCAACGCTGCAATACTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCgc                                                                                                              <  2:1054709/139‑2 (MQ=255)
                        atCGCAAATGTGCGGTTATCCATCGCTGCAGTATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               <  1:1002062/134‑1 (MQ=255)
                        atCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          <  1:1331497/139‑1 (MQ=255)
                        atCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGCATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               >  2:1002062/1‑134 (MQ=255)
                               aTGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              <  1:743969/128‑1 (MQ=255)
                               aTGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              >  2:743969/1‑128 (MQ=255)
                               aTGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTa                                                                                                   >  1:588963/1‑139 (MQ=255)
                                     ggTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTg                                                                                             <  1:1272667/139‑1 (MQ=255)
                                       ttATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATgag                                                                                                                          >  1:635092/1‑108 (MQ=255)
                                       ttATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATgag                                                                                                                          <  2:635092/108‑1 (MQ=255)
                                        tATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGt                                                                                          <  2:1379656/139‑1 (MQ=255)
                                         aTCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTg                                                                                         <  2:256506/139‑1 (MQ=255)
                                            cATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          >  1:1172752/1‑119 (MQ=255)
                                            cATCGCTGCAATATTCCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          <  2:1172752/119‑1 (MQ=255)
                                               cGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              <  1:1605186/112‑1 (MQ=255)
                                               cGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              >  2:1605186/1‑112 (MQ=255)
                                               cGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTc                                                                                   >  1:1353737/1‑139 (MQ=255)
                                                 cTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          >  2:78664/1‑114 (MQ=255)
                                                 cTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          <  1:78664/114‑1 (MQ=255)
                                                 cTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTGACGCTGCGTGGTTGTCgg                                                                                 >  1:1114634/1‑139 (MQ=255)
                                                 cTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCgg                                                                                 >  2:1041681/1‑139 (MQ=255)
                                                      atatTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATgcg                                                                            <  1:1041681/139‑1 (MQ=255)
                                                      atatTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATgcg                                                                            <  2:1351819/139‑1 (MQ=255)
                                                            cTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGa                                                                      <  2:171292/139‑1 (MQ=255)
                                                                  tGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          <  1:1378392/97‑1 (MQ=255)
                                                                  tGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTcaca                                                                                                          >  2:1378392/1‑97 (MQ=255)
                                                                   gTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCtt                                                               >  1:405631/1‑139 (MQ=255)
                                                                    tttttAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           <  2:526547/110‑1 (MQ=255)
                                                                    tttttAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCg                                                                                           >  1:526547/1‑110 (MQ=255)
                                                                      tttAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATc                                                            >  2:778028/1‑139 (MQ=255)
                                                                              aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGa                                                                      >  2:920111/1‑121 (MQ=255)
                                                                              aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGa                                                                      <  1:920111/121‑1 (MQ=255)
                                                                              aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGc                                                                  >  1:839510/1‑125 (MQ=255)
                                                                              aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGc                                                                  <  2:839510/125‑1 (MQ=255)
                                                                              aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTAc                                                    >  2:909345/1‑139 (MQ=255)
                                                                              aGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTAc                                                    >  1:426397/1‑139 (MQ=255)
                                                                                     gACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTATCGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAg                                             >  2:1447897/1‑139 (MQ=255)
                                                                                      acacGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAgg                                            <  1:1388998/139‑1 (MQ=255)
                                                                                          gacgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               <  1:638393/68‑1 (MQ=255)
                                                                                          gacgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               >  2:638393/1‑68 (MQ=255)
                                                                                          gacgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTa                                                     <  2:554968/126‑1 (MQ=255)
                                                                                          gacgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTa                                                     >  1:554968/1‑126 (MQ=255)
                                                                                            cgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACgaat                                                                                                                               >  1:859860/1‑47 (MQ=255)
                                                                                            cgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACgaat                                                                                                                               <  2:859860/49‑3 (MQ=255)
                                                                                            cgCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGAt                                      >  1:1415287/1‑139 (MQ=255)
                                                                                              ccAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTg                                                                       >  2:376656/1‑104 (MQ=255)
                                                                                              ccAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTg                                                                       <  1:376656/104‑1 (MQ=255)
                                                                                                aaaaGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATaaaa                                  <  1:778028/139‑1 (MQ=255)
                                                                                                          cgcgATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACcaaca                        <  1:909345/139‑1 (MQ=255)
                                                                                                           gcgATGTAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACcaacaa                       <  2:1415287/139‑1 (MQ=255)
                                                                                                                 tAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTATCGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTc                 <  1:1447897/139‑1 (MQ=255)
                                                                                                                 tAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATaa                                    <  1:664916/120‑1 (MQ=255)
                                                                                                                 tAAATAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATaa                                    >  2:664916/1‑120 (MQ=255)
                                                                                                                   aaTAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               <  2:1026072/43‑1 (MQ=255)
                                                                                                                   aaTAGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               >  1:1026072/1‑43 (MQ=255)
                                                                                                                      aGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCt            >  1:679392/1‑139 (MQ=255)
                                                                                                                      aGCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCt            <  2:588963/139‑1 (MQ=255)
                                                                                                                       gCGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTc           <  2:1496027/139‑1 (MQ=255)
                                                                                                                        cGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGc             <  2:803745/136‑1 (MQ=255)
                                                                                                                        cGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGc             >  1:803745/1‑136 (MQ=255)
                                                                                                                        cGAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCg          >  1:1060297/1‑139 (MQ=255)
                                                                                                                         gAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              >  2:1099236/1‑38 (MQ=255)
                                                                                                                         gAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGt                                                                                                              <  1:1099236/38‑1 (MQ=255)
                                                                                                                         gAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCg                                                          <  1:568710/90‑1 (MQ=38)
                                                                                                                         gAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCg                                                          >  2:568710/1‑90 (MQ=38)
                                                                                                                         gAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTa                                                     <  1:230711/95‑1 (MQ=38)
                                                                                                                         gAGGTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTa                                                     >  2:230711/1‑95 (MQ=255)
                                                                                                                           ggTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               <  1:1290853/35‑1 (MQ=255)
                                                                                                                           ggTCATGTTGTAACGAATGATGAGAATTGTCGGCg                                                                                                               >  2:1290853/1‑35 (MQ=255)
                                                                                                                           ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACg                                                                   >  1:1408262/1‑79 (MQ=39)
                                                                                                                           ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACg                                                                   <  2:1408262/79‑1 (MQ=39)
                                                                                                                           ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCa                                          >  2:314534/1‑104 (MQ=255)
                                                                                                                           ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCa                                          <  1:314534/104‑1 (MQ=255)
                                                                                                                           ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCt            >  2:1540063/1‑134 (MQ=255)
                                                                                                                           ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCt            <  1:1540063/134‑1 (MQ=255)
                                                                                                                           ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAg       >  2:1492823/1‑139 (MQ=255)
                                                                                                                           ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAg       >  2:1195123/1‑139 (MQ=255)
                                                                                                                           ggTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAg       >  2:332018/1‑139 (MQ=255)
                                                                                                                            gTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACcaac                         <  2:1627927/120‑1 (MQ=255)
                                                                                                                            gTCATGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACcaac                         >  1:1627927/1‑120 (MQ=255)
                                                                                                                               aTGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCg                                                 >  2:1566652/1‑93 (MQ=37)
                                                                                                                               aTGTTGTAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCg                                                 <  1:1566652/93‑1 (MQ=37)
                                                                                                                                tgttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCaa  <  2:679392/139‑1 (MQ=255)
                                                                                                                                tgttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCaa  <  2:1060297/139‑1 (MQ=255)
                                                                                                                                 gttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGc                                                                  >  2:23023/1‑74 (MQ=38)
                                                                                                                                 gttgtAACGAATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGc                                                                  <  1:23023/74‑1 (MQ=38)
                                                                                                                                       |                                                                                                                                   
AACCACCTTGTCGATCAATGAGATATCGCAAATGTGCGGTTATCCATCGCTGCAATATTTCTACTCTGTTTTTAAAAAAGCATATGACACGACGCCAAAAGAGTATCGCGATGTAAATAGCGAGGTCATGTTGTAGCGATGATGAGAATTGTCGGCGTCACATCAGGTAACGCTGCGTGGTTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACCCGCCAGGCATGATAAAACGCACCAACAACGCTTCAGGCTCGTAGCTCAA  >  NC_000913/3736022‑3736287

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