Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,635,457 T→C 5.1% intergenic (‑147/‑64) rob ← / → creA right oriC‑binding transcriptional activator, AraC family/putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,635,4570TC 5.1% 237.6 / 2.8 78intergenic (‑147/‑64)rob/creAright oriC‑binding transcriptional activator, AraC family/putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base C (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CTCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGA  >  NC_000913/4635319‑4635603
                                                                                                                                          |                                                                                                                                                  
cTCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTa                                                                                                                                                   <  6:182530/140‑1 (MQ=255)
cTCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTa                                                                                                                                                   >  5:182530/1‑140 (MQ=255)
cTCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAatat                                                                                                                                           <  3:168916/148‑1 (MQ=255)
      cTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTAtt                                                                                                                                     >  5:71514/1‑148 (MQ=255)
        ttCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTAc                                                                                                                                  >  4:60088/1‑149 (MQ=255)
             caacGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGa                                                                                                                             >  7:212452/1‑149 (MQ=255)
                cGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACg                                                                                                                           >  3:209824/1‑148 (MQ=255)
                cGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGa                                                                                                                          <  1:122736/149‑1 (MQ=255)
                 gAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAatat                                                                                                                                           >  3:202143/1‑131 (MQ=255)
                 gAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAatat                                                                                                                                           <  4:202143/131‑1 (MQ=255)
                       cTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAAtt                                                                                                                   >  4:110372/1‑149 (MQ=255)
                           cATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTc                                                                                                               <  5:198771/149‑1 (MQ=255)
                           cATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTc                                                                                                               <  2:103923/149‑1 (MQ=255)
                            aTCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTc                                                                                                               <  6:46041/148‑1 (MQ=255)
                            aTCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTc                                                                                                               >  5:46041/1‑148 (MQ=255)
                            aTCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCg                                                                                                              >  5:201607/1‑149 (MQ=255)
                             tCAGGTAATTGGATAATATCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCgt                                                                                                             <  5:102331/149‑1 (MQ=255)
                             tCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAAt                                                                                                                                                 <  1:102498/113‑1 (MQ=255)
                             tCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAAt                                                                                                                                                 >  2:102498/1‑113 (MQ=255)
                                ggTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTCAATCAATATg                                                                                                                                          >  2:134090/1‑117 (MQ=255)
                                ggTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTCAATCAATATg                                                                                                                                          >  2:69891/1‑117 (MQ=255)
                                ggTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTCAATCAATATg                                                                                                                                          <  1:69891/117‑1 (MQ=255)
                                ggTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTCAATCAATATg                                                                                                                                          <  1:134090/117‑1 (MQ=255)
                                 gTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTaa                                                                                                         <  4:8144/149‑1 (MQ=255)
                                 gTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTATTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTaa                                                                                                         >  5:210785/1‑149 (MQ=255)
                                  tAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTaa                                                                                                         <  3:60088/148‑1 (MQ=255)
                                  tAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTaa                                                                                                         <  5:190275/148‑1 (MQ=255)
                                  tAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAAt                                                                                                        <  2:30343/149‑1 (MQ=255)
                                            ataGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATaaaaa                                                                                               <  6:201607/148‑1 (MQ=255)
                                              aGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATgg                                                                                            >  7:110940/1‑149 (MQ=255)
                                               gCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATgg                                                                                            <  8:234872/148‑1 (MQ=255)
                                               gCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACAGTGACGAACTAATTGCTCGTGTAATAGATAAAAATgg                                                                                            >  1:135837/1‑148 (MQ=255)
                                                cATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGt                                                                                           >  5:63364/1‑148 (MQ=255)
                                                 aTTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTaa                                                                                         <  8:212452/149‑1 (MQ=255)
                                                 aTTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTaa                                                                                         >  8:230506/1‑149 (MQ=255)
                                                     tttGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACaa                                                                                      <  3:110372/148‑1 (MQ=255)
                                                     tttGCTGTTTTAGCATTCAGTGTTTTTATCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACaa                                                                                      <  2:206661/148‑1 (MQ=255)
                                                      ttGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAata                                                                                    >  1:182270/1‑149 (MQ=255)
                                                         cTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGa                                                                                 <  7:121095/149‑1 (MQ=255)
                                                               tAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATAc                                                                            <  6:71514/148‑1 (MQ=255)
                                                                aGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTGGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTaa                                                                                                                     <  2:98888/106‑1 (MQ=255)
                                                                aGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTGGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTaa                                                                                                                     >  1:98888/1‑106 (MQ=255)
                                                                    ttCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCAt                                                                      >  7:99008/1‑149 (MQ=255)
                                                                     tCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAAtt                                                                                                                   >  5:65639/1‑103 (MQ=255)
                                                                     tCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAAtt                                                                                                                   <  6:65639/103‑1 (MQ=255)
                                                                       aGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCAttt                                                                    >  5:170729/1‑148 (MQ=255)
                                                                           tttttttCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATc                                                                <  8:99008/148‑1 (MQ=255)
                                                                               tttCTTAGTAGAGTATGTTTTAGGGCTGGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATGTGATCCTGTc                                                           >  4:55335/1‑149 (MQ=255)
                                                                                ttCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCt                                                          >  3:161509/1‑149 (MQ=255)
                                                                                   ttAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCttt                                                        >  2:24145/1‑148 (MQ=255)
                                                                                      gtagAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTATCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGc                                                    <  1:113037/149‑1 (MQ=255)
                                                                                       tagAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGcc                                                   >  2:170590/1‑149 (MQ=255)
                                                                                       tagAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGcc                                                   <  6:42816/149‑1 (MQ=255)
                                                                                         gagTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCt                                                  >  4:17795/1‑148 (MQ=255)
                                                                                            tATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGAt                                               >  3:63063/1‑148 (MQ=255)
                                                                                               gTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTaa                                                                                         >  8:216900/1‑103 (MQ=255)
                                                                                               gTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTaa                                                                                         <  7:216900/103‑1 (MQ=255)
                                                                                               gTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAAt                                            >  1:90171/1‑148 (MQ=255)
                                                                                                ttttAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGc                                          >  3:233217/1‑149 (MQ=255)
                                                                                                 tttAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCt                                         >  6:78225/1‑149 (MQ=255)
                                                                                                   tAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTg                                        <  1:24145/148‑1 (MQ=255)
                                                                                                    aGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTgg                                       <  4:63063/148‑1 (MQ=255)
                                                                                                     gggCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTgggg                                     <  3:17795/149‑1 (MQ=255)
                                                                                                     gggCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTggg                                      <  1:170590/148‑1 (MQ=255)
                                                                                                     gggCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATGTGATCCTGTCTTTAAGCCTGATAATGCTgggg                                     >  4:48727/1‑149 (MQ=255)
                                                                                                      ggCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGc                                    <  4:161509/149‑1 (MQ=255)
                                                                                                           gACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCAtt                                <  4:209824/148‑1 (MQ=255)
                                                                                                            aCGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTgg                              >  5:100883/1‑149 (MQ=255)
                                                                                                            aCGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGATAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTgg                              >  4:34604/1‑149 (MQ=255)
                                                                                                              ggAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGATGGGGCCATTGGc                             >  2:80178/1‑148 (MQ=255)
                                                                                                               gAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTc                           <  2:208215/149‑1 (MQ=255)
                                                                                                                  aTAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATg                        <  7:230506/149‑1 (MQ=255)
                                                                                                                   tAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTACTTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGc                       <  5:78225/149‑1 (MQ=255)
                                                                                                                    aaaaGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTggg                                      <  7:154923/133‑1 (MQ=255)
                                                                                                                    aaaaGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTggg                                      >  8:154923/1‑133 (MQ=255)
                                                                                                                     aaaGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCag                     >  2:45732/1‑149 (MQ=255)
                                                                                                                                     gTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGt     >  3:230214/1‑149 (MQ=255)
                                                                                                                                       tCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTc    >  5:183544/1‑148 (MQ=255)
                                                                                                                                         ttAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGa  >  3:105278/1‑148 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
CTCATCCTTTCAACAACGAGCACCTGACATCAGGTAATTGGATAATAGCATTTTTTGCTGTTTTAGCATTCAGTGTTTTTTTCTTAGTAGAGTATGTTTTAGGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGA  >  NC_000913/4635319‑4635603

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