Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,274,740 G→A 7.1% intergenic (+79/+20) ssb → / ← yjcB single‑stranded DNA‑binding protein/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,274,7400GA 7.1% 215.7 / 3.7 68intergenic (+79/+20)ssb/yjcBsingle‑stranded DNA‑binding protein/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (34/29);  new base A (2/3);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GTGGTGGTCAGCCGCAGGGCGGTTGGGGTCAGCCTCAGCAGCCGCAGGGTGGCAATCAGTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACG  >  NC_000913/4274507‑4274988
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
ggagacagccgCCGCAGGGGGGGTGGGGGCAGCCTCAGCGTCGGCAGGGTGGCAATCAGTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCa                                                                                                                                                                                                                  <  1:62367/264‑1 (MQ=255)
       tcagcgcggggggggTTGGGGTCAGCTGAGGGGGCCGCAGGGTGGTAATCGGTGGAGCGGGGGCGCGCAGTCTGGCCCGAATCAGTCCGGTCCGGCAGCGCCGTCTAACGAGCCGACGATGGACTTTGATGATGACATTCCGTGCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCg                                                                                                                                                                                                          <  4:499355/266‑1 (MQ=255)
           ccGCAGGGCGGTTGGGGTCAGCCTCAGCAGCCGCAGGGTGGCAATCAGTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAgg                                                                                                                                                                                                                                  >  2:393775/1‑247 (MQ=255)
                   cGGTTGGGGTCAGCCTCAGCAGCCGCAGGGTGGCAATCAGTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAATCAGCTCCGGCAGCGCCGTCTAACGAGCCCCCGATGGACTTTGATGATTACATTCCGTTCTGATTNGTCATTAAAACATTAGGTTATATTGTTTGAAGGTGTCTGATTAAAGCATCTGCCAGCCATACAAAAGAAGCTCCNGTTATCGAGCCTTCTAAACCTGAGGACAAANNCATTGGCCACTGTGGCgaggcggtttc                                                                                                                                                                                                >  3:62367/1‑263 (MQ=255)
                                        gcCGCAGGGTGGCAATCAGTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCTTCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCTTCTCTTCTTTAATCCCCGACt                                                                                                                                                                         >  4:107526/1‑275 (MQ=255)
                                                     aaTCAGTTCAGCGGCGGCCCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTtgtg                                                                                                                                                               <  2:107527/272‑1 (MQ=255)
                                                          gTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATa                                                                                                                                                        <  2:109124/274‑1 (MQ=255)
                                                                                      cagcagTCCGCTCCTGCAGCTCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAg                                                                                                                           <  1:70027/275‑1 (MQ=255)
                                                                                         cagTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTAt                                                                                                                                                                                                                                                      <  5:360168/149‑1 (MQ=255)
                                                                                           gTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCAt                                                                                                                                                                                                                        >  3:95287/1‑177 (MQ=255)
                                                                                               gCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGc                                                                                                                                                                                                                    <  1:95287/177‑1 (MQ=255)
                                                                                                 tcaggCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCtt                                                                                                                                                                                                                                              <  8:345837/146‑1 (MQ=255)
                                                                                                 tCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCtt                                                                                                                                                                                                                                              >  5:110116/1‑149 (MQ=255)
                                                                                                 tCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGc                                                                                                                          >  1:37718/1‑265 (MQ=255)
                                                                                                              tCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGt                                                                                                                                                                                                                                 >  8:333744/1‑149 (MQ=255)
                                                                                                                 aaCGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCaa                                                                                                                                                                                                                              >  6:32622/1‑149 (MQ=255)
                                                                                                                    gAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAAt                                                                                                                                                                                                                           >  6:470459/1‑149 (MQ=255)
                                                                                                                      gccgccGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTAACATCACCGACTCCATGTTGTGCCAGATAGCGGTGATGTCCTCCATTGTAATGCCAAGCAACCCCATTGCCAACCACCAGGCGGAaaca                                                                                            >  4:495775/1‑274 (MQ=255)
                                                                                                                        cgccgATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCATATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAaaca                                                                                            >  2:481815/1‑272 (MQ=255)
                                                                                                                                        tgatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCTTc                                                                                                                                                                                                                                             >  7:467289/1‑111 (MQ=255)
                                                                                                                                        tgatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCATAACGGAGGCTTc                                                                                                                                                                                                                                             <  8:467289/111‑1 (MQ=255)
                                                                                                                                         gatgaCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCACTGTGGCGgtgt                                                                                                                                                                                                      >  5:235360/1‑149 (MQ=255)
                                                                                                                                             aCATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGa                                                                                                                                                                                                  <  5:8654/149‑1 (MQ=255)
                                                                                                                                               aTTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAtt                                                                                                                                                                                                <  5:470459/149‑1 (MQ=255)
                                                                                                                                                ttCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTc                                                                                                                                                                                               >  7:48841/1‑149 (MQ=255)
                                                                                                                                                   cGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtc                                                                                                                                                                                            <  7:333744/149‑1 (MQ=255)
                                                                                                                                                    gTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtct                                                                                                                                                                                           >  8:230628/1‑149 (MQ=255)
                                                                                                                                                     ttCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCAtctc                                                                                                                                                                                          <  6:201609/149‑1 (MQ=255)
                                                                                                                                                       cTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtctt                                                                                                                                                                                        <  6:235360/149‑1 (MQ=255)
                                                                                                                                                          aTTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTa                                                                                                                                                                                     <  6:110116/149‑1 (MQ=255)
                                                                                                                                                                 aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                                                                                                                              >  5:80005/1‑149 (MQ=255)
                                                                                                                                                                 aTTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAc                                                                                                                                                                              >  8:382649/1‑149 (MQ=255)
                                                                                                                                                                  ttAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCAcc                                                                                                                                                                             >  7:428162/1‑149 (MQ=255)
                                                                                                                                                                      aaCAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCg                                                                                                                                                                                                          >  5:141046/1‑116 (MQ=255)
                                                                                                                                                                      aaCAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCg                                                                                                                                                                                                          <  6:141046/116‑1 (MQ=255)
                                                                                                                                                                         aaTAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGaa                                        <  2:495780/275‑1 (MQ=255)
                                                                                                                                                                             ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttct                                                                                                                                                                                      <  7:370254/129‑1 (MQ=255)
                                                                                                                                                                             ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCtcttct                                                                                                                                                                                      >  8:370254/1‑129 (MQ=255)
                                                                                                                                                                             ggTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGAtgtt                                                                                                                                                                  <  8:48841/149‑1 (MQ=255)
                                                                                                                                                                                  tatTGTTTTAAGGTGGATGATTAAAGGATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGcc                                                                                                                                                             >  6:170626/1‑149 (MQ=255)
                                                                                                                                                                                   atTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCCCCCAAAGCCCGACATAAACACACCGACCAAAACGGAtttt                              >  4:403142/1‑275 (MQ=255)
                                                                                                                                                                                   atTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTAGTGCCa                                                                                                                                                            >  6:2413/1‑149 (MQ=255)
                                                                                                                                                                                       ttttAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAgg                          >  2:339434/1‑274 (MQ=255)
                                                                                                                                                                                         ttAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCgg                                                                                                 >  4:115404/1‑202 (MQ=255)
                                                                                                                                                                                           aaGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGAc                                                                                  <  2:26313/215‑1 (MQ=255)
                                                                                                                                                                                             ggTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAaac                                                                                             <  2:115405/202‑1 (MQ=255)
                                                                                                                                                                                             ggTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGa                                                          <  2:323020/237‑1 (MQ=255)
                                                                                                                                                                                                 gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGAAAGCCCGAAATAATCAGACCAACGAAAACGACGGTTGAGTGTGATTGCCa                >  1:14866/1‑275 (MQ=255)
                                                                                                                                                                                                 gatgatTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGAACACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCa                >  1:153684/1‑275 (MQ=255)
                                                                                                                                                                                                   tgatTAAAGGATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTcc                                                                                                                                            <  5:170626/149‑1 (MQ=255)
                                                                                                                                                                                                     atTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCaa                                                                                                        <  2:49500/183‑1 (MQ=255)
                                                                                                                                                                                                     atTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGc                                                              <  1:69762/225‑1 (MQ=255)
                                                                                                                                                                                                         aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAt                                                                                                                                                                                                 <  8:356006/90‑1 (MQ=255)
                                                                                                                                                                                                         aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGAt                                                                                                                                                                                                 >  7:356006/1‑90 (MQ=255)
                                                                                                                                                                                                         aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTcc                                                                                                                                            >  5:228866/1‑143 (MQ=255)
                                                                                                                                                                                                         aaGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTcc                                                                                                                                            <  6:228866/143‑1 (MQ=255)
                                                                                                                                                                                                             aTCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGc                                                                                                                                                              <  8:40428/121‑1 (MQ=255)
                                                                                                                                                                                                             aTCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGc                                                                                                                                                              >  7:40428/1‑121 (MQ=255)
                                                                                                                                                                                                              tCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                                                                                                                                 <  7:230628/149‑1 (MQ=255)
                                                                                                                                                                                                              tCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCCGACCGACGAAAACGGATGTTCAGTGTGCTGGCCAGAAACATCATTAc   >  2:51122/1‑275 (MQ=255)
                                                                                                                                                                                                                tGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGCCCGCTGCTGATCACCGCAAGCCCGACATAAGCAGACCGACGACAACGGATGTTGAGTCTGCTGGACAGAAACATCATTACg  >  1:137386/1‑274 (MQ=255)
                                                                                                                                                                                                                tGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACg  >  1:58252/1‑274 (MQ=255)
                                                                                                                                                                                                                         tAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAAcc                                                                                                                      >  7:265810/1‑149 (MQ=255)
                                                                                                                                                                                                                          aaaaaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCAtt                                                                                                                                      >  6:305464/1‑132 (MQ=255)
                                                                                                                                                                                                                          aaaaaaGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCAtt                                                                                                                                      <  5:305464/132‑1 (MQ=255)
                                                                                                                                                                                                                              aagaagCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCAtt                                                                                                                 <  5:32622/149‑1 (MQ=255)
                                                                                                                                                                                                                                 aagCCTCCATAACGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                                                                                                                                 <  7:189460/130‑1 (MQ=255)
                                                                                                                                                                                                                                 aagCCTCCATAACGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGGTGTGCCAGATAGCGGTGATGTCCGCCATTGTGAt                                                                                                                                 >  8:189460/1‑130 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
GTGGTGGTCAGCCGCAGGGCGGTTGGGGTCAGCCTCAGCAGCCGCAGGGTGGCAATCAGTTCAGCGGCGGCGCGCAGTCTCGCCCGCAGCAGTCCGCTCCGGCAGCGCCGTCTAACGAGCCGCCGATGGACTTTGATGATGACATTCCGTTCTGATTTGTCATTAAAACAATAGGTTATATTGTTTTAAGGTGGATGATTAAAGCATCTGCCAGCCATAAAAAAGAAGCCTCCGTTATGGAGGCTTCTACGTATCAGGTCAAAATCATTGGCCATTGTGGCGGTGTCTGATTCATCTCTTCTATCATCACCGACTCGATGTTGTGCCAGATAGCGGTGATGTCCGCCATTGTGATGCCAAGCAACCCCATTGCGAACCAACAGGCGGAAACAACGCCCAGACCGCTGCTGATCACCGCAAGCCCGACATAATCAGACCGACGAAAACGGATGTTGAGTGTGCTGGCCAGAAACATCATTACG  >  NC_000913/4274507‑4274988

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