Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,631,492 G→T 8.6% pseudogene (794/1284 nt) ydfJ ← pseudogene, MFS transporter family; interrupted by Qin prophage;Phage or Prophage Related; putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,631,4920GT 8.6% 187.5 / 3.9 63pseudogene (794/1284 nt)ydfJpseudogene, MFS transporter family; interrupted by Qin prophage;Phage or Prophage Related; putative transport protein
Reads supporting (aligned to +/- strand):  ref base G (32/25);  new base T (3/3);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGCAACCATAATGGTGCTCGGGGCATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTACCCGCCTGACCAAAACGCAGCCCTGTTGCCAGCCA  >  NC_000913/1631345‑1631640
                                                                                                                                                   |                                                                                                                                                    
aGTGCAACCATAATGGTGCTCGGGGCATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATgg                                                                                                                                                      <  4:87243/148‑1 (MQ=255)
aGTGCAACCATAATGGTGCTCGGGGCATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATgg                                                                                                                                                      <  3:134001/148‑1 (MQ=255)
  tGCAACCATAATGGTGCTCGGGGCATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCa                                                                                                                                                   >  3:23034/1‑149 (MQ=255)
     aaCCATAATGGTGCTCGGGGCATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATa                                                                                                                                                 >  3:129373/1‑148 (MQ=255)
      aCCATAATGGTGCTCGGGGCATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGTTCATa                                                                                                                                                 >  1:43657/1‑147 (MQ=255)
      aCCATAATGGTGCTCGGGGCATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGTTCATa                                                                                                                                                 >  1:48844/1‑147 (MQ=255)
        cATAATGGTGCTCGGGGCATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGc                                                                                                                                             >  7:163799/1‑149 (MQ=255)
          tAATGGTGCTCGGGGCATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGcc                                                                                                                                            >  8:49997/1‑148 (MQ=255)
                 gctcttggCATAGCTTTTATCTAAAATTATTGAAAGCATTGGCCATGCCATAACAATCGCGGAGGTATTAATAATAATATACGTGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGATAAtg                                                                                                                                     <  1:271576/142‑2 (MQ=255)
                     ggggCATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGaa                                                                                                                                  <  1:87706/147‑1 (MQ=255)
                         cATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGatc                                                                                                                            <  4:23034/149‑1 (MQ=255)
                             gCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGAtcatca                                                                                                                        <  7:57312/149‑1 (MQ=255)
                              cTTTTATCTACAATGATAGACAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCaa                                                                                                                       <  2:128466/149‑1 (MQ=255)
                              cTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCaa                                                                                                                       <  5:193838/149‑1 (MQ=255)
                              cTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTAGTCAAAATAATATACGGGATCCGGCGCCCAATTTTCTCGGTTAACCATCCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGAGCATCaa                                                                                                                       <  4:313535/149‑1 (MQ=255)
                               tcttATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAAt                                                                                                                      >  1:332246/3‑149 (MQ=255)
                                tttATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGTTCATAAAGCCGAGAATCGAACTGATCATCaaag                                                                                                                     >  2:208436/1‑147 (MQ=255)
                                       aCAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCAAAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGt                                                                                                              >  5:223012/1‑149 (MQ=255)
                                           tGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTgg                                                                                                           >  4:228423/1‑148 (MQ=255)
                                              tAGAAAGCATTGGCCATGCCAGCACAATCTCGGAGGTATTCATAATAATATACGGTATCCGGCGACCAATTTTATCTGCTAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAAtt                                                                                                       <  5:80523/149‑1 (MQ=255)
                                                gAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGc                                                                                                     >  1:94887/1‑149 (MQ=255)
                                                 aaaGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCt                                                                                                    >  4:282816/1‑149 (MQ=255)
                                                   aGCATTGGCCATGCCAGCACAATCGCGTAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCttt                                                                                                  <  2:277174/149‑1 (MQ=255)
                                                   aGCATTGGCCATGCCAGCACAATCGCGGAGGTATTAATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCttt                                                                                                  <  5:23106/149‑1 (MQ=255)
                                                    gCATTTGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTg                                                                                                 <  1:116203/149‑1 (MQ=255)
                                                    gCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTg                                                                                                 <  4:328755/149‑1 (MQ=255)
                                                      aTTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGTATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCTTAAATCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGtt                                                                                               <  3:60578/149‑1 (MQ=255)
                                                          gCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTaaa                                                                                            <  7:214196/148‑1 (MQ=255)
                                                          gCCATGCCAGCACAATAGAGGAGGTATTCATTATAATATACGGGATCCGGCGACCAATTTTTTCGGATAACCAAACAAGGAACGGAATGGTCATAAAGCCTAGAATCGATCTGATCATAATTGCATCTGTTGGAATTGGTTTGTTaaa                                                                                            >  4:182301/1‑148 (MQ=255)
                                                              tGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATa                                                                                       <  2:332246/149‑1 (MQ=255)
                                                              tGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATAGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAAt                                                                                        >  2:156031/1‑148 (MQ=255)
                                                                     acaATCGCGTAGGTATTCATAATAATATACGGGATCCGGAGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGTTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTg                                                                                <  3:185018/149‑1 (MQ=255)
                                                                         tCGCGGAGGTATTCATAATAATATACGGGATCTGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCttt                                                                                                  <  7:184168/127‑1 (MQ=255)
                                                                         tCGCGGAGGTATTCATAATAATATACGGGATCTGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCttt                                                                                                  >  8:184168/1‑127 (MQ=255)
                                                                         tCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGATTTGCTTTGTTAAACAATAACGTCTGCAc                                                                             >  5:186295/1‑148 (MQ=255)
                                                                         tCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCAc                                                                             >  4:109939/1‑148 (MQ=255)
                                                                            cgGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCgag                                                                                                                                         >  2:88061/1‑85 (MQ=255)
                                                                            cgGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTaaa                                                                         <  4:129373/149‑1 (MQ=255)
                                                                            cgGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTaaa                                                                         <  5:174289/149‑1 (MQ=255)
                                                                                    ttCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTg                                                                  >  3:228167/1‑148 (MQ=255)
                                                                                       ataataataTACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCaa                                                               <  6:223012/148‑1 (MQ=255)
                                                                                         aataataTACGGGATCCGGCGACCAATTTTATCGGATTACCAACCAAGGAACGGAATGGTCATAACGCCGAGAATCGAACTGCTCATCAATGCATCTGTTGGAATTGCTTAGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAgg                                                             >  3:18436/1‑148 (MQ=255)
                                                                                           taataTACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGa                                                            >  1:159252/1‑147 (MQ=255)
                                                                                                 aCGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTg                                                     >  3:208429/1‑148 (MQ=255)
                                                                                                    ggATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAAt                                                  >  7:302044/1‑148 (MQ=255)
                                                                                                      aTCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTaa                                               <  8:302044/149‑1 (MQ=255)
                                                                                                        ccGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAAcc                                             >  5:260580/1‑149 (MQ=255)
                                                                                                         cGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAAccc                                            <  7:49997/149‑1 (MQ=255)
                                                                                                          ggCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGTTCATAAAGCCGAGAATCTAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCg                                           <  2:48844/149‑1 (MQ=255)
                                                                                                          ggCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGTTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCg                                           <  2:43657/149‑1 (MQ=255)
                                                                                                            cGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGa                                         <  2:94887/149‑1 (MQ=255)
                                                                                                             gACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACGGATCATCAATGCATCTGTTGGAATTGCTGTGTTAAACAATAACGTCTGCACTAAATAGCCTGAAAGGAAAGTCTGAATTAACCCGGAg                                        >  7:317287/1‑149 (MQ=255)
                                                                                                                  aTTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTAcc                                   >  5:269635/1‑149 (MQ=255)
                                                                                                                   ttttATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTAccc                                  >  8:210353/1‑149 (MQ=255)
                                                                                                                        tCGGATAACCAACCAAGGAACGGAATGGGCATAAAGCCTAGAATCGAACTGATCATCAATGAATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGTAGTTACCCGCCt                              >  6:123705/1‑148 (MQ=255)
                                                                                                                          ggATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAGTCTGAATTAACCCGGAGTTACCCGCCTGAcc                          <  1:156031/149‑1 (MQ=255)
                                                                                                                           gATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTACCCGCCTGAc                           >  1:99660/1‑148 (MQ=255)
                                                                                                                           gATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAAAGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGATACCCGCCTGAc                           >  2:38517/1‑148 (MQ=255)
                                                                                                                              aaccaaccAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTACCCGCCTGACCaa                        >  2:117740/1‑148 (MQ=255)
                                                                                                                                 caaccaaGGAACGGAATGGTCATAAAGCCGTGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTACCCGCCTGACCAAAACg                    >  1:181446/1‑149 (MQ=255)
                                                                                                                                  aaccaaGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTACCCGCCTGACCAAAACg                    >  5:68591/1‑148 (MQ=255)
                                                                                                                                  aaccaaGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTACCCGCCTGACCAAAACg                    >  5:308644/1‑148 (MQ=255)
                                                                                                                                  aaccaaGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTACCCGCCTGACCAAAACGc                   >  6:238704/1‑149 (MQ=255)
                                                                                                                                   accaaGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTACCCGCCTGACCAAAACGCa                  <  7:151124/149‑1 (MQ=255)
                                                                                                                                            cGGAATGGTCATAAAGCCGTGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTACCCGCCTGACCAAAACGCAGCCCTGTTg         <  2:181446/149‑1 (MQ=255)
                                                                                                                                                   gTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTACCCGCCTGACCAAAACGCAGCCCTGTTgccagcca  <  3:228423/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGTGCAACCATAATGGTGCTCGGGGCATAGCTTTTATCTACAATGATAGAAAGCATTGGCCATGCCAGCACAATCGCGGAGGTATTCATAATAATATACGGGATCCGGCGACCAATTTTATCGGATAACCAACCAAGGAACGGAATGGTCATAAAGCCGAGAATCGAACTGATCATCAATGCATCTGTTGGAATTGCTTTGTTAAACAATAACGTCTGCACTAAATAGCCTGCAAGGAAAGTCTGAATTAACCCGGAGTTACCCGCCTGACCAAAACGCAGCCCTGTTGCCAGCCA  >  NC_000913/1631345‑1631640

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