Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,348,977 C→T 6.2% intergenic (‑233/+338) dcuB ← / ← dcuR C4‑dicarboxylate transporter, anaerobic; DcuS co‑sensor/response regulator in two‑component regulatory system with DcuS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,348,9770CT 6.2% 206.6 / 3.3 64intergenic (‑233/+338)dcuB/dcuRC4‑dicarboxylate transporter, anaerobic; DcuS co‑sensor/response regulator in two‑component regulatory system with DcuS
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base T (2/2);  total (32/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 = 2.72e-01

ATATCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGAC  >  NC_000913/4348838‑4349119
                                                                                                                                           |                                                                                                                                              
atatCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAACATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGAt                                                                                                                                       >  4:348385/1‑149 (MQ=255)
 tatCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATg                                                                                                                                      >  2:409118/1‑149 (MQ=255)
  atCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGt                                                                                                                                     >  4:139716/1‑149 (MQ=255)
  atCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGt                                                                                                                                     <  4:458020/149‑1 (MQ=255)
      gTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAAt                                                                                                                                 <  3:62842/149‑1 (MQ=255)
           gTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAAtgtttt                                                                                                                            <  1:355378/149‑1 (MQ=255)
            tCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTATTGCGCAGATGTTAAtgttttg                                                                                                                           >  1:534558/1‑149 (MQ=255)
            tCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTATTGCGCAGATGTTAAtgttttg                                                                                                                           >  1:569066/1‑149 (MQ=255)
            tCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAAtgttttg                                                                                                                           >  3:327714/1‑149 (MQ=255)
              gTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAAtgttttgtt                                                                                                                         <  3:284366/149‑1 (MQ=255)
              gTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAAtgttttgtt                                                                                                                         >  4:145027/1‑149 (MQ=255)
              gTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAAtgttttgtt                                                                                                                         <  2:170245/149‑1 (MQ=255)
              gTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAAtgttttgtt                                                                                                                         >  2:106455/1‑149 (MQ=255)
                  ccGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCt                                                                                                                     <  4:535381/149‑1 (MQ=255)
                   cGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTg                                                                                                                    >  2:104832/1‑149 (MQ=255)
                      cTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTa                                                                                                                 <  3:296365/149‑1 (MQ=255)
                           gCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTa                                                                                                            <  2:238255/149‑1 (MQ=255)
                           gCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTa                                                                                                            >  4:85215/1‑149 (MQ=255)
                             aaaTACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTAtt                                                                                                          <  1:366267/149‑1 (MQ=255)
                                   cTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGTGCAGATGTGAATGTTTTGTTTGCTGTTATCTTATTGATAtt                                                                                                    >  4:461838/1‑149 (MQ=255)
                                   cTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCa                                                                                                                                          >  3:88298/1‑111 (MQ=255)
                                   cTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCa                                                                                                                                          <  4:88298/111‑1 (MQ=255)
                                     cacaGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTa                                                                                                  >  2:254282/1‑149 (MQ=255)
                                      acaGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGata                                                                                                      <  4:20500/144‑1 (MQ=255)
                                      acaGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTCGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGata                                                                                                      >  3:20500/1‑144 (MQ=255)
                                       caGTGAATATTGGCTGATTAAAGCAGCGGCAAACTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATc                                                                                                <  4:375306/149‑1 (MQ=255)
                                        aGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCg                                                                                               <  4:496042/149‑1 (MQ=255)
                                         gTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCgg                                                                                              <  4:242756/149‑1 (MQ=255)
                                         gTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCgg                                                                                              >  1:218543/1‑149 (MQ=255)
                                         gTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACGGCGCAGATGCTAATGTTTTGTTTGCTGTTATATTATTGATATTTATCgg                                                                                              >  4:478035/1‑149 (MQ=255)
                                                ttGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTAttt                                                                                       <  2:218543/149‑1 (MQ=255)
                                                   gCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTAtttttt                                                                                    <  4:433807/149‑1 (MQ=255)
                                                   gCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTAtttttt                                                                                    <  4:502694/149‑1 (MQ=255)
                                                   gCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTAtttttt                                                                                    >  4:120058/1‑149 (MQ=255)
                                                      gATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGtt                                                                                 >  4:504739/1‑149 (MQ=255)
                                                          aaaGCAGCGGCAAATTTATGCACTGTGATAAGCGTCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTGATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCa                                                                             <  1:479142/149‑1 (MQ=255)
                                                            agcagcagcAAATTTATTCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAgg                                                                           <  1:299268/149‑1 (MQ=255)
                                                             gcagcGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAggg                                                                          >  3:4732/1‑149 (MQ=255)
                                                                      aaTTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCa                                                                 <  2:3918/149‑1 (MQ=255)
                                                                                cTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTtat                                                       <  4:415224/149‑1 (MQ=255)
                                                                                 tgtgATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTtata                                                      <  2:359141/149‑1 (MQ=255)
                                                                                     aTAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGAtt                                                  >  4:41151/1‑149 (MQ=255)
                                                                                         gCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTGATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATc                                              >  4:336157/1‑149 (MQ=255)
                                                                                           ggCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCGTTAGGTTTATAGATTAATCTg                                            >  3:119543/1‑149 (MQ=255)
                                                                                             cTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGAt                                          <  1:321149/149‑1 (MQ=255)
                                                                                                  tCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTAcc                                     >  3:508001/1‑149 (MQ=255)
                                                                                                       gcgaAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCAttt                                <  3:461838/147‑1 (MQ=255)
                                                                                                          aaaGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTtgtg                             >  4:377230/1‑149 (MQ=255)
                                                                                                             gACTCTCTGGCGATGGGTTTATAGATTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGt                          >  2:455537/1‑149 (MQ=255)
                                                                                                                tctctGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTaaa                       >  3:141621/1‑149 (MQ=255)
                                                                                                                  tctGGCGATGGGTTTATAGTTTGTATTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTaaaaa                     <  2:534558/149‑1 (MQ=255)
                                                                                                                  tctGGCGATGGGTTTATAGTTTGTATTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTaaaaa                     <  2:569066/149‑1 (MQ=255)
                                                                                                                  tctGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAAGTTTATAGATTAATCTGATCTACCCATTTGTGGGTaaaaa                     <  4:327714/149‑1 (MQ=255)
                                                                                                                     ggCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATac                  <  3:560754/149‑1 (MQ=255)
                                                                                                                       cGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATacac                <  1:326716/149‑1 (MQ=255)
                                                                                                                       cGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATacac                <  1:254282/149‑1 (MQ=255)
                                                                                                                        gATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATAcaca               >  3:560278/1‑149 (MQ=255)
                                                                                                                         aTGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACAt              >  2:128556/1‑149 (MQ=255)
                                                                                                                              tttATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGc         >  3:332636/1‑149 (MQ=255)
                                                                                                                              tatataGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATAGGTCTATTTTTTGTTTTCAGGTAAGTTGTCTTTAGGTTTATAGATTATTTTGATTTACCCAGTTGTGGGTAAAAATACTCATagtta         >  4:159154/3‑145 (MQ=255)
                                                                                                                                tataGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCgg       <  4:332636/149‑1 (MQ=255)
                                                                                                                                 ataGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCggg      <  1:228671/149‑1 (MQ=255)
                                                                                                                                     tttGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGAc  >  4:187491/1‑149 (MQ=255)
                                                                                                                                     tttGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGAc  >  1:18855/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                              
ATATCCGTTTGGTCGTTACCGGCTTTAGCAAATACCTCACAGTGAATATTGGCTGATTAAAGCAGCGGCAAATTTATGCACTGTGATAAGCGGCTTTTTCAGAGAGAAAGACTCTCTGGCGATGGGTTTATAGTTTGTACTGCGCAGATGTTAATGTTTTGTTTGCTGTTATCTTATTGATATTTATCGGTCTATTTTTTGTTCTCAGGGCAGATGTCATTAGGTTTATAGATTAATCTGATCTACCCATTTGTGGGTAAAAATACACATAATGCGGGTGAC  >  NC_000913/4348838‑4349119

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