Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 9.2% 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 9.2% 248.6 / 11.6 86intergenic (+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 (37/40);  new base G (4/4);  total (41/45)
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

ACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAG  >  NC_000913/3608548‑3608839
                                                                                                                                                    |                                                                                                                                                 
acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                  >  3:139364/1‑149 (MQ=255)
acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctc                                                                                                                                                   >  2:432853/1‑148 (MQ=255)
   cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                               <  4:816949/149‑1 (MQ=255)
     gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                              >  1:754217/1‑148 (MQ=255)
     gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                             <  3:265978/149‑1 (MQ=255)
      ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                <  4:130433/145‑1 (MQ=255)
      ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                >  3:130433/1‑145 (MQ=255)
        aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGag‑ag                                                                                                                                         <  3:580578/149‑1 (MQ=255)
         gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCTCCTCTGGGGAG‑Agg                                                                                                                                        <  3:67490/149‑1 (MQ=255)
         gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGag‑ag                                                                                                                                         <  1:319121/148‑1 (MQ=255)
         gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGag‑ag                                                                                                                                         <  2:226194/148‑1 (MQ=255)
           gATGACCGGGTTGTGGTTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑Aggg                                                                                                                                       >  1:580439/1‑148 (MQ=255)
           gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGt                                                                                                                                      >  3:30833/1‑149 (MQ=255)
             tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGtt                                                                                                                                     <  2:510813/148‑1 (MQ=255)
              gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCG‑Aggg                                                                                                                                       <  4:903360/145‑1 (MQ=255)
               aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAgg                                                                                                                                  <  4:470791/149‑1 (MQ=255)
               aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAgg                                                                                                                                  >  4:112758/1‑149 (MQ=255)
                 cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAggg                                                                                                                                 <  2:148691/148‑1 (MQ=255)
                 cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGt                                                                                                                                <  3:644412/149‑1 (MQ=255)
                  gggTTGTGGCTGGCAGTGCTGGCAGATGCGGGGGCGACGGTGCTGGTGAAAACGAATGCGATAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GAGCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGACAggga                                                                                                                                >  2:689923/1‑147 (MQ=255)
                  gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGt                                                                                                                                >  2:197909/1‑148 (MQ=255)
                    gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAg                                                                                                                             <  3:781725/149‑1 (MQ=255)
                    gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAg                                                                                                                             <  3:238879/149‑1 (MQ=255)
                    gTTGTGGCTGGCAGTGCTGGCACATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGCTAGGGTGa                                                                                                                              <  1:783568/148‑1 (MQ=255)
                     ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAgg                                                                                                                            <  4:785094/149‑1 (MQ=255)
                     ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAgg                                                                                                                            <  3:309811/149‑1 (MQ=255)
                           ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGaaa                                                                                                                       <  1:805179/148‑1 (MQ=255)
                             ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAgg                                                                                                                                  <  3:296398/135‑1 (MQ=255)
                             ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAgg                                                                                                                                  >  4:296398/1‑135 (MQ=255)
                             ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAgg                                                                                                                    >  3:224344/1‑149 (MQ=255)
                                    cTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATc                                                                                                             >  3:365480/1‑149 (MQ=255)
                                        cAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGaa                                                                                                          >  2:737947/1‑148 (MQ=255)
                                            taggggggCGAGGGTGCGGGTAACAGCGAAGGCGTAAAGATTGTTGCCACGGAGAAAAGGCAACCG‑GGCAGCAACATTCGGCGCGACCGGCCCCCCCGCCCCTCTGGGGGG‑GGGGTTCGGGTGAGGGAAAAAGGCGGCATCGAAGCCaa                                                                                                     <  4:365480/146‑1 (MQ=255)
                                                   gggagggtggtGGGGACAGCGAAGGGGTTAAGGTTGTTGCGAGGGAAATAAGGAAAACC‑AACCCAAACATGGAGGGCATACGTTCCCCTCGCCCCGCTGGGGAG‑AGGGTGGGGTTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                              <  4:895378/144‑1 (MQ=255)
                                                    cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                              >  1:366333/1‑148 (MQ=255)
                                                     gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                            >  4:667718/1‑149 (MQ=255)
                                                       cGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                           >  1:648394/1‑148 (MQ=255)
                                                              ggTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTaagatcgg                                                                                    >  1:664310/1‑141 (MQ=255)
                                                               gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                   >  1:564041/1‑148 (MQ=255)
                                                                   cggaaaaTGGGTTAAGTTCTTTGCGGGGGGGA‑AAGGAAAACCAGATCACAAAATTGAGCGTGATCGGTCCCCTCGCCCCGC‑GGGGAGAAGGGTTGGGCGGGGAGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                              <  3:112758/144‑1 (MQ=255)
                                                                      caaaTGCTTTAAGATTGTGGCGAAGGAGAAAAGGAAAACC‑GATCGCAACATTGGGCGCGACCGCCCCCCCCCCCCTTGGGGGGAG‑GGGGTTGGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                            <  1:737947/146‑1 (MQ=255)
                                                                       gAATGCTTTAAGTTTGTGGCGCAGGAGAAAAGGAAACCC‑GACCGAAAAATGGAGCGCGATCGGCCCCCCCGCCCCTCGGGGGAG‑AGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                           <  1:432853/148‑1 (MQ=255)
                                                                        aaTGCGTTAAGTTTGTTGCGAAGGAAATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGCCCCCCCGCCCCTCTGGGGAG‑GGGGTTAGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                         <  4:823348/149‑1 (MQ=255)
                                                                              ttAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAgg                                                                                                                    >  1:651202/1‑100 (MQ=255)
                                                                              ttAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAgg                                                                                                                    <  2:651202/100‑1 (MQ=255)
                                                                                   aTTTTTGCGAAGAAAATAAGGCAAACC‑GATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑GGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAt                                                              <  3:884661/149‑1 (MQ=38)
                                                                                       ttgCGCAGGAGATAAGGAAAACG‑GAGCGAAACATTGCGCGCGATCGGCCCCCCCGCCCCTCGGGGGAG‑GGGGTTGGGGGGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                          <  4:139364/149‑1 (MQ=255)
                                                                                       tcgggaaGGAGAAAAGAAAAACC‑AAACAAAATATGGCGCGCGAGCGCCCCCCCCGGCCCTCGGGGGGG‑GGGGTGAGGGTGAGGAGAAAAGGCGGCCTCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAg                                                           <  1:831349/142‑1 (MQ=255)
                                                                                           gcAGGAAAGAAGGAAAACC‑GATCGAAACATGGAGCGCGATCGGCCCCCTCGCCCCTCGGGGGAG‑AGGGTTGGGGTGGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                      <  4:30833/149‑1 (MQ=255)
                                                                                             aGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTc                                                     >  2:36134/1‑148 (MQ=255)
                                                                                                   tAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                              >  3:893964/1‑149 (MQ=255)
                                                                                                   tAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGCCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCATTTACGCATCTAATAACGATAAGGCAGTCCATCCGt                                              >  4:560097/1‑149 (MQ=39)
                                                                                                    aaGGCAAACC‑GATCGCAACATTGGGCGCGATCGGCCCCCTCGCCCTTCGGGGGAG‑AGGTTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                              <  1:60639/148‑1 (MQ=255)
                                                                                                        cAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                          <  2:754217/148‑1 (MQ=255)
                                                                                                          aaCC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                        >  1:458674/1‑148 (MQ=255)
                                                                                                          aaCC‑GATCGAAACATTGAGCGGGATCGGCCCCCCCGCCCCTCGGGGGGG‑GGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                        <  2:366333/148‑1 (MQ=255)
                                                                                                           ccc‑GATCAAAATATTGGGGGCGAGCGGCCCCCCGGCCCCTGTGGGGAG‑GGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                       <  2:315070/147‑1 (MQ=255)
                                                                                                                aTCGCACAATGGAGCGCGATCGGCCCCCCCGCCCCTCGGGGGAG‑AGGGTTAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                  <  3:560097/149‑1 (MQ=255)
                                                                                                                    cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                              >  3:192878/1‑149 (MQ=255)
                                                                                                                      actttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                             <  2:564041/145‑1 (MQ=255)
                                                                                                                        caatgttgcgATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCg                           >  1:845521/8‑148 (MQ=255)
                                                                                                                            aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGaaa                                                                                                                       >  2:910059/1‑52 (MQ=255)
                                                                                                                            aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGaaa                                                                                                                       <  1:910059/52‑1 (MQ=255)
                                                                                                                             gcgcGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                     >  3:538243/1‑149 (MQ=255)
                                                                                                                                 gATCGGTCCCCTCGCCCCTCTGGGGGG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                 >  4:479188/1‑149 (MQ=255)
                                                                                                                                 gATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                 >  3:659032/1‑149 (MQ=255)
                                                                                                                                  aTCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                 >  2:280899/1‑148 (MQ=255)
                                                                                                                                  aTCGGTCCCCTCGCCCCTCAGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAGTAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCACTTCGTGTTCCATaa                 >  2:81827/1‑148 (MQ=255)
                                                                                                                                   gcGGTCCCCTCGCCCCCATGGGGAA‑GGGGTTAGGGATAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                <  1:689923/147‑1 (MQ=255)
                                                                                                                                     ctaacccTCTCCCCAGAGGGGCG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg              >  1:381634/5‑148 (MQ=255)
                                                                                                                                      taacccTCTCCCCAGAGGGGCG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt             >  2:379286/4‑148 (MQ=255)
                                                                                                                                      taacccTCTCCCCAGAGGGGCG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa            >  3:722052/4‑149 (MQ=255)
                                                                                                                                       tCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa            >  1:908689/1‑148 (MQ=255)
                                                                                                                                        acccTCTCCCCAGAGGGGCG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc           >  2:587681/2‑148 (MQ=255)
                                                                                                                                         cccTCTCCCCAGAGGGGCG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                <  4:821723/142‑1 (MQ=255)
                                                                                                                                          ccTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                        >  3:75240/1‑101 (MQ=255)
                                                                                                                                          ccTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                        <  4:75240/101‑1 (MQ=255)
                                                                                                                                          ccTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa         >  1:907364/1‑148 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa       <  3:667718/149‑1 (MQ=255)
                                                                                                                                             cGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc      >  1:387282/1‑148 (MQ=255)
                                                                                                                                               ccccTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAg                                                       >  4:379707/1‑97 (MQ=255)
                                                                                                                                               ccccTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAg                                                       <  3:379707/97‑1 (MQ=255)
                                                                                                                                                cccTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                >  3:119342/1‑119 (MQ=255)
                                                                                                                                                cccTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                <  4:119342/119‑1 (MQ=255)
                                                                                                                                                cccTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                       >  4:615293/1‑128 (MQ=255)
                                                                                                                                                cccTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                       <  3:615293/128‑1 (MQ=255)
                                                                                                                                                cccTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg  >  4:173713/1‑149 (MQ=255)
                                                                                                                                                 ccTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg  <  2:907364/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAG  >  NC_000913/3608548‑3608839

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