Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 681,694 A→G 8.0% intergenic (+31/+29) djlC → / ← hscC J domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913681,6940AG 8.0% 226.6 / 8.9 77intergenic (+31/+29)djlC/hscCJ domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base G (4/2);  total (43/34)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-03

AGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAGCGCCCGCGCCCACAACTCTTCCAGTTTGGCTTTAAAGGTGCGATTGATGAGCATATCGCGCGGATAAATTTTCAATGCGGATAACCGCGTCCGACTCTCTTCAATCTGCTG  >  NC_000913/681430‑681965
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                 
aGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGAGGTCGATGAGCCAAAACCTGAAAATCATGCAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAAAGGTCCTGGTGATAGGTTATGCTCCTGGCAAAATCGCGAGTTTGTTTAGATGACGACGTGTTTTTTGCGCCAATGCGGGATCAAGCGTTTCcgc                                                                                                                                                                                                                                                                         >  4:207314/1‑275 (MQ=255)
 ggCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGCGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgcc                                                                                                                                                                                                                                                                        >  2:235969/1‑275 (MQ=255)
                          tttggaggcagatGTCAGCCGGAGGCTTTTGACGAAAGTGAAACTCAGGCCGAAGGATGAAGATGTAGATGCGCCAATACATGAAAATCATTAAGAAAAAAATGATGAGGGTGAAAAACCACATAGCATTACCAGAATGATCCAAATCAGTTTAACGGTCCTGGTGATATGTTATGCTCTCGGCAAAATCGCGATTTTGTTTAGCTGACTATGTGTTCTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGACAGCCGTGTTTATggg                                                                                                                                                                                                                                                <  4:12313/262‑1 (MQ=255)
                          aTTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                                                                                                                                                                                                <  2:307596/274‑1 (MQ=255)
                                   cATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACCCCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGTTGCAGCGTTTACGCCGCATCCGGTAGCCGTTTTTATGGGATTCCAATg                                                                                                                                                                                                                                       >  2:99794/1‑274 (MQ=255)
                                     ttGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGc                                                                                                                                                                                                                                     >  2:189499/1‑274 (MQ=255)
                                      tGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCt                                                                                                                                                                                                                                    >  2:349572/1‑274 (MQ=255)
                                                                           ccGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAAGGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcg                                                                                                                                                                                                <  2:245964/273‑1 (MQ=255)
                                                                               aggaTGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtct                                                                                                                                                                                                                                            >  2:210980/1‑225 (MQ=255)
                                                                                   tGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAACAAAATTACGTGGGTGAAAAACCACATAGCATTCCCAGCATTATCAACATCTGTTTAACTGTACGGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGGTTTTTAGCTGACGATGTGTTTTTTTCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                                                                                                                                                                                                        <  4:210977/225‑1 (MQ=255)
                                                                                             gATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATcc                                                                                                                                                                                                   >  2:23697/1‑252 (MQ=255)
                                                                                             gATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTg                                                                                                                                                                            >  2:308822/1‑275 (MQ=255)
                                                                                                 aGCCAAAACCTGAAAAACAGGAAAAAAAAAATGATGGGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGTTGCTGGTGTTAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcgc                                                                                                                                                                                               <  4:23697/252‑1 (MQ=255)
                                                                                                             aaaaTCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtct                                                                                                                                                                                                                                            >  1:58921/1‑195 (MQ=255)
                                                                                                                 tCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                                                                                                                                                                                                        <  3:58921/195‑1 (MQ=255)
                                                                                                                    tGAAGAAAAAAATGATGAGGGTGAAAAACCACATGGCATTACCAGCATTGTCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGa                                                                                                                                                     <  2:396575/275‑1 (MQ=255)
                                                                                                                    tGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAAATTCATCCACGCGGGTCATATCGTTTGACTGCAACGCCGCATCAAAATCCGTg                                                                                                                                                      >  2:56096/1‑274 (MQ=255)
                                                                                                                    tGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAGCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCNNNTGACTGCAACGCCGCaaacaaaaaagta                                                                                                                                                      >  3:152666/1‑262 (MQ=255)
                                                                                                                        aaaaaaaaTGAGGAGGGAGCAAAACCCAAGCTCCTAACCATCATGATAAAAATACGCTTAACGGTGCGGGTAAGAGGTTTTTCTGTCGGCACAATCGCGATGTTGTTTAGTTGACGATGTGTTTCTTGCGCCGATGCCGGATGCAGCTTTTAAGCCGCATTCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCCCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCAc                                                                                                                                                 <  4:189497/274‑1 (MQ=255)
                                                                                                                            aaaaTGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAc                                                                                                                                                                                                                                                                           <  7:212445/149‑1 (MQ=255)
                                                                                                                                 gatgaGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAAcgctgc                                                                                                                                                                                                                                                                      <  8:90461/149‑1 (MQ=255)
                                                                                                                                       ggTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcg                                                                                                                                                                                                                                                                          >  5:347684/1‑139 (MQ=255)
                                                                                                                                       ggTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcg                                                                                                                                                                                                                                                                          <  6:347684/139‑1 (MQ=255)
                                                                                                                                          gAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAg                                                                                                                                                                                                                                                             >  5:369884/1‑149 (MQ=255)
                                                                                                                                                           tACCAGCAGTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCCCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCTTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAgcgccc                                                                                                             <  2:472778/275‑1 (MQ=255)
                                                                                                                                                                     aTCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATcc                                                                                                                                                                                                   >  1:39399/1‑180 (MQ=255)
                                                                                                                                                                          aaTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaataat                                                                                                                                                                                                                             >  7:129114/1‑149 (MQ=255)
                                                                                                                                                                           aTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtc                                                                                                                                                                                                                            >  6:393967/1‑149 (MQ=255)
                                                                                                                                                                           aTCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTAGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCTGGATGTCGCGTTTACGCCGCATCCGGCAGCCgt                                                                                                                                                                                                                                                         >  4:363007/1‑120 (MQ=255)
                                                                                                                                                                               gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGttt                                                                                                                                                                                                                                                     <  2:363012/120‑1 (MQ=255)
                                                                                                                                                                                tttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATcc                                                                                                                                                                                                                                                                  >  7:30986/1‑106 (MQ=38)
                                                                                                                                                                                tttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATcc                                                                                                                                                                                                                                                                  <  8:30986/106‑1 (MQ=38)
                                                                                                                                                                                  tAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCACCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAGCGCCCGCGCCCCCAACTCTTCCAGTTTg                                                                                     >  2:478261/1‑275 (MQ=255)
                                                                                                                                                                                          cTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcgc                                                                                                                                                                                               >  2:445935/1‑163 (MQ=255)
                                                                                                                                                                                              ggATAGGTTATGCTTTCGGCAAAACCGCGATGTTGTTTAGCTGACGATGTGTGTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGc                                                                                                                                                                                           <  4:445930/162‑1 (MQ=255)
                                                                                                                                                                                               taTAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGAGGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAGCGCCCGCGCCCACAACTCTTCCAGTTTGGCTTTAAAgg                                                                           <  2:337288/271‑1 (MQ=255)
                                                                                                                                                                                               gATAGGTTATGCTCTCGGCAATATCGCGATGTTGTTTAGCTGACGATGTGTTGTTTGCGTCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAGCGCCCGCGCCCACAACTCTTCCAGTTTGGCTTTAAAgg                                                                           <  2:270280/272‑1 (MQ=255)
                                                                                                                                                                                               gATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACtt                                                                                                                                                                                                        >  8:243586/1‑149 (MQ=255)
                                                                                                                                                                                               gATAGGTTATGCTCTCGGCAAAAGCGCGATGTTGTTTAGCGGACGATGTGTTTTTTGCGCCGAGGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAGCGCCCGCGCCCACAACTCTTCCAGTTTGGCTTTAAAgg                                                                           <  1:140191/272‑1 (MQ=255)
                                                                                                                                                                                                  aGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTTTTTATGGGATCTCAATGGCTAAATACTCGCTCGCCCGCCGCCGACCTTCATCAACGCGGGCCATATCGTTTGACTGCAACGCCGCAACACAATCGGTGATCACCCGGCCGATCACTTCCAGCCCGTACACACGCCCCCGCGCCACACAACTCTCCAGTTTGGCTTTAAAGGttcct                                                                      >  4:115398/1‑270 (MQ=255)
                                                                                                                                                                                                  aGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAc                                                                                                                                                                                                 >  2:282766/1‑153 (MQ=255)
                                                                                                                                                                                                  aGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCACGCCCGCCGCCGAACTTCATCCACGCGGGCCCTATCGTTTGCCGTCAACTCCGCATCAAAATCGGTGAACACCCGGCCGCTCCCCTCCCGCCCCGCCCCCAGCGCCCGCCCCCACAACTCTTCCCGTTTTGCCTTAACGGtccta                                                                      >  4:326357/1‑270 (MQ=255)
                                                                                                                                                                                                   ggTTATGCTCTCGGCAAAATGGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAGCGCCCGCGCCCACAACTCTTCCAGTTTGGCTTTAAAGGTGCGAtt                                                                    <  2:207317/275‑1 (MQ=255)
                                                                                                                                                                                                   ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATc                                                                                                                                                                                                    >  8:6208/1‑149 (MQ=255)
                                                                                                                                                                                                      tATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCgg                                                                                                                                                                                             <  4:282762/153‑1 (MQ=255)
                                                                                                                                                                                                       aTGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAcg                                                                                                                                                                                                >  7:422225/1‑149 (MQ=255)
                                                                                                                                                                                                       aTGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCCAAATCGGTTATCACCCTGCCGATCCCTTCTCGCCAGTCACCCAACGCCCGCGCCCCCACCTCTTCCAGTTTGGCTTTAGATGTACGACtgatg                                                                >  4:172833/1‑275 (MQ=255)
                                                                                                                                                                                                        tGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGa                                                                                                                                                                                                                                               <  6:206397/101‑1 (MQ=255)
                                                                                                                                                                                                        tGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGa                                                                                                                                                                                                                                               >  5:206397/1‑101 (MQ=255)
                                                                                                                                                                                                           tctcGGCAGAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg                                                                                                                                                                                            >  7:437950/1‑149 (MQ=255)
                                                                                                                                                                                                           tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg                                                                                                                                                                                            >  5:253914/1‑149 (MQ=255)
                                                                                                                                                                                                           tctcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg                                                                                                                                                                                            >  6:97140/1‑149 (MQ=255)
                                                                                                                                                                                                             tcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCa                                                                                                                                                               >  2:409326/1‑176 (MQ=255)
                                                                                                                                                                                                              cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACt                                                                                                                                                                                                         >  3:54687/1‑133 (MQ=255)
                                                                                                                                                                                                              cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAGCGCCCGCGCCCACAACTCTTCCAGTTTGGCTTTAAAGGTGCGATTGATGAGCatat                                                         >  2:436561/1‑275 (MQ=255)
                                                                                                                                                                                                               ggCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAGCGCCCGCGCCCACAACTCTTCCAGTTTGGCTTTAAAGGTGCGATCGATGAGCATATc                                                        >  2:192515/1‑275 (MQ=255)
                                                                                                                                                                                                                 cAAAAGCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAAt                                                                                                                                                           <  4:409321/176‑1 (MQ=255)
                                                                                                                                                                                                                  aaaaTCGCGATGTTGTTTAGCGGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCAt                                                                                                                                                                                                     <  1:54687/133‑1 (MQ=255)
                                                                                                                                                                                                                                 tttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAAcgc                                                                                                                                                                      >  5:185437/1‑149 (MQ=255)
                                                                                                                                                                                                                                       tGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATc                                                                                                                                                                <  5:438638/149‑1 (MQ=255)
                                                                                                                                                                                                                                        gACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAGCGCCCGCGCCCACAACTCTTCCAGTTTGGCTTTAAAGGTGCGATTGATGAGCCTATCGCGCGGATAAATTTTCAATGCGGa                                >  2:5065/1‑274 (MQ=255)
                                                                                                                                                                                                                                                gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGt                                                                                                                                                       <  8:129114/149‑1 (MQ=255)
                                                                                                                                                                                                                                                gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGt                                                                                                                                                       <  5:393967/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   ttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGAt                                                                                                                                                    <  6:253914/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCa                                                                                                                                                  <  7:6208/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       gcgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCAcc                                                                                                                                                >  6:284479/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         gcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGACCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATTACCTGGCCGATCTCTTCTCGCCCGTCACTCAGCGCCCGCGCCCACAACTCTTACAGTTTTGCTTTAAAGGTTCGATTGATGCGCATACCCTGCTGACAAATTTTCACTTCGGATAACCGCGTCCGGttctc              >  4:159673/1‑275 (MQ=255)
                                                                                                                                                                                                                                                           ggctGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCCCGCTCGTCACCCAGCGCCCGCGCCCACAACTCTTTCAGTTTGGCTTTAAAGGTGCGATTGCTGGGCATACCGCGCGGATAAATTTTCAATGCGGATAACCGCGTTCGActctct             >  2:57065/4‑274 (MQ=255)
                                                                                                                                                                                                                                                              tGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCAAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGa                                                                                                                                         <  6:369884/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               gCCGGATGGTGCGGTTACTCAGCAGCCGGCAGCCGGTTGTATGGGAGCTCAATGGCTACATAGTCGCCCGCCAGCCGCCGAACGTAATCCAGGCGAGCCATATCGTTTGACTGGAACGCAGAAGCAAAATCGGTGATCACCCTGCCGATCTCTTCTCGCTCGTCTCCCAGCGCCCGCGCCCACACCTCTTCCAGTTTGGCTTTAAAGGTGCGATTGATGAGCATATCGCGCGGATAAATTTTCAATGCGGATAACCGCGTCCGACTCTCTTCAAt        <  4:235965/275‑1 (MQ=255)
                                                                                                                                                                                                                                                               gCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAt                                                                                                                                        <  5:97140/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                ccGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGAtc                                                                                                                                       <  5:284479/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   gATGCGGCGTACACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATctct                                                                                                                                    >  7:279855/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   gATGCGGCGTAAACGCTGTATCCGGCATCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATctct                                                                                                                                    >  7:183480/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   gATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATctct                                                                                                                                    >  6:13310/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     tGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCtcttc                                                                                                                                  <  8:422225/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gCAGCGTTTCAGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAGCGCCCGCGCCCACAACTCTTCCAGTTTGGCTTTAAAGGTGCGATTGATGAGCATATCGCGCGGATAAATTTTCAATGCGGATAACCGCGTCCGACTCTCTTCAATctgctg  <  2:115399/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                 
AGGCTGGAAGGCGGAGATGCTCATCCATTAGAGGCCATTGTCACCCGGATGCTTTTGACGAAAGTGAAACTCGAGCCGGAGGATGAAGATGTCGATGAGCCAACACCTGAAAATCATGAAGAAAAAAATGATGAGGGTGAAAAACCACAGAGCATTACCAGCATTATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTGACTGCAACGCCGCATCAAAATCGGTGATCACCCGGCCGATCTCTTCTCGCTCGTCACCCAGCGCCCGCGCCCACAACTCTTCCAGTTTGGCTTTAAAGGTGCGATTGATGAGCATATCGCGCGGATAAATTTTCAATGCGGATAACCGCGTCCGACTCTCTTCAATCTGCTG  >  NC_000913/681430‑681965

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