Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,037,018 G→T 5.3% P363P (CCC→CCA recJ ← ssDNA exonuclease, 5' ‑‑> 3'‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,037,0180GT 5.3% 297.8 / 2.1 84P363P (CCC→CCArecJssDNA exonuclease, 5' ‑‑> 3'‑specific
Reads supporting (aligned to +/- strand):  ref base G (36/43);  new base T (3/2);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCA  >  NC_000913/3036879‑3037154
                                                                                                                                           |                                                                                                                                        
cACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCg                                                                                                                                          <  2:1034644/140‑1 (MQ=255)
cACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCg                                                                                                                                          >  2:1546630/1‑140 (MQ=255)
cACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCg                                                                                                                                          >  1:1085017/1‑140 (MQ=255)
cACGCATATGCAGCCCCTGAATGGAACGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCg                                                                                                                                          >  1:1085774/1‑140 (MQ=255)
  cGCATAGGCAGCCCCTGAATGGTGCGACCGGAACCTTTCAGCGTACCGTCACCGGCGGGCGCAAACGCGATAACCGGACGGTGTAACCGCTCTTTGATGCGCGAAGCCAGAAAACCGACAACGCCCTGATGCCAGTCggg                                                                                                                                        <  2:365157/140‑1 (MQ=255)
    cATTTGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCTACACCGCCCTGATGCCATTCGGGGt                                                                                                                                      <  1:518326/140‑1 (MQ=255)
      tatGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGa                                                                                                                                    <  1:771164/140‑1 (MQ=255)
       atGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGAt                                                                                                                                   >  1:95011/1‑140 (MQ=255)
       atGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGAt                                                                                                                                   >  1:95010/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCTGGGTGATACATTg                                                                                                                             >  1:163280/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                             >  2:575280/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                             >  2:179687/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                             >  1:163271/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                             >  1:1440817/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                             >  1:1410998/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                             >  1:1058745/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                             >  1:1039111/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                             >  1:1059425/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                             >  1:1122430/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                             >  1:114260/1‑140 (MQ=255)
             ccccTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGAAGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTg                                                                                                                             >  2:680530/1‑140 (MQ=255)
                cTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCa                                                                                                                          <  1:1646497/140‑1 (MQ=255)
                cTGAATGGACCGACCGGAAGCATTCAGGGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCACTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTAGGGGTGATACATTGCCa                                                                                                                          <  2:1659665/140‑1 (MQ=255)
                  gAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGc                                                                                                                        >  2:68295/1‑140 (MQ=255)
                      ggAGCGACCGGAACCTTTCAGCGTACCGTCTCCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgc                                                                                                                    <  2:847328/140‑1 (MQ=255)
                      ggAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgc                                                                                                                    <  2:847425/140‑1 (MQ=255)
                      ggAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgc                                                                                                                    <  2:1291901/140‑1 (MQ=255)
                      ggAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgc                                                                                                                    <  2:1291045/140‑1 (MQ=255)
                      ggAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTAGGGGTGATACATTGCCAGCccgc                                                                                                                    <  2:847006/140‑1 (MQ=255)
                       gAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCTGGGTGATACATTGCCAGCccgcc                                                                                                                   >  1:754661/1‑140 (MQ=255)
                       gAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCTGGGTGATACATTGCCAGCccgcc                                                                                                                   >  1:753691/1‑140 (MQ=255)
                       gAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgcc                                                                                                                   >  1:93596/1‑140 (MQ=255)
                       gAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgcc                                                                                                                   >  2:1182399/1‑140 (MQ=255)
                       gAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgcc                                                                                                                   >  2:890713/1‑140 (MQ=255)
                        aGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCccgccg                                                                                                                  >  2:721392/1‑140 (MQ=255)
                             ccGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAg                                                                                                             <  1:225938/140‑1 (MQ=255)
                             ccGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAg                                                                                                             <  1:82727/140‑1 (MQ=255)
                             ccGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAg                                                                                                             <  1:82881/140‑1 (MQ=255)
                             ccGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAg                                                                                                             <  2:114260/140‑1 (MQ=255)
                             ccGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAg                                                                                                             >  2:875526/1‑140 (MQ=255)
                             ccGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAg                                                                                                             <  2:1059425/140‑1 (MQ=255)
                                  tccTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCgtgt                                                                                                        <  2:175192/139‑1 (MQ=255)
                                  aCCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCgtgt                                                                                                        <  2:175178/140‑1 (MQ=255)
                                     tttCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCAc                                                                                                     >  2:250092/1‑140 (MQ=255)
                                      ttCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACg                                                                                                    >  1:160462/1‑140 (MQ=255)
                                       tCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACgg                                                                                                   >  2:1256851/1‑140 (MQ=255)
                                       tCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACgg                                                                                                   >  2:1256724/1‑140 (MQ=255)
                                           cGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTgc                                                                                               >  1:838297/1‑140 (MQ=255)
                                           cGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTgc                                                                                               >  1:118831/1‑140 (MQ=255)
                                             tACCGTCACCTGCTTGCGCAAAATCTATAAAAGGACGGAGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTgcgc                                                                                             >  1:947346/1‑140 (MQ=255)
                                             tACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTgcgc                                                                                             >  2:710578/1‑140 (MQ=255)
                                                cGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTcc                                                                                          >  1:514246/1‑140 (MQ=255)
                                                 gTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCa                                                                                         <  2:1312644/140‑1 (MQ=255)
                                                     ccTGCTGGCGGATACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGATTACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGttt                                                                                     <  2:318396/140‑1 (MQ=255)
                                                       tgctgGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTtct                                                                                   <  1:593197/140‑1 (MQ=255)
                                                           ggCGCAAACGCGATAACCGGACGGTGAATACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGca                                                                               <  1:1445132/140‑1 (MQ=255)
                                                           ggCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGca                                                                               <  1:1041933/140‑1 (MQ=255)
                                                           ggCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGca                                                                               <  2:774289/140‑1 (MQ=255)
                                                           ggCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGca                                                                               <  1:850261/140‑1 (MQ=255)
                                                           ggCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGca                                                                               <  1:698292/140‑1 (MQ=255)
                                                           ggCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGca                                                                               <  1:1370066/140‑1 (MQ=255)
                                                           ggCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGca                                                                               <  1:1408861/140‑1 (MQ=255)
                                                           ggCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGca                                                                               <  2:655065/140‑1 (MQ=255)
                                                            gcgcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGcac                                                                              <  2:1181704/140‑1 (MQ=255)
                                                            gcgcAAACGCGATAACCGGACGGTGAAAACGCTCATTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGcac                                                                              <  1:1042682/140‑1 (MQ=255)
                                                             cgcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCTGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGcaca                                                                             <  2:754661/140‑1 (MQ=255)
                                                             cgcAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCTGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGcaca                                                                             <  2:753691/140‑1 (MQ=255)
                                                                   cgcgATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCa                                                                       <  2:313779/140‑1 (MQ=255)
                                                                       aTAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGc                                                                   >  2:68292/1‑140 (MQ=255)
                                                                        tAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCt                                                                  <  1:537012/140‑1 (MQ=255)
                                                                             ggACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAAt                                                             <  1:189439/140‑1 (MQ=255)
                                                                                cGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTg                                                          <  2:1600150/140‑1 (MQ=255)
                                                                                      aaaCGCTCTTTGATGCGCGACGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTATCTCGCACAGGGTCAGGGCTTCAATTGGCATTcc                                                    <  2:1627662/140‑1 (MQ=255)
                                                                                      aaaCGCTCTTTGATGCGCGAATCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACTTTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTcc                                                    <  2:152258/140‑1 (MQ=255)
                                                                                             ctTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCg                                             <  1:1632917/140‑1 (MQ=255)
                                                                                             ctTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCg                                             <  1:147659/140‑1 (MQ=255)
                                                                                                     cgcgAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCttt                                     <  1:1073097/140‑1 (MQ=255)
                                                                                                      gcgAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCtttt                                    >  1:1188309/1‑140 (MQ=255)
                                                                                                      gcgAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCttta                                    >  2:1568752/1‑139 (MQ=255)
                                                                                                       cgAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCTGTGGGTCAGGTATTAAATTTGCATTCCTTGTTCGATCTCtttta                                   <  1:949739/140‑2 (MQ=255)
                                                                                                       cgAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTc                                   <  1:318902/140‑1 (MQ=255)
                                                                                                           gCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGt                               <  2:285844/140‑1 (MQ=255)
                                                                                                                       aCAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGAGCTCTTTTCGCGTCTGGTTTAgcgc                   >  1:128937/1‑140 (MQ=255)
                                                                                                                        cAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCa                  <  1:575280/140‑1 (MQ=255)
                                                                                                                        cAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCa                  <  1:1415980/140‑1 (MQ=255)
                                                                                                                          aCGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTACAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATc                >  2:1387369/1‑140 (MQ=255)
                                                                                                                                        ttCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCa  >  2:1185700/1‑140 (MQ=255)
                                                                                                                                        ttCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCa  >  1:94671/1‑140 (MQ=255)
                                                                                                                                        ttCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCa  >  1:1244796/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                        
CACGCATATGCAGCCCCTGAATGGAGCGACCGGAACCTTTCAGCGTACCGTCACCTGCTGGCGCAAACGCGATAACCGGACGGTGAAAACGCTCTTTGATGCGCGAAGCCAGAATACCGACAACGCCCTGATGCCATTCGGGGTGATACATTGCCAGCCCGCCGGGTAGCGTGTCACGGCTGCGCTCCAGTTTCTCGCACAGGGTCAGGGCTTCAATTTGCATTCCTTGTTCGATCTCTTTTCGCGTCTGGTTTAGCGCATCGAGTTCATTTGCCA  >  NC_000913/3036879‑3037154

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