Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,689 G→T 7.0% intergenic (+40/+62) 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,6890GT 7.0% 242.1 / 5.1 81intergenic (+40/+62)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base T (4/2);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAC  >  NC_000913/3608547‑3608835
                                                                                                                                               |                                                                                                                                                   
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                                                                                                                                            >  3:130433/1‑145 (MQ=255)
       ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                            <  4:130433/145‑1 (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:238879/149‑1 (MQ=255)
                     gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAg                                                                                                                         <  3:781725/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                                                                                                                <  2:651202/100‑1 (MQ=255)
                                                                               ttAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAgg                                                                                                                >  1:651202/1‑100 (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)
                                                                                                                                               |                                                                                                                                                   
CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAC  >  NC_000913/3608547‑3608835

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: ‑