Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,236,156 T→C 5.5% intergenic (+8/‑106) yjbE → / → yjbF extracellular polysaccharide production threonine‑rich protein/extracellular polysaccharide production lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,236,1560TC 5.5% 240.4 / 3.0 72intergenic (+8/‑106)yjbE/yjbFextracellular polysaccharide production threonine‑rich protein/extracellular polysaccharide production lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (33/35);  new base C (2/2);  total (35/37)
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.01e-01

ACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAG  >  NC_000913/4236008‑4236303
                                                                                                                                                    |                                                                                                                                                   
aCGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtat                                                                                                                                                     <  2:279051/149‑1 (MQ=255)
aCGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtat                                                                                                                                                     <  4:254507/149‑1 (MQ=255)
  gTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtatgt                                                                                                                                                   >  4:508081/1‑149 (MQ=255)
   tCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtatgta                                                                                                                                                  >  1:30357/1‑149 (MQ=255)
    cGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAgtatgtat                                                                                                                                                 <  1:84051/149‑1 (MQ=255)
      gTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATcc                                                                                                                                               <  1:394402/149‑1 (MQ=255)
        ttCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATcccc                                                                                                                                             <  3:511886/149‑1 (MQ=255)
                ggAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAaataat                                                                                                                                     <  2:455457/149‑1 (MQ=255)
                ggAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAaataat                                                                                                                                     >  1:15549/1‑149 (MQ=255)
                  aGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTc                                                                                                                                   <  1:376478/149‑1 (MQ=255)
                    ttCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGa                                                                                                                                 <  3:251203/149‑1 (MQ=255)
                    ttCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGa                                                                                                                                 <  1:435304/149‑1 (MQ=255)
                    ttCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGa                                                                                                                                 >  4:662614/1‑149 (MQ=255)
                      cctccGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGt                                                                                                                               >  3:27090/1‑149 (MQ=255)
                       ctccGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTc                                                                                                                              >  2:394871/1‑149 (MQ=255)
                       ctccGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTc                                                                                                                              >  4:318861/1‑149 (MQ=255)
                        tccGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCa                                                                                                                             >  4:164709/1‑149 (MQ=255)
                          cgcgACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCAtt                                                                                                                           >  2:422406/1‑149 (MQ=255)
                           gcgACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTg                                                                                                                          >  1:690344/1‑149 (MQ=255)
                                  gCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCtgt                                                                                                                   <  3:508081/149‑1 (MQ=255)
                                    gTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTgg                                                                                                                 >  1:232619/1‑149 (MQ=255)
                                        gCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAg                                                                                                             >  3:203601/1‑149 (MQ=255)
                                         cACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGa                                                                                                            <  1:422406/149‑1 (MQ=255)
                                          aCCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGaa                                                                                                           <  2:690344/149‑1 (MQ=255)
                                          aCCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGaa                                                                                                           >  2:317295/1‑149 (MQ=255)
                                          aCCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGaa                                                                                                           <  2:15549/149‑1 (MQ=255)
                                          aCCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGaa                                                                                                           >  1:8251/1‑149 (MQ=255)
                                            cGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCtgt                                                                                                                   >  4:644624/1‑139 (MQ=255)
                                            cGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCtgt                                                                                                                   <  3:644624/139‑1 (MQ=255)
                                                 gCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGaa                                                                                                    <  4:31562/149‑1 (MQ=255)
                                                 gCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGaa                                                                                                    <  4:45549/149‑1 (MQ=255)
                                                     gtcggTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGa                                                                                                <  2:199576/147‑1 (MQ=255)
                                                     ggcggTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGa                                                                                                <  3:662614/149‑1 (MQ=255)
                                                         gtgggtgTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATtaa                                                                                            <  4:65609/149‑1 (MQ=255)
                                                          tgggtgTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATtaac                                                                                           <  4:27090/149‑1 (MQ=255)
                                                            ggtgTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATtaacca                                                                                         <  2:252019/149‑1 (MQ=255)
                                                            ggtgTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATtaacca                                                                                         <  3:318861/149‑1 (MQ=255)
                                                               gtCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATtaaccataa                                                                                      <  4:376984/149‑1 (MQ=255)
                                                                          ggCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCgg                                                                           >  1:422718/1‑149 (MQ=255)
                                                                          ggCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCgg                                                                           <  1:394871/149‑1 (MQ=255)
                                                                          ggCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTACGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCgg                                                                           >  4:126610/1‑149 (MQ=255)
                                                                          ggCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTACGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCgg                                                                           >  4:538334/1‑149 (MQ=255)
                                                                           gCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGgt                                                                          >  1:681037/1‑149 (MQ=255)
                                                                             aaCCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGgtgt                                                                        >  3:691622/1‑149 (MQ=255)
                                                                                 ggTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGtt                                                                    <  4:691622/149‑1 (MQ=255)
                                                                                   tGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACCGGAGTGTGGTTAt                                                                  <  2:8251/149‑1 (MQ=255)
                                                                                    ggCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTAtt                                                                 >  4:81472/1‑149 (MQ=255)
                                                                                     gCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACCGGAGTGTGGTTAtgg                                                                <  1:317295/149‑3 (MQ=255)
                                                                                           aTGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGcccc                                                          >  2:311353/1‑149 (MQ=255)
                                                                                             ggTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTCGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCctc                                                        >  1:119860/1‑149 (MQ=255)
                                                                                              gTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCctct                                                       >  3:88845/1‑149 (MQ=255)
                                                                                                 cTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGga                                                    >  1:101087/1‑149 (MQ=255)
                                                                                                 cTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGAATAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGga                                                    >  2:652204/1‑149 (MQ=255)
                                                                                                   aaTACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGgaga                                                  >  3:334869/1‑149 (MQ=255)
                                                                                                    aTACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGgaga                                                  >  3:100414/1‑148 (MQ=255)
                                                                                                    aTACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGgaga                                                  <  4:100414/148‑1 (MQ=255)
                                                                                                    aTACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGagaa                                                 <  1:311353/149‑1 (MQ=255)
                                                                                                          ggACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGga                                                    <  1:554739/140‑1 (MQ=255)
                                                                                                          ggACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGga                                                    >  2:554739/1‑140 (MQ=255)
                                                                                                              caccacAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGaa                                       <  2:30357/149‑1 (MQ=255)
                                                                                                                    aaCCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGAc                                  <  2:204593/148‑1 (MQ=255)
                                                                                                                    aaCCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGAc                                  >  1:204593/1‑148 (MQ=255)
                                                                                                                      ccACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTg                               >  4:231239/1‑149 (MQ=255)
                                                                                                                      ccACGACCACCAGTACCCAGTAATAAAGTATGTATCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATtaa                                                                                            <  2:601057/87‑1 (MQ=255)
                                                                                                                      ccACGACCACCAGTACCCAGTAATAAAGTATGTATCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATtaa                                                                                            >  1:601057/1‑87 (MQ=255)
                                                                                                                           accaccAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTc                          <  4:500988/149‑1 (MQ=255)
                                                                                                                             caccaGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCAt                        >  3:265801/1‑149 (MQ=255)
                                                                                                                              accaGTACCCAGTAATAAAGTACGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCAtt                       <  3:126610/149‑1 (MQ=255)
                                                                                                                              accaGTACCCAGTAATAAAGTACGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCAtt                       <  3:538334/149‑1 (MQ=255)
                                                                                                                                  gtccccAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGaa                                       >  2:136711/4‑129 (MQ=255)
                                                                                                                                  gTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGaa                                       <  1:136711/129‑1 (MQ=255)
                                                                                                                                  gTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTa                   <  2:119860/149‑1 (MQ=255)
                                                                                                                                  gTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTa                   <  4:334869/149‑1 (MQ=255)
                                                                                                                                      ccAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTg               >  2:515446/1‑149 (MQ=255)
                                                                                                                                         gTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCt            >  1:677348/1‑149 (MQ=255)
                                                                                                                                            ataaAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGtt         <  3:81472/149‑1 (MQ=255)
                                                                                                                                            ataaAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATAAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGtt         >  4:472915/1‑149 (MQ=255)
                                                                                                                                                  tatgtatCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACa   >  1:180586/1‑149 (MQ=255)
                                                                                                                                                   atgtatCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAg  <  4:482429/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGTCGGTTTCGGCGGGGAGTTCCTCCGCGACCAGCGTCAGCACCGTAAGCTCGGCGGTGGGTGTCGCGCTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAG  >  NC_000913/4236008‑4236303

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: