Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,375,067 A→G 10.8% A35A (GCA→GCG ycjQ → putative Zn‑dependent NAD(P)‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,375,0670AG10.8% 256.1 / 18.3 83A35A (GCA→GCGycjQputative Zn‑dependent NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base G (5/4);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GCTAAAGCGGCCAGTTTCTTTACAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATG  >  NC_000913/1374928‑1375214
                                                                                                                                           |                                                                                                                                                   
gCTAAAGCGGCCAGTTTCTTTACAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAaca                                                                                                                                             <  4:476783/148‑1 (MQ=255)
   aaaGCGGCCAGTTTCTTTACAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACacgg                                                                                                                                          >  1:509490/1‑148 (MQ=255)
               ttCTTTACAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGtcga                                                                                                                              >  1:92864/1‑145 (MQ=255)
                        aGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTCGACCACTTCcg                                                                                                                     <  2:92864/148‑1 (MQ=255)
                         ggAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACacgga                                                                                                                                         <  2:478119/127‑1 (MQ=255)
                         ggAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACacgga                                                                                                                                         >  1:478119/1‑127 (MQ=255)
                         ggAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCcgc                                                                                                                    >  2:986288/1‑148 (MQ=255)
                              tttAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGTGCGCCGaaa                                                                                                                                               <  1:305301/116‑1 (MQ=255)
                              tttAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGTGCGCCGaaa                                                                                                                                               >  2:305301/1‑116 (MQ=255)
                                   aTGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGcc                                                                                                          >  2:131709/1‑148 (MQ=255)
                                        aaaGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGttt                                                                                                     >  4:757101/1‑148 (MQ=255)
                                         aaGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGTGCGCCGaaa                                                                                                                                               <  2:651814/105‑1 (MQ=255)
                                         aaGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGTGCGCCGaaa                                                                                                                                               >  1:651814/1‑105 (MQ=255)
                                             tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGcc                                                                                                          <  1:456549/138‑1 (MQ=255)
                                             tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGcc                                                                                                          >  2:456549/1‑138 (MQ=255)
                                             tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtga                                                                                                <  1:986288/148‑1 (MQ=255)
                                             tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtga                                                                                                <  3:757101/148‑1 (MQ=255)
                                             tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtga                                                                                                >  4:368867/1‑148 (MQ=255)
                                             tagtagCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtga                                                                                                <  1:131709/148‑1 (MQ=255)
                                                  gCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGCGGTCGACTTCCGCGCCGCCAGCCCGTTTATtg                                                                                                 >  4:222971/1‑142 (MQ=255)
                                                  gCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGCGGTCGACTTCCGCGCCGCCAGCCCGTTTATtg                                                                                                 <  3:222971/142‑1 (MQ=255)
                                                   ccACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGcc                                                                                                          >  4:956502/1‑132 (MQ=255)
                                                   ccACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGcc                                                                                                          <  3:956502/132‑1 (MQ=255)
                                                     acaGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTCGACCACTTcc                                                                                                                      >  1:849876/1‑118 (MQ=255)
                                                     acaGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTCGACCACTTcc                                                                                                                      <  2:849876/118‑1 (MQ=255)
                                                      caGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGaa                                                                                                                                                >  4:515696/1‑91 (MQ=255)
                                                      caGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGaa                                                                                                                                                <  3:515696/91‑1 (MQ=255)
                                                      caGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACtt                                                                                       >  3:61554/1‑148 (MQ=255)
                                                          aCCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAAc                                                                                   <  2:509490/148‑1 (MQ=255)
                                                            cgcgTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCg                                                                                                                            <  1:514727/105‑1 (MQ=255)
                                                            cgcgTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCg                                                                                                                            >  2:514727/1‑105 (MQ=255)
                                                                   tgctgcGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATgg                                                                          <  4:61554/148‑1 (MQ=255)
                                                                    gctgcGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGc                                                                         >  1:684755/1‑148 (MQ=255)
                                                                       gcgcTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGa                                                                      <  3:368867/148‑1 (MQ=255)
                                                                            ggTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTc                                                                 <  2:684755/148‑1 (MQ=255)
                                                                              ttGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtg                                                                                                 >  1:356778/1‑114 (MQ=255)
                                                                              ttGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtg                                                                                                 <  2:356778/114‑1 (MQ=255)
                                                                              ttGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGa                                                                              >  3:524133/1‑133 (MQ=255)
                                                                              ttGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGa                                                                              <  4:524133/133‑1 (MQ=255)
                                                                                      gAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGt                                                                                                                              >  4:1003371/1‑77 (MQ=255)
                                                                                      gAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGt                                                                                                                              <  3:1003371/77‑1 (MQ=255)
                                                                                      gAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGa                                                                              <  4:822938/125‑1 (MQ=255)
                                                                                      gAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGa                                                                              >  3:822938/1‑125 (MQ=255)
                                                                                          aTCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTa                                                                                     <  2:975317/114‑1 (MQ=255)
                                                                                          aTCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTa                                                                                     >  1:975317/1‑114 (MQ=255)
                                                                                           tCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACacggaacgga                                                                                                                                    >  4:795753/1‑66 (MQ=255)
                                                                                           tCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACacggaacgga                                                                                                                                    <  3:795753/66‑1 (MQ=255)
                                                                                           tCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCa                                                  >  3:377067/1‑148 (MQ=255)
                                                                                             gggCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCcg                                                           <  3:1026069/137‑1 (MQ=255)
                                                                                             gggCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCcg                                                           >  4:1026069/1‑137 (MQ=255)
                                                                                              ggCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACacgga                                                                                                                                         <  3:562975/58‑1 (MQ=255)
                                                                                              ggCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACacgga                                                                                                                                         >  4:562975/1‑58 (MQ=255)
                                                                                                cGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCa                                                                <  2:76934/129‑1 (MQ=255)
                                                                                                cGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCa                                                                >  1:76934/1‑129 (MQ=255)
                                                                                                  aTTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATgg                                                                          >  2:37049/1‑117 (MQ=255)
                                                                                                  aTTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATgg                                                                          <  1:37049/117‑1 (MQ=255)
                                                                                                  aTTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATgg                                                                          >  4:1137018/1‑117 (MQ=255)
                                                                                                  aTTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATgg                                                                          <  3:1137018/117‑1 (MQ=255)
                                                                                                           aaTGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATg                                                                           >  3:502654/1‑107 (MQ=255)
                                                                                                           aaTGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATg                                                                           <  4:502654/107‑1 (MQ=255)
                                                                                                            aTGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCAt                                 <  4:377067/148‑1 (MQ=255)
                                                                                                             tgaagtgaagATCCGCGTGCGTTTCGGCGCACCGAAACacgga                                                                                                                                         >  2:338910/1‑43 (MQ=255)
                                                                                                             tgaagtgaagATCCGCGTGCGTTTCGGCGCACCGAAACacgga                                                                                                                                         <  1:338910/43‑1 (MQ=255)
                                                                                                             tgaagtgaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCAtt                                >  2:110050/1‑148 (MQ=255)
                                                                                                              gaagtgaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTg                               >  4:128576/1‑148 (MQ=255)
                                                                                                                  tgaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCcg                                                           <  3:739188/116‑1 (MQ=255)
                                                                                                                  tgaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCcg                                                           >  4:739188/1‑116 (MQ=255)
                                                                                                                  tgaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGtt                           <  1:110050/148‑1 (MQ=255)
                                                                                                                   gaagATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGttt                          <  3:128576/148‑1 (MQ=255)
                                                                                                                     agATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACt                                                                                        >  3:623651/1‑84 (MQ=255)
                                                                                                                     agATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACt                                                                                        <  4:623651/84‑1 (MQ=255)
                                                                                                                     agATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTgg                        >  1:971261/1‑148 (MQ=255)
                                                                                                                       aTCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtgatga                                                                                             <  2:1108934/77‑1 (MQ=255)
                                                                                                                       aTCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATtgatga                                                                                             >  1:1108934/1‑77 (MQ=255)
                                                                                                                                tgtTTCGGTGCGCCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGa                                                                      <  1:880386/90‑1 (MQ=255)
                                                                                                                                tgtTTCGGTGCGCCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGa                                                                      >  2:880386/2‑91 (MQ=255)
                                                                                                                                   ttCGGTGCGCCGAAACACGGAACGGAAGTGGTCGACTTCCGCgccgcc                                                                                                              >  4:425938/1‑48 (MQ=38)
                                                                                                                                   ttCGGTGCGCCGAAACACGGAACGGAAGTGGTCGACTTCCGCgccgcc                                                                                                              <  3:425938/48‑1 (MQ=38)
                                                                                                                                   ttCGGTGCGCCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTg          >  4:690592/1‑148 (MQ=255)
                                                                                                                                    tCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTgg         <  2:971261/148‑1 (MQ=255)
                                                                                                                                          cACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGttt                          >  2:956633/1‑125 (MQ=255)
                                                                                                                                          cACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACAt   >  3:367869/1‑148 (MQ=255)
                                                                                                                                          cACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAAttt                          <  1:956633/125‑1 (MQ=255)
                                                                                                                                           aCCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACAtg  <  4:367869/148‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GCTAAAGCGGCCAGTTTCTTTACAAGGAGTTTTAAATGAAAAAGTTAGTAGCCACAGCACCGCGTGTTGCTGCGCTGGTTGAGTATGAAGATCGGGCGATTTTAGCTAATGAAGTGAAGATCCGCGTGCGTTTCGGCGCACCGAAACACGGAACGGAAGTGGTCGACTTCCGCGCCGCCAGCCCGTTTATTGATGAAGACTTTAACGGCGAATGGCAGATGTTCACTCCGCGTCCGGCAGATGCGCCGCGCGGCATTGAGTTTGGCAAATTCCAGCTTGGCAACATG  >  NC_000913/1374928‑1375214

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