Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,449,803 G→C 7.5% A516G (GCT→GGT) ‡ tynA ← tyramine oxidase, copper‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,449,8030GC 7.5% 324.1 / 11.5 102A516G (GCT→GGT) ‡tynAtyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base G (42/52);  new base C (2/6);  total (44/58)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 105 total aligned reads displayed.

CATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATACTCC  >  NC_000913/1449654‑1449952
                                                                                                                                                     |                                                                                                                                                     
cATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAg                                                                                                                                                       <  2:1349701/150‑1 (MQ=255)
cATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAg                                                                                                                                                       >  2:1280293/1‑150 (MQ=255)
 aTCCAGATCGAGGCGGAAATTATAAATATGTTGGTGGGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGc                                                                                                                                                      <  2:1713693/150‑1 (MQ=255)
    cAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCaac                                                                                                                                                   >  2:143741/1‑148 (MQ=255)
    cAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGAGTGTCATCTTTCGCAGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTAACCGCTTCGATGCCCGTAGCAcc                                                                                                                                                   >  1:1188855/1‑150 (MQ=255)
     aGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCg                                                                                                                                                  <  2:938984/150‑1 (MQ=255)
     aGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCg                                                                                                                                                  >  2:1735220/1‑150 (MQ=255)
        tCGAGTCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCa                                                                                                                                               <  2:734184/150‑1 (MQ=255)
        tCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCa                                                                                                                                               <  1:1075903/150‑1 (MQ=255)
          gaggagGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTGCACCGCTTCCATGCCCGTAGCCCCGGCATc                                                                                                                                             >  2:227299/1‑150 (MQ=255)
            ggCGGAAATTATAAATATCTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGa                                                                                                                                           <  2:1564103/150‑1 (MQ=255)
               gtaaaTTATAAATATGTTCGTGTGTAGTCCCCACGATATTGTGATGTTTAAGCGTTCCGTAGCGCGTGTCATCTTTCGACGTCTCATCGCGCATGGTTTACGCGTTAACACGTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGc                                                                                                                                        <  2:165293/148‑1 (MQ=255)
               ggAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGAAGc                                                                                                                                        <  1:561563/150‑1 (MQ=255)
               ggAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGAAGc                                                                                                                                        <  1:1324659/150‑1 (MQ=255)
               ggAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTAAACACCTTTAAACGCTTCGATGCCCGTAGCACCGGCATCGATGc                                                                                                                                        >  2:335955/1‑150 (MQ=255)
               ggAAATAATAAATATGTTTGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTTATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTGCGATGCCGGTGCTACGGGCATCGAAGc                                                                                                                                        <  2:1172531/150‑1 (MQ=255)
                  aaTTATAAATATGTTGGTGTTTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCATCGATGCCGGTGCTACGGGCATCGAAGcggt                                                                                                                                     <  1:146007/150‑4 (MQ=255)
                  aaTTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGAAGcggt                                                                                                                                     <  2:526944/150‑4 (MQ=255)
                     tataAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAg                                                                                                                                  >  2:256128/1‑150 (MQ=255)
                       taAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTg                                                                                                                                >  1:478146/1‑150 (MQ=255)
                       taAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTg                                                                                                                                >  1:1801283/1‑150 (MQ=255)
                        aaaTATGTTGGTGTTTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGc                                                                                                                               <  1:907289/150‑1 (MQ=255)
                        aaaTATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGc                                                                                                                               <  2:653645/150‑1 (MQ=255)
                         aaTATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGcc                                                                                                                              <  1:1678781/150‑1 (MQ=255)
                           tatGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGt                                                                                                                            <  1:1646230/150‑1 (MQ=255)
                                 ggggtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCAt                                                                                                                      <  1:994096/147‑1 (MQ=255)
                                 ggTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGGCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCAt                                                                                                                      <  2:300056/150‑1 (MQ=255)
                                  gtgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATg                                                                                                                     <  1:767028/150‑1 (MQ=255)
                                   tgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATgg                                                                                                                    >  2:1800244/1‑150 (MQ=255)
                                           aCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCa                                                                                                            >  2:840486/1‑150 (MQ=255)
                                           aCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTGTCGCCGTCTCATCGTGCATGGTTTTCGCTTTACCACCTTTCCCCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATccc                                                                                                            >  2:1188382/1‑149 (MQ=255)
                                              cACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTTTCATCTTTCTCCGTCTCATCGTGCATTGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTc                                                                                                         <  1:1447424/150‑1 (MQ=255)
                                              cACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTc                                                                                                         >  1:487668/1‑150 (MQ=255)
                                               aCGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCa                                                                                                        <  2:247021/150‑1 (MQ=255)
                                                  atatTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCaaaa                                                                                                     >  1:172571/1‑150 (MQ=255)
                                                   tatTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTTCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCaataa                                                                                                    <  2:277624/150‑1 (MQ=255)
                                                     ttGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATg                                                                                                  >  1:780830/1‑150 (MQ=255)
                                                     ttGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATg                                                                                                  <  2:1021137/150‑1 (MQ=255)
                                                       gtgATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGTTGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTa                                                                                                <  2:1598594/150‑1 (MQ=255)
                                                       gtgATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTa                                                                                                >  2:936013/1‑150 (MQ=255)
                                                        tgATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCGCATAGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAg                                                                                               <  1:1159321/150‑1 (MQ=255)
                                                        tgATCGATAAGCGTGCCGTAGCGCGGGTCAGCTTTCGCCGGCTCATCGTGCATGGTTTTCGCTTTATCACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAACTGTAg                                                                                               <  1:1185245/150‑1 (MQ=255)
                                                         gATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGt                                                                                              >  2:1244814/1‑150 (MQ=255)
                                                         gATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGt                                                                                              <  2:1238042/150‑1 (MQ=255)
                                                            cGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCAt                                                                                           >  2:777091/1‑150 (MQ=255)
                                                             gATTAGCGTGCCGTAGCGCGTGTCTTCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTAATGGAAGATCCAGTCAAAAATGTAGTCATa                                                                                          <  2:1155313/150‑1 (MQ=255)
                                                             gATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATa                                                                                          <  2:750727/150‑1 (MQ=255)
                                                                 aGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTa                                                                                      <  1:826327/150‑1 (MQ=255)
                                                                  gCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTAc                                                                                     <  2:1276936/150‑1 (MQ=255)
                                                                  gCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTAc                                                                                     <  1:1735220/150‑1 (MQ=255)
                                                                   ggTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTAcc                                                                                    <  1:1280293/149‑1 (MQ=255)
                                                                    gagCCGTATCGCGTGTCATCTTTCGCCGTCTCATCTTGCATGGTTTTCGCTTTAACACCTTTCACCTCTTCTATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATACAGTCAAAAATGTAGTCATAGTTAccc                                                                                   <  1:33516/148‑1 (MQ=255)
                                                                    gTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTAccc                                                                                   >  2:126347/1‑150 (MQ=255)
                                                                        cGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTg                                                                               >  1:1732139/1‑150 (MQ=255)
                                                                         gTATCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGt                                                                              <  1:1147978/150‑1 (MQ=255)
                                                                         gTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGTAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGt                                                                              <  1:1412845/150‑1 (MQ=255)
                                                                          tAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTa                                                                             >  2:921474/1‑150 (MQ=255)
                                                                           atcgcgTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCATCGATGCCCTTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTAc                                                                            <  1:378565/148‑1 (MQ=255)
                                                                           aGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTAc                                                                            <  1:1773873/150‑1 (MQ=255)
                                                                           aGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGGTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTAc                                                                            >  1:99096/1‑150 (MQ=255)
                                                                                gtgtCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATc                                                                       >  2:1635646/1‑150 (MQ=255)
                                                                                  gtCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCa                                                                     >  1:1788177/1‑150 (MQ=255)
                                                                                  gtCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCa                                                                     <  2:613427/150‑1 (MQ=255)
                                                                                   tCATCTTTCGCCGTCCCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAg                                                                    <  1:224720/150‑1 (MQ=255)
                                                                                      tCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAgcgc                                                                 >  2:794001/1‑150 (MQ=255)
                                                                                      tCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAgcgc                                                                 >  2:153342/1‑150 (MQ=255)
                                                                                       cTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGca                                                                <  1:1645957/150‑1 (MQ=255)
                                                                                         ttCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGcacc                                                              <  1:1544615/150‑1 (MQ=255)
                                                                                            gCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCCCt                                                           >  1:1114666/1‑150 (MQ=255)
                                                                                             ccGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGCTGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTa                                                          >  1:1721835/1‑150 (MQ=255)
                                                                                              cGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCACGTAGCACCGGCATCTATGCCAAAAGTCAATTTATCATGGAAGATCCAGAAAAAAATGTAGTCATAGTTACCCACTGTACTGATCCATCGCACCACTaa                                                         >  1:1373082/1‑150 (MQ=255)
                                                                                                   cATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCCCCACTAACTccc                                                    >  1:1224691/1‑150 (MQ=255)
                                                                                                      cGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGc                                                 <  2:636690/150‑1 (MQ=255)
                                                                                                        tGCATGGTTTTCGCTTTAACACCTTTCCCCGCTTCGATGCCCGTAGCACCTGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTGATAGTTACCCACTGTACTGATACAGCGCACCACTAACTCCCGGCGt                                               >  1:925339/1‑150 (MQ=255)
                                                                                                          cATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTc                                             >  2:83393/1‑150 (MQ=255)
                                                                                                             ggTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGt                                          >  2:52671/1‑150 (MQ=255)
                                                                                                             ggTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGt                                          <  1:1019842/150‑1 (MQ=255)
                                                                                                             ggTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGt                                          <  2:1801669/150‑1 (MQ=255)
                                                                                                             ggTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGt                                          <  1:256128/150‑1 (MQ=255)
                                                                                                                    gCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCCCCACTAACTCCCGGCGTTCGGTACTGACg                                   >  2:1099471/1‑150 (MQ=255)
                                                                                                                            acCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTg                           >  2:1799467/1‑150 (MQ=255)
                                                                                                                             ccTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATACAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTgg                          >  2:1124625/1‑150 (MQ=255)
                                                                                                                              cTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGc                         <  1:1799467/150‑1 (MQ=255)
                                                                                                                               tttCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCATTAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCCGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGcc                        >  1:1195712/1‑150 (MQ=255)
                                                                                                                               tttCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGcc                        <  1:82051/150‑1 (MQ=255)
                                                                                                                                 tCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGcccc                      <  1:1721469/150‑2 (MQ=255)
                                                                                                                                 tCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCa                      <  2:1211166/150‑1 (MQ=255)
                                                                                                                                     cGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTc                  <  1:1419488/150‑1 (MQ=255)
                                                                                                                                     cGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCCAGCCCaacgt                  <  1:936013/150‑5 (MQ=255)
                                                                                                                                      gcatCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAATATGTAGTCAAAGTTACCCAATGTACTTAGGCAGCGCACCACTACCTCCCGGCGTTCGGTACTGACGTTGGGCAGGCCCATTTca                 >  2:697205/4‑149 (MQ=255)
                                                                                                                                       catCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt                >  1:1191319/3‑150 (MQ=255)
                                                                                                                                          cGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGAt             <  1:1440911/150‑1 (MQ=255)
                                                                                                                                          cGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGAt             <  2:487668/150‑1 (MQ=255)
                                                                                                                                          cGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTAATCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTCCCCACTGTACTGATCCAGCGGACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCACATTTCCttct             >  2:79720/1‑147 (MQ=255)
                                                                                                                                           gATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATg            <  2:1735388/150‑1 (MQ=255)
                                                                                                                                             tGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCt          <  2:780830/150‑1 (MQ=255)
                                                                                                                                               cccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCCGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTta        >  1:903692/1‑150 (MQ=255)
                                                                                                                                               cccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTta        >  1:1357624/1‑150 (MQ=255)
                                                                                                                                               cccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTta        <  2:1011511/150‑1 (MQ=255)
                                                                                                                                                     gCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACAACTAAATCCCGGCTTTCGGTACTGAAGTTGGGCTGGCCCATTTACTGATGCTTATACtaa  >  1:1737551/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATACTCC  >  NC_000913/1449654‑1449952

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