Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,584,183 A→G V860A (GTG→GCG)  hsdR ← type I restriction enzyme EcoKI endonuclease component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,584,1830AG97.6% 282.7 / ‑2.5 82V860A (GTG→GCG) hsdRtype I restriction enzyme EcoKI endonuclease component
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (35/45);  total (36/46)
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.90e-01

CGCTGGAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCT  >  NC_000913/4584048‑4584316
                                                                                                                                       |                                                                                                                                     
cGCTGGAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCgg                                                                                                                                   >  1:1045107/1‑140 (MQ=255)
 gCTGGAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGa                                                                                                                                  >  2:286365/1‑140 (MQ=255)
    ggAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGc                                                                                                                               >  1:930657/1‑140 (MQ=255)
     gAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCa                                                                                                                              >  2:12664/1‑140 (MQ=255)
     gAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCa                                                                                                                              >  1:1290025/1‑140 (MQ=255)
             ccaccaGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCg                                                                                                                      <  2:1090817/140‑1 (MQ=255)
             ccaccaGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCg                                                                                                                      <  1:865844/140‑1 (MQ=255)
             ccaccaGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCg                                                                                                                      <  1:766376/140‑1 (MQ=255)
               accaGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGAc                                                                                                                    <  2:862921/140‑1 (MQ=255)
                ccaGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACg                                                                                                                   <  1:609496/140‑1 (MQ=255)
                    ttgttCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTCCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGctct                                                                                                               <  2:1223860/140‑1 (MQ=255)
                       ttCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCa                                                                                                            <  2:1190474/140‑1 (MQ=255)
                       ttCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCa                                                                                                            <  1:880878/140‑1 (MQ=255)
                       ttCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCa                                                                                                            <  1:580573/140‑1 (MQ=255)
                       ttCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCa                                                                                                            <  2:706136/140‑1 (MQ=255)
                          aTGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCg                                                                                                         <  2:1146520/140‑1 (MQ=255)
                          aTGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCg                                                                                                         >  2:416045/1‑140 (MQ=255)
                           tGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGt                                                                                                        >  1:1284026/1‑140 (MQ=255)
                            ggCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTg                                                                                                       <  2:12680/140‑1 (MQ=255)
                            ggCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTg                                                                                                       <  2:1288747/140‑1 (MQ=255)
                            ggCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTg                                                                                                       >  1:141780/1‑140 (MQ=255)
                            ggCTGTACTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTg                                                                                                       >  2:1260565/1‑140 (MQ=255)
                             gCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGc                                                                                                      <  2:769333/140‑1 (MQ=255)
                                 tgCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGc                                                                                                          >  2:724138/1‑132 (MQ=255)
                                 tgCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGc                                                                                                          <  1:724138/132‑1 (MQ=255)
                                 tgCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTa                                                                                                  <  1:1075148/140‑1 (MQ=255)
                                 tgCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTa                                                                                                  <  1:1126315/140‑1 (MQ=255)
                                 tgCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTa                                                                                                  <  1:1203839/140‑1 (MQ=255)
                                 tgCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTa                                                                                                  >  1:11299/1‑140 (MQ=255)
                                   cTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGa                                                                                                >  2:1328218/1‑140 (MQ=255)
                                        cAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGc                                                                                                          >  2:481779/1‑125 (MQ=255)
                                        cAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGc                                                                                                          <  1:481779/125‑1 (MQ=255)
                                        cAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCg                                                                                           >  1:839669/1‑140 (MQ=255)
                                        cAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCg                                                                                           >  1:49843/1‑140 (MQ=255)
                                           aaCTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGaca                                                                                        <  2:444497/140‑1 (MQ=255)
                                            aCTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGacac                                                                                       <  1:105319/140‑1 (MQ=255)
                                                cGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATAGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGATGTCGACCCTCTCGAGCGTGCTGTAGATATCGACACAGTc                                                                                   >  2:1236832/1‑140 (MQ=255)
                                                  gCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCaa                                                                                 >  1:445775/1‑140 (MQ=255)
                                                     aTCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCaaaaa                                                                              >  2:556631/1‑140 (MQ=255)
                                                            tCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTaaa                                                                       <  2:668126/140‑1 (MQ=255)
                                                             cGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAg                                                                      >  2:430976/1‑140 (MQ=255)
                                                                 gATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTgg                                                                  <  2:11299/140‑1 (MQ=255)
                                                                  aTTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGt                                                                 >  1:289262/1‑140 (MQ=255)
                                                                        taGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTAtt                                                           >  1:342785/1‑140 (MQ=255)
                                                                           gTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCAc                                                        <  1:393902/140‑1 (MQ=255)
                                                                             ttCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCt                                                      <  2:342785/140‑1 (MQ=255)
                                                                             ttCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCt                                                      <  2:1045107/140‑1 (MQ=255)
                                                                              tCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTc                                                     >  1:140798/1‑140 (MQ=255)
                                                                              tCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGGGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTc                                                     >  2:1392731/1‑140 (MQ=255)
                                                                              tCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTACACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTc                                                     >  1:316435/1‑140 (MQ=255)
                                                                                tGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTcc                                                    <  2:295967/139‑1 (MQ=255)
                                                                                tGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTcc                                                    >  1:295967/1‑139 (MQ=255)
                                                                                    tCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCa                                               <  1:1392731/140‑1 (MQ=255)
                                                                                     cGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCAt                                              >  2:110710/1‑140 (MQ=255)
                                                                                      ggTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATa                                             >  2:844691/1‑140 (MQ=255)
                                                                                            ttCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGcgcg                                       <  2:141780/140‑1 (MQ=255)
                                                                                             tCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGt                                      >  1:688159/1‑140 (MQ=255)
                                                                                                 tGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGgcg                                  <  2:289262/140‑1 (MQ=255)
                                                                                                  gACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGgcgc                                 >  2:82052/1‑140 (MQ=255)
                                                                                                  gACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGgcgc                                 <  1:1236832/140‑1 (MQ=255)
                                                                                                  gACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGgcgc                                 <  1:1183298/140‑1 (MQ=255)
                                                                                                     cAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGc                              >  2:541148/1‑140 (MQ=255)
                                                                                                     cAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGc                              >  1:358910/1‑140 (MQ=255)
                                                                                                      aGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGcc                             <  2:140798/140‑1 (MQ=255)
                                                                                                      aGCGTTTGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGcc                             >  2:1039426/1‑140 (MQ=255)
                                                                                                            tGCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCAt                       <  2:964835/140‑1 (MQ=255)
                                                                                                             gCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATc                      >  2:983166/1‑140 (MQ=255)
                                                                                                             gCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATc                      <  2:316435/140‑1 (MQ=255)
                                                                                                             gCAGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATc                      <  1:286365/140‑1 (MQ=255)
                                                                                                               aGTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTg                    <  1:12664/140‑1 (MQ=255)
                                                                                                                gTTCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGt                   <  1:541148/140‑1 (MQ=255)
                                                                                                                  tCCACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTc                 >  2:304927/1‑140 (MQ=255)
                                                                                                                   ccACCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCg                >  1:221986/1‑140 (MQ=255)
                                                                                                                     aCCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTa              <  1:882186/140‑1 (MQ=255)
                                                                                                                     aCCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTa              <  2:688159/140‑1 (MQ=255)
                                                                                                                     aCCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTa              <  1:491198/140‑1 (MQ=255)
                                                                                                                     aCCTTCGGGCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTa              <  1:1039426/140‑1 (MQ=255)
                                                                                                                           gggCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAAt        >  2:373849/1‑140 (MQ=255)
                                                                                                                            ggCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATg       <  1:1260565/140‑1 (MQ=255)
                                                                                                                            ggCGCACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGACATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATg       <  2:839669/140‑1 (MQ=255)
                                                                                                                               gcACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCgg    <  1:249695/140‑1 (MQ=255)
                                                                                                                               gcACCACCGCCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCgg    <  1:1042577/140‑1 (MQ=255)
                                                                                                                                 accaccgccGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCt  <  2:1284026/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CGCTGGAGCTTCGCCACCAGTTGTTCATGGCTGTGCTCGGCAAAACTGCGGCCATCCGCTTCGGTGATTTTATAGGTTTCTGAATCGGTAATTTCATTGACCAGCGTTTGCAGTTCCACCTTCGGGCGCACCACCACCGGACGCATGGTGTCGACGCTCTCCAGCGTGCTGTAGATATCGACACAGTCAAAAATCTTAAAGCTGGTTTTATTCACCTCCGGGCATAAGCGCGTGGCGCGGCCTTTCATCTGTTCGTACAGAATGCGGCT  >  NC_000913/4584048‑4584316

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