Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 8.0% intergenic (+46/+56) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,6950CG 8.0% 241.6 / 6.3 92intergenic (+46/+56)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (40/40);  new base G (2/5);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAA  >  NC_000913/3608548‑3608842
                                                                                                                                                      |                                                                                                                                                    
acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCctct                                                                                                                                                     >  1:429278/1‑149 (MQ=255)
 cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTg                                                                                                                                                    <  1:214172/149‑1 (MQ=255)
       tAGGGATGACCGGGTTGTGGCTGGCAGTGCAGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                              <  2:804158/149‑1 (MQ=255)
        aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                             <  1:765177/149‑1 (MQ=255)
        aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGAAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                             <  1:623139/149‑1 (MQ=255)
         gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            <  3:638717/149‑1 (MQ=255)
         gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            <  2:742232/149‑1 (MQ=255)
          ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                           >  4:399208/1‑149 (MQ=255)
          ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                           >  4:1007721/1‑149 (MQ=255)
          ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                           >  1:374568/1‑149 (MQ=255)
           gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                          >  1:388584/1‑149 (MQ=255)
              gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑g                                                                                                                                      <  1:572187/149‑1 (MQ=255)
               aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCTCCTCTGGGGAGAGGGTTA‑gg                                                                                                                                     >  3:32236/1‑149 (MQ=255)
                ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑ggg                                                                                                                                    >  3:705863/1‑149 (MQ=255)
                 cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATAGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCGCTGGGGAGAGGGTTA‑GGGt                                                                                                                                   <  4:890360/149‑1 (MQ=255)
                     ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAgg                                                                                                                               <  4:736660/149‑1 (MQ=255)
                     ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAgg                                                                                                                               <  1:883717/149‑1 (MQ=255)
                     ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAgg                                                                                                                               <  3:780526/149‑1 (MQ=255)
                       gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAgggg                                                                                                                             >  4:556300/1‑149 (MQ=255)
                          gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            >  1:480730/1‑132 (MQ=255)
                          gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            <  2:480730/132‑1 (MQ=255)
                            tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTATGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAg                                                                                                                        >  2:679203/1‑149 (MQ=255)
                              gcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAggc                                                                                                                      >  4:709775/1‑149 (MQ=255)
                              gcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAggc                                                                                                                      <  2:429278/149‑1 (MQ=255)
                               cAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAggcg                                                                                                                     >  3:759769/1‑149 (MQ=255)
                               cAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGCCCCCTCGCCCCTCGGGGGAGAGGGTTA‑GGGTGAGGGGAAAAggcg                                                                                                                     <  1:57548/149‑1 (MQ=255)
                                       ttccgatcGGGGGCGACGGTGCTGGAGACAGGAAATGCTTTAAGATTGTGGCGCAGGAGAAAA‑GGC‑AAACCGATCGAAACATTGAGCGC‑GACCGGCCCCCCCGCCCCTGTGGGGAGGGGGTTG‑GGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                             <  2:916748/142‑1 (MQ=255)
                                       gCAGATGGGGGGGGGCGGGTGCTGGTGACAGCAAATGCTTAAAGATTGTTGCGAAGGAGAAAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGCCCCCGCGCCCCTCGGGGGAGGGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                             <  1:432473/149‑1 (MQ=255)
                                       gCAGATGGGGGGGCGAGGGTGGGGGTGAAAGCAAATGCGTTAGGATGGTTGCGAAGGAAAAAA‑GGC‑AACCCGATCGAAACATTGAGCGC‑GATCGGCCCCCCCCCCCCTCGGGGGAGGGGTTTA‑GGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                             <  2:487056/149‑1 (MQ=255)
                                       gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCTCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                             <  4:32236/149‑1 (MQ=255)
                                             cggggggCGAGGGTGCTGGAGACAGAAATGGCGTTAGGATTGTGGCGCGGGAAATAA‑GGC‑AACCCGATGGCAAATTGGGGCGA‑GACGGGCCCCCCCGCCCCTGGGGGGAGGGGGTTG‑GGTAGGGGGGAAAAGGCGGCATCGAAGCCAAt                                                                                                       <  2:96874/147‑1 (MQ=255)
                                                   gggcgggTGCTGGTACAAGCGAATGCTTTAAGTTTGTTGCGCAGGAGATAA‑GGC‑AAACGGATCGCAACATTGGGCGC‑GATCGGCCCCCTCGCCCCTCTGGGGAGGGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                 <  2:927460/144‑1 (MQ=37)
                                                   gCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                 >  4:482368/1‑149 (MQ=255)
                                                     gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                               >  1:315881/1‑149 (MQ=255)
                                                       cGGGGCTGGTGACAGGAAATGCGTTAAGATTGTTGGGGAGGAGATAA‑GGC‑AACCCGATCAAAAAATTGAGCGC‑GACCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTAt                                                                                             <  4:795563/149‑1 (MQ=37)
                                                           gggggggACGGCAAATGGTTTAGGTTTGTTGGGAGGGAGAAAG‑GGA‑AAACGGAGCGAAACATTGAGCGC‑GATGGGCCCCCCCGCCCCTTGGGGGGGAGGGTTA‑GGGGGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                         <  4:409650/143‑1 (MQ=255)
                                                                tGACAGCAAATGCGTTAAGTTTGTGGCGCAGAAGATAA‑GGC‑AAACCGATCGCAACATGGAGCGC‑GACCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                    <  2:656265/149‑1 (MQ=38)
                                                                ggaaaGCAAATGCTTTAAGATTTTGGCGCA‑GAGAAAAGGGC‑AACCCGATCGAAACATGGAGCGG‑GATGGGTCCCCTCGCCCCTTGGGGAAGAGGTTTA‑GGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                    <  2:435172/145‑1 (MQ=255)
                                                                 gACAGCGAATGCGTTAAATTTGTGGCGGAGGAGATAA‑GAC‑AAACCGTTCGAAACTTTGAGCGC‑GCTCGGTCCCCGCGCCCCTGTGGGGAG‑GGGTTAGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                   <  3:378/149‑1 (MQ=255)
                                                                  gcAGCAAATGCTTTAAGGTGTTTGCCAAGGAGAAAG‑GGC‑AACCCAACCGCAACATGGA‑CGCGGGGCGGGCCGCTCGCCCCTCGGGGAAGAGGGTTG‑GGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                  <  3:363632/148‑1 (MQ=255)
                                                                    ggagagtGCTTTAAGGTTGTCGCGCAGGGGATAA‑GGA‑AAACGGATCGAAGCATAGAGCGT‑GACCGCCCCCCCCGCCCCTTTGGGGGGAGGGTGG‑GGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                <  3:116238/143‑1 (MQ=255)
                                                                       gAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑gg                                                                                                                                     >  1:975560/1‑93 (MQ=255)
                                                                       gAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑gg                                                                                                                                     <  2:975560/93‑1 (MQ=255)
                                                                       gAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGATA‑GGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCATTACGaa                                                                             >  2:12391/1‑149 (MQ=255)
                                                                            cGTTAAGATTGTTGCGAAGGAGATAA‑GGA‑AAACGGATCGAAACATGGAGCGC‑GATCGGCCCCCTCCCCCCTCGGGGGAGGGGGTTA‑GGGGGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaaa                                                                        <  2:374568/149‑1 (MQ=255)
                                                                             gtaaggtttgttgCGCAGGAGAAAA‑GGC‑AAACCGATCGAAACATGGAGCGC‑GACCGGCCCCCCCGCCCCTCGGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                       <  4:100273/142‑1 (MQ=38)
                                                                             gTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                       >  3:100273/1‑149 (MQ=255)
                                                                              ttAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                      >  3:141048/1‑149 (MQ=255)
                                                                              ttAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                      >  3:405136/1‑149 (MQ=255)
                                                                               aaagttggttgCGAGGGAAATAA‑GAA‑AACCCCATGGAAAAATTGAGCGC‑GGCCGTCCCCCGCCCCCCTCGGGGGGGGGGGTTG‑GGTGGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                     <  3:482368/142‑1 (MQ=255)
                                                                               aaaGTTTGTTGCGAAGGAAATAA‑GGC‑AAACCGATCGAAACTTTGGGCGC‑GCTCGCCCCCCCCCCCCCTCGGGGGAGAGGGTTA‑GGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                     <  3:399208/148‑1 (MQ=255)
                                                                                aaGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAAc                                                                    >  3:374976/1‑149 (MQ=255)
                                                                                  gATTGTTGGGGAGAAGAAAA‑GGC‑AAGCCAATCGAAACTTTGGCGGC‑AACCGGCCCCTGCGCCCTTCGGGGGAGGGGTTGG‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                  <  3:400000/149‑1 (MQ=255)
                                                                                  gATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAGGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                  >  1:951434/1‑149 (MQ=255)
                                                                                  gATTGTTGCGCAGGAAAAAA‑GGA‑AAACCGATCGAAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                  <  2:973549/149‑1 (MQ=38)
                                                                                      gttgCGCAGGAAAAAA‑GCC‑AACCCGATCGAAACATGGAGCGC‑GATCGGCCCCCCCGCCCCTCTGGGGAGAGGGTTG‑GGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                              <  2:55712/149‑1 (MQ=255)
                                                                                       tggcgaAGGAGA‑AA‑GGCAAAACCGATCGAAACTTGGAGCGC‑GATCGACCCCCTCACCCGTCGGGGGGGGGGGTTA‑GGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                             <  1:512426/147‑1 (MQ=255)
                                                                                           gcAGGAGAAAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGCCCCCCTCCCCCCTCTGGGGAGGGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                         <  4:374976/149‑1 (MQ=39)
                                                                                             aGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTcc                                                       >  1:13207/1‑149 (MQ=255)
                                                                                                 gaTAA‑GGC‑AACCGGATCGAAAAAGTGGGCGC‑GATCGGCCCCCCCGCCCCTCGGGGGAGGGGGTAG‑GGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                   <  4:363686/149‑1 (MQ=255)
                                                                                                  aTAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                                                  >  3:426260/1‑149 (MQ=255)
                                                                                                        gC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                             >  1:458256/1‑149 (MQ=255)
                                                                                                           aaaCCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                           >  3:825907/1‑149 (MQ=255)
                                                                                                           aaaCCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                           >  4:442079/1‑149 (MQ=255)
                                                                                                            acccGATCGCAAAATTGGGCGG‑GCTCGCTCCCCCCGCCCCTGTGGGGAGGGGGTTG‑GGGGGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                          <  3:403094/147‑1 (MQ=255)
                                                                                                             aCCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                         >  2:989065/1‑149 (MQ=255)
                                                                                                              ccGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGc                                        >  3:957668/1‑149 (MQ=255)
                                                                                                               cGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTACATCCGTCTCTTTAGCa                                       >  4:883927/1‑149 (MQ=255)
                                                                                                                   cgaaaCATTGGGCGG‑GATCGGTCCCCTCCCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                   <  4:957668/146‑1 (MQ=255)
                                                                                                                   cgaaaCATGGGGCGC‑GATCGGCCCCCTCGCCCCTCGGGGGAGAGGGTTA‑GGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                   <  4:836046/146‑1 (MQ=255)
                                                                                                                     cAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                                 >  2:376929/1‑149 (MQ=255)
                                                                                                                         aTTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                             >  3:267262/1‑149 (MQ=255)
                                                                                                                         aTTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                             >  2:287471/1‑149 (MQ=255)
                                                                                                                             aGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                         >  2:442925/1‑149 (MQ=255)
                                                                                                                              gccc‑gATCGGCCCCCCCGCCCCTCGGGGAGGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                        <  1:679203/146‑1 (MQ=255)
                                                                                                                                   gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                    >  2:675189/1‑149 (MQ=255)
                                                                                                                                    aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                   >  1:925116/1‑149 (MQ=255)
                                                                                                                                    aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                   >  3:29420/1‑149 (MQ=255)
                                                                                                                                     tCGGTCCCCTCGCCTGCTTGGGGAGGGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                  <  3:355003/149‑1 (MQ=255)
                                                                                                                                     tCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                  >  2:313951/1‑149 (MQ=255)
                                                                                                                                     tCGGTCCCCCCCCCTCTCTGGGGAGGGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                  <  1:1005637/149‑1 (MQ=255)
                                                                                                                                       ggTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                <  4:405136/149‑1 (MQ=255)
                                                                                                                                       ctaacccTCTCCCCAGAGGGGCGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                  >  1:252319/5‑147 (MQ=255)
                                                                                                                                       ctaacccTCTCCCCAGAGGGGCGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                >  1:412217/5‑149 (MQ=255)
                                                                                                                                         tCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc              >  2:746202/1‑149 (MQ=255)
                                                                                                                                             cccGCCCCTTGGGGGGGAGGGTTA‑GGGAGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa          <  1:456443/147‑1 (MQ=255)
                                                                                                                                             cTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa          >  2:532382/1‑149 (MQ=255)
                                                                                                                                              tCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc         >  3:656539/1‑149 (MQ=255)
                                                                                                                                                gCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                    <  3:368138/72‑1 (MQ=255)
                                                                                                                                                gCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                    >  4:368138/1‑72 (MQ=255)
                                                                                                                                                 ccccTGTGGGAGGAGGGATA‑GGTTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa      <  4:1022359/149‑1 (MQ=255)
                                                                                                                                                 ccccTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa      <  4:759769/149‑1 (MQ=255)
                                                                                                                                                 ccccTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa      <  3:86261/149‑1 (MQ=255)
                                                                                                                                                     tctGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa  >  2:346342/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAA‑GGC‑AAACCGATCGCAACATTGAGCGC‑GATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTA‑GGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAA  >  NC_000913/3608548‑3608842

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