Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,686,141 T→C L103L (TTG→CTG)  hmp → nitric oxide dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,686,1410TC99.7% 1231.8 / NA 294L103L (TTG→CTG) hmpnitric oxide dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (129/164);  total (129/165)
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

Only 100 of 302 total aligned reads displayed.

ACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATCGCGA  >  NC_000913/2686043‑2686238
                                                                                                  |                                                                                                 
  ctgcctgcGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGc                                                                                               >  1:4072198/1‑101 (MQ=255)
   tgcctgcGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCa                                                                                              >  2:4072199/1‑101 (MQ=255)
    gcctgcGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCaa                                                                                             <  2:4071630/101‑1 (MQ=255)
    gcctgcGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCaa                                                                                             >  2:4072204/1‑101 (MQ=255)
     cctgcGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAAc                                                                                            >  1:4072205/1‑101 (MQ=255)
       tgcGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCa                                                                                              >  2:4072210/1‑97 (MQ=255)
       tgcGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGc                                                                                          <  1:4071654/101‑1 (MQ=255)
        gcgcTGCTGCCAGCGGTAGAAAAAAGCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCt                                                                                         <  1:4071702/101‑1 (MQ=255)
          gctgctGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGc                                                                                               >  1:4072213/1‑93 (MQ=255)
           ctgctgCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGa                                                                                      <  1:4071750/101‑1 (MQ=255)
            tgctgcCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGAc                                                                                     >  2:4072215/1‑101 (MQ=255)
            tgctgcCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGAc                                                                                     >  1:4072217/1‑101 (MQ=255)
                gcCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGaaa                                                                                 >  2:4072221/1‑101 (MQ=255)
                  cAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATg                                                                               <  2:4071763/101‑1 (MQ=255)
                  cAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATg                                                                               >  1:4072225/1‑101 (MQ=255)
                   aGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGt                                                                              >  2:4072227/1‑101 (MQ=255)
                      ggTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCa                                                                           >  1:4072231/1‑101 (MQ=255)
                        tAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGc                                                                         >  2:4072235/1‑101 (MQ=255)
                          gAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGccc                                                                       <  1:4072030/101‑1 (MQ=255)
                          gAAAAAATAGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGccc                                                                       >  2:4072238/1‑101 (MQ=255)
                           aaaaaaTCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCg                                                                      >  1:4072239/1‑101 (MQ=255)
                             aaaaTCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCggg                                                                    <  1:4071965/101‑1 (MQ=255)
                              aaaTCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCgggg                                                                   >  1:4072242/1‑101 (MQ=255)
                                aTCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCa                                                                 <  2:4071845/101‑1 (MQ=255)
                                   gcgcAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCgg                                                                     >  2:4072248/1‑94 (MQ=255)
                                   gcgcAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGa                                                              >  1:4072247/1‑101 (MQ=255)
                                    cgcAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGaa                                                             <  2:4071482/101‑1 (MQ=255)
                                    cgcAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGaa                                                             <  1:4071795/101‑1 (MQ=255)
                                       aGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACTAAATGTTCAGCCCGGGGCAGGAAGTg                                                          <  1:4071640/101‑1 (MQ=255)
                                       aGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTg                                                          <  2:4071692/101‑1 (MQ=255)
                                         aaGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCt                                                        >  1:4072254/1‑101 (MQ=255)
                                         aaGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCt                                                        >  2:4072253/1‑101 (MQ=255)
                                          aGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTg                                                       <  1:4071862/101‑1 (MQ=255)
                                          aGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTg                                                       <  1:4071899/101‑1 (MQ=255)
                                          aGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTg                                                       <  1:4071812/101‑1 (MQ=255)
                                            cacaCCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGa                                                     >  2:4072258/1‑101 (MQ=255)
                                             acacCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGAc                                                    >  1:4072260/1‑101 (MQ=255)
                                             acacCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGAc                                                    <  1:4071743/101‑1 (MQ=255)
                                             acacCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGAc                                                    <  1:4071777/101‑1 (MQ=255)
                                             acacCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGAc                                                    <  1:4071844/101‑1 (MQ=255)
                                               acCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGAcgc                                                  >  1:4072264/1‑101 (MQ=255)
                                                ccAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTgg                                                      <  2:4071714/96‑1 (MQ=255)
                                                ccAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGAcgcg                                                 >  2:4072266/1‑101 (MQ=255)
                                                 cAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGt                                                <  2:4071972/101‑1 (MQ=255)
                                                  aGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTg                                               <  1:4071417/101‑1 (MQ=255)
                                                   gCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTgg                                              <  2:4071657/101‑1 (MQ=255)
                                                   gCTTCCAGATCAAACCGGAACAGTACAACATCGACGGTGAACACCTGCTGGCa                                                                                              >  2:4072270/1‑53 (MQ=255)
                                                    cTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTggg                                             <  1:4071839/101‑1 (MQ=255)
                                                     ttCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACACTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTgggg                                            >  2:4072274/1‑101 (MQ=255)
                                                      tCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGt                                           <  1:4072139/101‑1 (MQ=255)
                                                       ccAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTa                                          >  1:4072275/1‑101 (MQ=255)
                                                       ccAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTa                                          <  2:4071908/101‑1 (MQ=255)
                                                        cAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTaa                                         >  2:4072278/1‑101 (MQ=255)
                                                        cAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTaa                                         <  2:4072134/101‑1 (MQ=255)
                                                         aGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTaaa                                        <  1:4071875/101‑1 (MQ=255)
                                                         aGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGCGCTGGACGCGTgggg                                            >  2:4072281/1‑97 (MQ=255)
                                                           aTCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGc                                      <  1:4071669/101‑1 (MQ=255)
                                                           aTCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGc                                      <  1:4071738/101‑1 (MQ=255)
                                                           aTCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGc                                      <  2:4071924/101‑1 (MQ=255)
                                                             cAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAg                                       <  2:4071539/98‑1 (MQ=255)
                                                             cAAACCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCt                                    >  2:4072285/1‑101 (MQ=255)
                                                              aaaCCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTa                                   <  2:4071799/101‑1 (MQ=255)
                                                              aaaCCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTa                                   <  1:4071998/101‑1 (MQ=255)
                                                                aCCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATg                                 <  2:4072028/101‑1 (MQ=255)
                                                                aCCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATg                                 <  2:4071854/101‑1 (MQ=255)
                                                                aCCGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATg                                 <  2:4071728/101‑1 (MQ=255)
                                                                 ccGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTgg                                                                                       >  2:4072288/1‑46 (MQ=255)
                                                                  cGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGgt                               <  1:4071944/101‑1 (MQ=255)
                                                                  cGGAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGgt                               <  2:4071957/101‑1 (MQ=255)
                                                                   ggAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTgg                                              >  2:4072293/1‑85 (MQ=255)
                                                                   ggAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGgtg                              <  2:4071737/101‑1 (MQ=255)
                                                                   ggAACAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGgtg                              >  1:4072291/1‑101 (MQ=255)
                                                                     aaCAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTa                            <  2:4072016/101‑1 (MQ=255)
                                                                       cAGTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACt                          >  2:4072296/1‑101 (MQ=255)
                                                                         gTACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTgg                        <  2:4071874/101‑1 (MQ=255)
                                                                          tACAACATCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGc                       <  1:4071857/101‑1 (MQ=255)
                                                                           acaacaTCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAg                                                            <  2:4071791/63‑1 (MQ=255)
                                                                           acaacaTCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCt                      <  2:4071787/101‑1 (MQ=255)
                                                                            caacaTCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTa                     <  1:4071832/101‑1 (MQ=255)
                                                                            caacaTCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTa                     >  1:4072298/1‑101 (MQ=255)
                                                                              acaTCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAAt                   <  2:4071886/101‑1 (MQ=255)
                                                                               caTCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATg                  >  1:4072302/1‑101 (MQ=255)
                                                                                aTCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGt                 <  1:4071765/101‑1 (MQ=255)
                                                                                aTCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGt                 >  1:4072306/1‑101 (MQ=255)
                                                                                aTCGTCGGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGt                 >  2:4072307/1‑101 (MQ=255)
                                                                                 tcgtcgGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGCGGGGTAAAGCCTATGGTGTACTGGCTAATGTa                <  2:4072190/101‑1 (MQ=255)
                                                                                  cgtcgGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTAt               >  2:4072310/1‑101 (MQ=255)
                                                                                     cgGTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTa            >  2:4072314/1‑101 (MQ=255)
                                                                                      ggTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTAt           >  1:4072316/1‑101 (MQ=255)
                                                                                      ggTGAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTAt           <  2:4071902/101‑1 (MQ=255)
                                                                                         gAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCaa        >  2:4072320/1‑101 (MQ=255)
                                                                                         gAACACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCaa        <  2:4072071/101‑1 (MQ=255)
                                                                                          aaCACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAAt       <  2:4071986/101‑1 (MQ=255)
                                                                                          aaCACCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAAt       <  1:4071918/101‑1 (MQ=255)
                                                                                           acacCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATc      >  1:4072323/1‑101 (MQ=255)
                                                                                            cacCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATcg     <  2:4071864/101‑1 (MQ=255)
                                                                                            cacCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATcg     <  2:4071775/101‑1 (MQ=255)
                                                                                             acCTGCTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATcgc    <  2:4071831/101‑1 (MQ=255)
                                                                                               ctgctgGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATCGCGa  <  1:4071822/101‑1 (MQ=255)
                                                                                               ctgctgGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATCGCGa  <  1:4072059/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
ACCTGCCTGCGCTGCTGCCAGCGGTAGAAAAAATCGCGCAGAAGCACACCAGCTTCCAGATCAAACCGGAACAGTACAACATCGTCGGTGAACACCTGTTGGCAACGCTGGACGAAATGTTCAGCCCGGGGCAGGAAGTGCTGGACGCGTGGGGTAAAGCCTATGGTGTACTGGCTAATGTATTTATCAATCGCGA  >  NC_000913/2686043‑2686238

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