Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,701 A→C 8.7% intergenic (+52/+50) 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,7010AC 8.7% 211.6 / 16.1 82intergenic (+52/+50)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (34/37);  new base C (4/3);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-01

GTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATAT  >  NC_000913/3608553‑3608845
                                                                                                                                                     |                                                                                                                                                 
gTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                   <  3:265978/149‑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: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                                                                                                                                  <  3:309811/149‑1 (MQ=255)
                ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAgg                                                                                                                                  <  4:785094/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                                                                                                                          <  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                                                                                                                             <  1:910059/52‑1 (MQ=255)
                                                                                                                       aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGaaa                                                                                                                             >  2:910059/1‑52 (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)
                                                                                                                                           ccctggggggAG‑AGGGTTGGGGAGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa         <  1:529593/142‑1 (MQ=255)
                                                                                                                                           cccagaggggCG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa         <  1:379286/142‑1 (MQ=255)
                                                                                                                                           cccTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                      <  4:119342/119‑1 (MQ=255)
                                                                                                                                           cccTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                      >  3:119342/1‑119 (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)
                                                                                                                                               ggggggAG‑GGGGTTGGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa     <  2:458674/146‑1 (MQ=255)
                                                                                                                                                  gggAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGTAatat  <  2:648394/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACC‑GATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAG‑AGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATAT  >  NC_000913/3608553‑3608845

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