Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,694 T→A 7.3% intergenic (+45/+57) 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,6940TA 7.3% 250.4 / 9.3 86intergenic (+45/+57)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (37/40);  new base A (4/2);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAG  >  NC_000913/3608547‑3608839
                                                                                                                                                    |                                                                                                                                                  
cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCct                                                                                                                                                    >  2:588417/1‑148 (MQ=255)
 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:112758/1‑149 (MQ=255)
                aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAgg                                                                                                                                  <  4:470791/149‑1 (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                                                                                                                                  >  4:296398/1‑135 (MQ=255)
                              ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAgg                                                                                                                                  <  3:296398/135‑1 (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                                                       <  3:379707/97‑1 (MQ=255)
                                                                                                                                                ccccTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAg                                                       >  4:379707/1‑97 (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)
                                                                                                                                                    |                                                                                                                                                  
CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAG  >  NC_000913/3608547‑3608839

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