Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 816,956 C→A 5.4% intergenic (‑309/‑88) ybhK ← / → moaA putative CofD superfamily transferase/molybdopterin biosynthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913816,9560CA 5.4% 352.4 / 5.0 102intergenic (‑309/‑88)ybhK/moaAputative CofD superfamily transferase/molybdopterin biosynthesis protein A
Reads supporting (aligned to +/- strand):  ref base C (51/45);  new base A (2/4);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

TTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGCCTGTCGAT  >  NC_000913/816810‑817102
                                                                                                                                                  |                                                                                                                                                  
ttCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGACTCTGCACCTGGGTCAACTGATActg                                                                                                                                                  >  6:303985/1‑147 (MQ=255)
   ccTTATTCTCTCTCGAAATAGCACGATCATGACGCTATATTCATGATTACATAGAGAAATTGTGGATTTTAAAAATCTCATTTACACGCTATTATCGGAATAACCACTAAACACGCTCGCCTCTGCACCTGGTTCAAATGATACGGTGc                                                                                                                                               <  6:315913/149‑1 (MQ=255)
       aTTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTg                                                                                                                                           >  7:239801/1‑149 (MQ=255)
        ttCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTgg                                                                                                                                          <  8:211405/149‑1 (MQ=255)
           tGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCg                                                                                                                                       >  4:240049/1‑149 (MQ=255)
             tAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGt                                                                                                                                      <  3:338158/148‑1 (MQ=255)
             tAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGt                                                                                                                                      <  3:240049/148‑1 (MQ=255)
               gCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGa                                                                                                                                    >  7:266672/1‑148 (MQ=255)
                    aTAGCACGATAATGACGCTATATACAATATTACATAGCGACAGTGTCGATTGTAAAAATCTCATTTAAACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTATGGCCGTGACAATGc                                                                                                                              <  5:295335/149‑1 (MQ=255)
                          cGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTa                                                                                                                         >  4:74404/1‑148 (MQ=255)
                                        tataCATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAg                                                                                                           <  8:239801/148‑1 (MQ=255)
                                          taCATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTgg                                                                                                                                          <  3:78806/115‑1 (MQ=255)
                                          taCATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGAAACGGTGCTTTgg                                                                                                                                          >  4:78806/1‑115 (MQ=255)
                                           aCATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGc                                                                                                        >  5:158647/1‑148 (MQ=255)
                                             aTGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCgaa                                                                                                     <  2:338797/149‑1 (MQ=255)
                                             aTGATTACATAGCGAAAGTGAGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGAACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCAACAGGGCg                                                                                                       >  2:101529/1‑147 (MQ=255)
                                              tGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCgaag                                                                                                    <  6:158647/149‑1 (MQ=255)
                                              tGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCgaag                                                                                                    <  7:188825/149‑1 (MQ=255)
                                               gATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCgaag                                                                                                    >  2:192969/1‑148 (MQ=255)
                                                    cATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGgaaga                                                                                               >  8:147426/1‑148 (MQ=255)
                                                     aTAGCGAAAGTGTGGATTGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGaaa                                                                                             <  1:192969/149‑1 (MQ=255)
                                                       aGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAgg                                                                                                          <  4:310776/134‑1 (MQ=255)
                                                       aGCGAAAGTGTGGATGGTAAAAATCTCAATTAAACGCGAGTATCTGCATAACCACAAAACACTCTAGCCTCTGCAACTGGGTCAACTGATACGGTGCTTTGGCAGTGACAATGCTCGTAAAGCTAGCGACAAGGGCGAAGGAAGAAAt                                                                                            >  5:95474/1‑148 (MQ=255)
                                                         cGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGAc                                                                                         >  8:395161/1‑149 (MQ=255)
                                                         cGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGAc                                                                                         >  4:101266/1‑149 (MQ=255)
                                                         cGAAAGTGTGGATGGTAAAAATCGCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGGCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGAc                                                                                         <  1:139624/149‑1 (MQ=255)
                                                            aaGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTc                                                                                      >  8:325658/1‑149 (MQ=255)
                                                                   gATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTccc                                                                               >  4:370396/1‑149 (MQ=255)
                                                                   gATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTccc                                                                               <  5:266427/149‑1 (MQ=255)
                                                                       gTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTAt                                                                           >  1:303395/1‑149 (MQ=255)
                                                                         aaaaaTCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTc                                                                                      >  4:221626/1‑136 (MQ=255)
                                                                         aaaaaTCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTc                                                                                      <  3:221626/136‑1 (MQ=255)
                                                                          aaaaTCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCt                                                                         <  8:266672/148‑1 (MQ=255)
                                                                               ctcATATACACGCTAGTATCGGCAAAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTCATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGACACCAGGGCGAAGGAAGAAATGACTTCGCCTCCAGTATCTGGAAAg                                                                   >  8:233015/1‑149 (MQ=255)
                                                                                  aTTTACCCGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAg                                                                   <  4:96271/146‑1 (MQ=255)
                                                                                  aTTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAg                                                                   <  3:370396/146‑1 (MQ=255)
                                                                                  aTTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAg                                                                   <  3:74404/146‑1 (MQ=255)
                                                                                  aTTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGgtg                                                                >  2:167046/1‑149 (MQ=255)
                                                                                  aTTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGgt                                                                 <  1:135320/148‑1 (MQ=255)
                                                                                  aTTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGgt                                                                 <  7:147426/148‑1 (MQ=255)
                                                                                      acacGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTc                                                                                                                            <  6:383423/85‑1 (MQ=255)
                                                                                      acacGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTc                                                                                                                            >  5:383423/1‑85 (MQ=255)
                                                                                        acGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTccc                                                                               >  5:356245/1‑128 (MQ=255)
                                                                                        acGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATg                                                          >  3:152857/1‑149 (MQ=255)
                                                                                        acGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATg                                                          >  6:144387/1‑149 (MQ=255)
                                                                                          gCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGc                                                        <  6:389379/149‑1 (MQ=255)
                                                                                          gCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGc                                                        >  7:270270/1‑149 (MQ=255)
                                                                                              gaaTCTGCATAACCACTATACACTCTAGCCGCTGGACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTGCCGAATCTGGAAAGGGGTACATGGCtt                                                      <  6:95474/145‑1 (MQ=255)
                                                                                              gTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCtt                                                      <  2:303395/147‑1 (MQ=255)
                                                                                                 tCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACaa                                                 >  1:354518/1‑149 (MQ=255)
                                                                                                      aTAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACGTCGCCTCCCGTATCTGGAAAGTTGTACATGGCTTCACAactgac                                            >  4:264854/1‑149 (MQ=255)
                                                                                                       tAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAactga                                             <  7:383603/147‑1 (MQ=255)
                                                                                                       tAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAactga                                             >  8:383603/1‑147 (MQ=255)
                                                                                                          ccACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGAt                                        >  3:295895/1‑149 (MQ=255)
                                                                                                            aCTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATg                                       >  6:125238/1‑148 (MQ=255)
                                                                                                            aCTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGc                                      <  7:325658/149‑1 (MQ=255)
                                                                                                            aCTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGc                                      <  7:82538/149‑1 (MQ=255)
                                                                                                             cTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTAGCTGGAAAGGTGTACATGGCTTCACAACTGACTGCTGCa                                     >  6:41366/1‑149 (MQ=255)
                                                                                                                aaCACTCTAGCCTCTGCCCCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCCCCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCAttt                                  >  5:32681/1‑149 (MQ=255)
                                                                                                                aaCACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGAt                                        >  3:204919/1‑143 (MQ=255)
                                                                                                                aaCACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCAttt                                  <  4:295895/149‑1 (MQ=255)
                                                                                                                aaCACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCAt                                    <  7:395161/147‑1 (MQ=255)
                                                                                                                 acacTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTg                                 >  4:238488/1‑149 (MQ=255)
                                                                                                                  cacTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAg                                                                   <  1:84450/114‑1 (MQ=255)
                                                                                                                  cacTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAg                                                                   >  2:84450/1‑114 (MQ=255)
                                                                                                                   acTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCACGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTgcg                               >  1:271078/1‑149 (MQ=255)
                                                                                                                    ctctAGCCTCTGCACCTGGGTCAACTGATTCGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTgcgc                              >  2:311194/1‑149 (MQ=255)
                                                                                                                     tctAGCCTCTGCACCTGGTTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCAttt                                  <  1:329031/144‑1 (MQ=255)
                                                                                                                     tctAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATTCTCGTAAAGATTGCACCCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGcgcg                             >  1:364149/1‑149 (MQ=255)
                                                                                                                     tctAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTcac                                                   >  8:243904/1‑127 (MQ=255)
                                                                                                                     tctAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTcac                                                   <  7:243904/127‑1 (MQ=255)
                                                                                                                     tctAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGcgcg                             >  8:343835/1‑149 (MQ=255)
                                                                                                                     tctAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGcgcg                             >  6:282453/1‑149 (MQ=255)
                                                                                                                     tctAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGACTACATGGCTTCACAACTGACTGATTCAttagatat                             >  2:329031/1‑143 (MQ=255)
                                                                                                                       tAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGc                                                                                    <  6:43948/92‑1 (MQ=255)
                                                                                                                       tAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGc                                                                                    >  5:43948/1‑92 (MQ=255)
                                                                                                                       tAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGt                            >  1:51221/1‑148 (MQ=255)
                                                                                                                          ccTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTc                                                                                 >  6:92568/1‑92 (MQ=255)
                                                                                                                          ccTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTc                                                                                 <  5:92568/92‑1 (MQ=255)
                                                                                                                            tctGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGttt                      >  1:219339/1‑149 (MQ=255)
                                                                                                                              tGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAATAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGttt                      <  1:101529/147‑1 (MQ=255)
                                                                                                                              tGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGttt                      <  7:343835/147‑1 (MQ=255)
                                                                                                                              tGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGttt                      <  5:282453/147‑1 (MQ=255)
                                                                                                                               gCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTa                                                                            <  8:304866/92‑1 (MQ=255)
                                                                                                                               gCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTa                                                                            >  7:304866/1‑92 (MQ=255)
                                                                                                                                cACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGATTGATGCATTTGCGCGTAAGTTTta                    <  7:233015/147‑1 (MQ=255)
                                                                                                                                cACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTtac                   >  7:331522/1‑148 (MQ=255)
                                                                                                                                cACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTta                    <  1:167046/147‑1 (MQ=255)
                                                                                                                                cACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTta                    <  6:32681/147‑1 (MQ=255)
                                                                                                                                  ccTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTtactac                <  2:51221/149‑1 (MQ=255)
                                                                                                                                   cTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTactact               >  6:344775/1‑149 (MQ=255)
                                                                                                                                   cTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTactact               <  5:303985/149‑1 (MQ=255)
                                                                                                                                       gTCAACTGATAAGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCAt                                    <  1:212138/124‑1 (MQ=255)
                                                                                                                                       gTCAACTGATAAGCTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCAt                                    >  2:212138/1‑124 (MQ=255)
                                                                                                                                        tCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTgcgc          <  6:209890/149‑1 (MQ=255)
                                                                                                                                            cTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGATAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGCCTGt      >  8:119493/1‑149 (MQ=255)
                                                                                                                                            cTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGCCTGt      >  6:131093/1‑149 (MQ=255)
                                                                                                                                               aTACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTtac                   >  5:318835/1‑133 (MQ=255)
                                                                                                                                               aTACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTtac                   <  6:318835/133‑1 (MQ=255)
                                                                                                                                                tACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGCCTGTCGAt  >  8:273948/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTCCCTTATTCTGTAGCGAAATAGCACGATCATGACGCTATATACATGATTACATAGCGAAAGTGTGGATGGTAAAAATCTCATTTACACGCTAGTATCGGCATAACCACTAAACACTCTAGCCTCTGCACCTGGGTCAACTGATACGGTGCTTTGGCCGTGACAATGCTCGTAAAGATTGCCACCAGGGCGAAGGAAGAAATGACTTCGCCTCCCGTATCTGGAAAGGTGTACATGGCTTCACAACTGACTGATGCATTTGCGCGTAAGTTTTACTACTTGCGCCTGTCGAT  >  NC_000913/816810‑817102

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