Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,982 G→C 7.6% intergenic (+28/+32) yhcG → / ← yhcH DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily/DUF386 family protein, cupin superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,368,9820GC 7.6% 234.7 / 9.7 66intergenic (+28/+32)yhcG/yhcHDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily/DUF386 family protein, cupin superfamily
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base C (3/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCGG  >  NC_000913/3368836‑3369130
                                                                                                                                                  |                                                                                                                                                    
gTTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtg                                                                                                                                                    >  4:784248/1‑149 (MQ=255)
 tttGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAgtgt                                                                                                                                                   >  3:443459/1‑149 (MQ=255)
      cTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCa                                                                                                                                              <  4:384416/149‑1 (MQ=255)
        aaCACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGa                                                                                                                                            >  4:368083/1‑149 (MQ=255)
          cacCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAAt                                                                                                                                          >  1:153929/1‑149 (MQ=255)
           acCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATc                                                                                                                                         >  1:643562/1‑149 (MQ=255)
              gTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCttt                                                                                                                                      <  1:793199/149‑1 (MQ=255)
                     cAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGa                                                                                                                               <  4:443459/149‑1 (MQ=255)
                      aaGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGaa                                                                                                                              >  3:44996/1‑149 (MQ=255)
                            atatAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGt                                                                                                                        >  4:68283/1‑149 (MQ=255)
                             tataAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGtt                                                                                                                       <  1:563246/149‑1 (MQ=255)
                                  gATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAg                                                                                                                  >  1:809096/1‑149 (MQ=255)
                                  gATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAg                                                                                                                  >  2:628814/1‑149 (MQ=255)
                                   aTGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCttt                                                                                                                                      >  3:762196/1‑128 (MQ=255)
                                   aTGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCttt                                                                                                                                      <  4:762196/128‑1 (MQ=255)
                                   aTGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGc                                                                                                                 <  2:438707/149‑1 (MQ=255)
                                   aTGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGc                                                                                                                 >  3:417281/1‑149 (MQ=255)
                                       aaCTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAAt                                                                                                                                          >  2:384965/1‑120 (MQ=255)
                                       aaCTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAAt                                                                                                                                          <  1:384965/120‑1 (MQ=255)
                                       aaCTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCAtt                                                                                                             >  2:184994/1‑149 (MQ=255)
                                                atgaGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGcc                                                                                                    >  1:579650/1‑149 (MQ=255)
                                                    gagaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTaa                                                                                                >  3:542429/4‑149 (MQ=255)
                                                    gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTaa                                                                                                >  1:385013/1‑149 (MQ=255)
                                                    gaAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTaa                                                                                                >  1:790829/1‑149 (MQ=255)
                                                     aaaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                               <  1:115767/149‑1 (MQ=255)
                                                     aaaCGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAc                                                                                               <  3:88419/149‑1 (MQ=255)
                                                          aCTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTca                                                                                          <  2:153929/149‑1 (MQ=255)
                                                            tCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcacc                                                                                        >  2:404899/1‑149 (MQ=255)
                                                               cAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACcaca                                                                                     >  3:635567/1‑149 (MQ=255)
                                                                  aTGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAAcc                                                                                  >  2:49934/1‑149 (MQ=255)
                                                                      aCTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCtttt                                                                              <  1:404899/149‑1 (MQ=255)
                                                                       cTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCttttt                                                                             <  4:417281/149‑1 (MQ=255)
                                                                          gTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAt                                                                          <  2:643562/149‑1 (MQ=255)
                                                                                 cacaGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCa                                                                   <  3:784248/149‑1 (MQ=255)
                                                                                               gagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcac                                                                                         >  1:526832/1‑113 (MQ=255)
                                                                                               gagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTcac                                                                                         <  2:526832/113‑1 (MQ=255)
                                                                                               gagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAAc                                                     <  1:163386/149‑1 (MQ=255)
                                                                                               gagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAAc                                                     <  2:790829/149‑1 (MQ=255)
                                                                                               gagGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAAc                                                     <  3:368083/149‑1 (MQ=255)
                                                                                                 ggAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGc                                                   >  4:772209/1‑149 (MQ=255)
                                                                                                        tataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTATATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTgg                                            >  4:22101/1‑149 (MQ=255)
                                                                                                        tataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTgg                                            <  1:184994/149‑1 (MQ=255)
                                                                                                         ataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtctc                                                                       >  1:349581/1‑121 (MQ=255)
                                                                                                         ataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAAtctc                                                                       <  2:349581/121‑1 (MQ=255)
                                                                                                         ataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATTCTg                                             <  4:801599/147‑1 (MQ=255)
                                                                                                         ataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATTCTg                                             >  3:801599/1‑147 (MQ=255)
                                                                                                         ataGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGt                                           <  4:635567/149‑1 (MQ=255)
                                                                                                            gACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTa                                        <  1:358656/149‑1 (MQ=255)
                                                                                                               gCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAcc                                                                                              <  1:527674/92‑1 (MQ=255)
                                                                                                               gCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAAcc                                                                                              >  2:527674/1‑92 (MQ=255)
                                                                                                               gCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGc                                     >  4:732853/1‑149 (MQ=255)
                                                                                                                 cGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCgg                                   <  1:765349/149‑1 (MQ=255)
                                                                                                                 cGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCgg                                   >  4:60367/1‑149 (MQ=255)
                                                                                                                  gTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGt                                  <  1:628814/149‑1 (MQ=255)
                                                                                                                   ttAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGtt                                 >  1:260922/1‑149 (MQ=255)
                                                                                                                        ttCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCt                            <  3:772209/149‑1 (MQ=255)
                                                                                                                        ttCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCt                            <  2:579650/149‑1 (MQ=255)
                                                                                                                           aGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGc                         >  3:95803/1‑149 (MQ=255)
                                                                                                                                   aagaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTaa                                                                           <  2:729518/88‑1 (MQ=255)
                                                                                                                                   aagaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTaa                                                                           >  1:729518/4‑91 (MQ=255)
                                                                                                                                     gaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGc                                                            <  1:658812/103‑1 (MQ=255)
                                                                                                                                     gaTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGc                                                            >  2:658812/2‑104 (MQ=255)
                                                                                                                                      aTTCTGCAGACACTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGc              >  2:290151/1‑149 (MQ=255)
                                                                                                                                              gACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCt                            <  1:283056/127‑1 (MQ=255)
                                                                                                                                              gACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCt                            >  2:283056/1‑127 (MQ=255)
                                                                                                                                              gACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTc      >  2:110651/1‑149 (MQ=255)
                                                                                                                                               aCAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTcc     >  3:686000/1‑149 (MQ=255)
                                                                                                                                               aCAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTcc     <  1:49934/149‑1 (MQ=255)
                                                                                                                                                  gtgtCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCgg  >  1:479772/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTTTGCCTAACACCGTTCTGGCAAGCGAATATAAGATGCAACTACCTGATGAGAAACGACTCGCAGATGAACTCGTTCGAACACAGGCGGTGCTAGAGGAAGGCTATAGACGCCGTTAATTTCAGCGGTTTCGAATTCTGCAGACAGTGTCTGCAGAATCTTTAAAACGAACTAAGTTTTAAGCCATTAAATCAGCCTTAACCTTCACCACAACCTTTTTAATCTCTCCAGGCTCGCCGACAACGCATCCTGGTTTATGCGGTTCACCTGGCATAAACACGGCGAACATTCCCGG  >  NC_000913/3368836‑3369130

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