Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,624,391 C→G A133A (GCG→GCC)  ytjB ← putative inner membrane protein Smp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,624,3910CG90.0% 161.6 / 3.3 60A133A (GCG→GCC) ytjBputative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base C (5/1);  new base G (32/22);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGC  >  NC_000913/4624143‑4624633
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
ggTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTTCTGCCAGCGGGTAAGTTTTCCCTGTAGCAGCGGGCGGGTTAGCACTACGCCGATTGCCAGTGAAAACAGCAACAACAGGCGTAAAATGTTAATGGTGTTATTCACCTGTTTGGCTTCGGTGGCGAGCGTATGGGGGTCGGGCGTCAGGGGCCGGTAGCCGAGCGGGCCGTTTTTTCCCgg                                                                                                                                                                                                                                                   >  1:78341/1‑250 (MQ=255)
         tttttCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCGGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTGAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTggtt                                                                                                                                                                                                                                          <  1:327081/250‑3 (MQ=255)
         tttttCCCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAATGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCTAGCGGTCCGTTTTTCCCGG‑CAATTggtt                                                                                                                                                                                                                                          <  2:387897/250‑3 (MQ=255)
                        tcctcttccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTAAGTTTACCCTGTAGCAGCGTGCGGGGCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGGGTCGAGCGTCAGGCGCAGATAGCCGAGCGGGCCGTTTTTCCCCG‑CAATTGGGTCGACAATCTGCTgggt                                                                                                                                                                                                                          >  1:397622/1‑249 (MQ=255)
                         cctcttccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCCGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCCACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGGCCGTTTTTCCCCG‑CAATTGGTTCGACAATCTGCTGGttt                                                                                                                                                                                                                         >  1:177086/1‑250 (MQ=255)
                           tcttccGGCACCGGTTTGCTGGCGGTTAACAGGCAGGGCGATTTCTGCCAGCGGGTACGTTTACCCTGTAGCAACGTGCGGGGCAGCACTACGCCGATTTCCAGTGAAAGCAGCAAAATCAGGCGTAAAATGTTTGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGGGTCGAGCGTCAGGCGCAGATAGCCGCGCGGTCCGTTTTTTCCGG‑CAATTTGGTCGACAAATTGCTGGttta                                                                                                                                                                                                                        >  1:8934/1‑248 (MQ=255)
                           tcttccGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAAACTGCTGGTTaa                                                                                                                                                                                                                        >  1:409397/1‑250 (MQ=255)
                                        gTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTTTTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCt                                                                                                                                                                                                          <  1:188445/251‑1 (MQ=255)
                                                  gggttaCAGGAATGGCGTTTGCAGCCAGCGGGTCCGTTTACCCTGTAGCAGCGTGGGGGTCAGCACTCCTCCGATTGCCAGGGAAAGCAGCAACATCAGGCGTTAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATTGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTTGACAATCTGCTTGTTAAAATAGCCGCCCGCTTTTTTACCGt                                                                                                                                                                                                <  2:238658/246‑1 (MQ=255)
                                                      aaCAGGAAGGGCGATTGCTGCCAGCGGGTACGTTTACCCGGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGa                                                                                                                                                                                             <  2:136551/250‑1 (MQ=255)
                                                              gggCGATTGCTGCCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGcc                                                                                                                                                                                    >  1:228406/1‑251 (MQ=255)
                                                                         gcCAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCa                                                                                                                                                                           <  1:292156/249‑1 (MQ=255)
                                                                          ccAGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACtt                                                                                                                                                                        >  1:446837/1‑251 (MQ=255)
                                                                           cAGCGGTTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTc                                                                                                                                                                       <  1:467263/251‑1 (MQ=255)
                                                                            aGCGGGTACGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGAAAAATCAGGCGTAAAATGTTCGTGGGGTTATCCACCTGTTGGGCTTCGGTGGCGCGCGTATGGGTGTCGAGCGTCAGGCGCCGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAAATGGTTCGCCAATCTGCTGGTTCAAATAGCCGCCCGCTTTTTTACCGTCGACCGCCAGACGGTCGCGCACTTCg                                                                                                                                                                      >  2:179921/1‑251 (MQ=255)
                                                                               gggTACGTTTACCCTGTAACAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGCGCGTATGGGTGTCGAGCGTCAGGCGCAGCTAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACACTCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACg                                                                                                                                                                   >  1:32713/1‑251 (MQ=255)
                                                                                   ccGTTTGCCCTGTAGCCGCGTGCGGGTCAGCACTGCGCGGATTGCAAGTGAAAGGAGCCACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTGGGTGGGGGGGGTATGGGTGGCGAGCGTCAGGCGCAGATAGCCTAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGTTTTTTTACCGTCGAGCGCCAGCCGGTCGCGAACTTCGACGCttt                                                                                                                                                               <  2:188910/250‑1 (MQ=255)
                                                                                     gTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAACTGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCCGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCTCTTTTTTCCCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCg                                                                                                                                                             >  2:288855/1‑251 (MQ=255)
                                                                                      tttACCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGCCGCTTTCGc                                                                                                                                                            >  2:467163/1‑251 (MQ=255)
                                                                                         aCCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTCCCgac                                                                                                                                                                                               >  2:124168/1‑211 (MQ=255)
                                                                                         aCCCTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCgac                                                                                                                                                                                               <  1:124168/213‑3 (MQ=255)
                                                                                            cTGTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATAAGGAGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCCCCAGAAc                                                                                                                                                      >  1:311513/1‑251 (MQ=255)
                                                                                              gTAGCAGCGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACgg                                                                                                                                                    <  1:172211/251‑2 (MQ=255)
                                                                                                agcagcGTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACgggc                                                                                                                                                  <  2:177086/251‑4 (MQ=255)
                                                                                                      gTGCGGGTCAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATa                                                                                                                                              >  1:104233/1‑249 (MQ=255)
                                                                                                              cAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTa                                                                                                                                                                                                                         <  2:62348/166‑1 (MQ=255)
                                                                                                              cAGCACTACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTa                                                                                                                                                                                                                         >  1:62348/1‑166 (MQ=255)
                                                                                                              cAGCACGACGCCGATTGCCAGTGAAAGCAGCAACATCAGGCTCAAAATGTTAGTGGTGTTATCCACCTGTTGGGATTCGGTGGCGCGCGTATGGGGGTCGAGCGCCAGGCGCAGATAGCCGATCGGTACGTTTTT‑CCCGACAATAGGTTAGACAATCAGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGGGCGCCAGCCGGTAGCGCACTTCGACGCTTTCACCAGAACGGGCGATACGATCGccgg                                                                                                                                    >  2:96660/1‑249 (MQ=255)
                                                                                                                           aTTGCCAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCCTacac                                                                                                                       >  2:51677/1‑251 (MQ=255)
                                                                                                                               ccAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGCACGGGCGATAAGATCGCCTTGTTCGTCATACACGccc                                                                                                                   >  1:227987/1‑251 (MQ=255)
                                                                                                                               ccAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATacccgccc                                                                                                                   >  1:456749/1‑245 (MQ=255)
                                                                                                                               ccAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGGGGTGTTATCCACCTGTTGGGCTTAGGTGGCGAACGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTAGACAATCTGATGGTTAAAATAGCCGCCCGCTTTTTTACCGTCATGCGCCAGCCGGTCGCGCACTTAGACGCTTTCTCCAGAACGGGCGATGAGATCGCCTTGTTCGCCATacacgcac                                                                                                                   >  2:346384/1‑247 (MQ=255)
                                                                                                                                cAGTGAAAGCAGCAACATGAGGTGTAAAATGTTAGTGGTGTTATCCACCGGTTGGGCTTCGGTGGCGATCGTATGGGTTTCGAGCGGCAGGCGCAGATAGTCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCcg                                                                                                                  <  2:399119/251‑1 (MQ=255)
                                                                                                                                cAGTGAAAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCcg                                                                                                                  <  1:48515/251‑1 (MQ=255)
                                                                                                                                        gcagcaAAATCAGGGGGAAATTGTTTGTGGGTTTTTCCCCCTGTTGGGCTTCGGTGGCGGGTGGATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTTGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTGACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTGCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGGCGAgg                                                                                                          <  2:67870/251‑1 (MQ=255)
                                                                                                                                             aaCATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGTGCTTCTGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGCGCGGTCCGTTTTTCCCGG‑CACTTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCCCTTCGACGCTTTCGCCAGAACGGGCGTTAAGATCGCCTTGTTCGTCATACCCGCCCCCGTCggagttacc                                                                                                     >  2:146651/1‑243 (MQ=255)
                                                                                                                                                      cGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTAGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCACGCGTCGAGGATACGGATTTCATc                                                                                            >  1:298187/1‑251 (MQ=255)
                                                                                                                                                        tAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGCTACGGCTTTCATCCg                                                                                          >  1:479738/1‑251 (MQ=255)
                                                                                                                                                        tAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCg                                                                                          >  1:80668/1‑251 (MQ=255)
                                                                                                                                                           aaTGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCCGGCGCAGCTAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACACTCTGCTGGTTAACATAGACGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGCTCGCCTTCTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTTATTCGttt                                                                                       >  1:61068/1‑250 (MQ=255)
                                                                                                                                                               ttAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGCGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGCCCGGGCGATAAGATCGCCTTGTTCGTCCTACACGACCGCCTCGAGGGTAGGGCtttacttcgttaaatg                                                                                   >  2:300441/1‑238 (MQ=255)
                                                                                                                                                                 aGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGAGCACTTCGCCGCTTTCGCCAGAACGGGCGATCAGATCGCCTTGATCGTCATACACGCCCGCGTCGAGGAAACGGCTTTCACCCGCTAACTgct                                                                                 >  2:145508/1‑249 (MQ=255)
                                                                                                                                                                        ttATCCACCTGTTGGGCTCCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGa                                                                           <  1:126835/250‑1 (MQ=255)
                                                                                                                                                                         tATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTCCCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGCACGGGCGATAAGATCGCCTTGTTCGTCATACAAGCCCGCGTCGAGGATACGCCTTTCATCCGTTagatgttcgggggtc                                                                         >  2:350471/1‑237 (MQ=255)
                                                                                                                                                                          aTCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATc                                                                         <  1:386374/250‑1 (MQ=255)
                                                                                                                                                                          aTCCACCGGGTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATc                                                                         <  1:146651/250‑1 (MQ=255)
                                                                                                                                                                              aCCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCTGCCGGTCGCGCACTTCGACGATTTCGCCAGAACGGGCGATAAGATCGCCTTGGTCGTCATACACGCCCGCGTCGAGGGTACGACTTTCATCCGTTAACTgttagtggctcgcctg                                                                    >  2:29241/1‑236 (MQ=255)
                                                                                                                                                                                              gggCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCTCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATcc                                                    <  1:315481/250‑1 (MQ=255)
                                                                                                                                                                                               ggCGAGCGTATGGGTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCg                                                   >  1:387025/1‑251 (MQ=255)
                                                                                                                                                                                                          gggTGTCGAGCGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGATCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTa                                        <  1:349370/251‑1 (MQ=255)
                                                                                                                                                                                                                    cgtCAGGCGCAGATAGCCGAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATACGATCGCCTTGTTCGTCATACACGCCCGCGCCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCg                              >  1:59572/2‑251 (MQ=255)
                                                                                                                                                                                                                        aGGCGCAGATAGCCGAGCTGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGc                          <  1:168960/251‑1 (MQ=255)
                                                                                                                                                                                                                              aGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCTTCGAGGATACGGCTTTCATCCGTTAACTGATCg                                                                               >  1:403405/1‑192 (MQ=255)
                                                                                                                                                                                                                                      gAGCGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGtta            <  1:350471/251‑3 (MQ=255)
                                                                                                                                                                                                                                         cGGTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGTTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCt         <  1:482196/251‑1 (MQ=255)
                                                                                                                                                                                                                                          ggTCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTg        >  1:499212/1‑251 (MQ=255)
                                                                                                                                                                                                                                            tCCGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTTACCTgcc      >  2:315476/1‑249 (MQ=255)
                                                                                                                                                                                                                                             ccGTTTTTCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGCGAGTCGGTACGCATCCGCGGTGCAACGTTCAGCGTTACCTGAcc     >  2:190731/1‑250 (MQ=255)
                                                                                                                                                                                                                                                tttttCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCAACCGGCGAGGTGGTCCGCATGAGCGGAGCacagttcaacgttaactgaggcgc  >  2:328490/1‑228 (MQ=255)
                                                                                                                                                                                                                                                tttttCCCGG‑CAATTGGTTCGACAATCTGCTGGTTAAAATAGCCGCCCGCTTTTTTACCGTCGAGCGCCAGACGGTCGCGCACTTCGACGCTTTCGCCAGAACGGGCGATAAGATCGCCTTGTTCGTCATACACGCCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATACGGCGAGTCGGAACGCATCAGCGGTGAAACGTTCAGCGTTACCTGACgcgc  >  1:25757/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
GGTCACTCTTTTTTCTCGCTTTCTTCCTCTTCCGGCACCGGTTTGCTGGCGGTTAACAGGAAGGGCGATTGCTGCCAGCGGGTGCGTTTACCCTGTAGCAGCGTGCGGGTCAGCACTACACCGATTGCCAGTGAGAGCAGCAACATCAGGCGTAAAATGTTAGTGGTGTTATCCACCTGTTGGGCTTCGGTGGCGAGCGTATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCG‑CAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAGTCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGC  >  NC_000913/4624143‑4624633

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑