Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,449,802 A→G 5.4% A516A (GCT→GCC) ‡ tynA ← tyramine oxidase, copper‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,449,8020AG 5.4% 288.8 / 4.9 93A516A (GCT→GCC) ‡tynAtyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base A (44/43);  new base G (3/2);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

CATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATA  >  NC_000913/1449654‑1449948
                                                                                                                                                    |                                                                                                                                                  
cATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTa                                                                                                                                                    >  2:757023/1‑149 (MQ=255)
cATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTa                                                                                                                                                    >  1:116144/1‑149 (MQ=255)
 aTCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAg                                                                                                                                                   >  1:121750/1‑149 (MQ=255)
     agctCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCAcc                                                                                                                                               <  1:36386/146‑1 (MQ=255)
        tCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGc                                                                                                                                            >  1:697967/1‑149 (MQ=255)
          gAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCAc                                                                                                                                                >  2:1370625/1‑143 (MQ=255)
          gAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCAc                                                                                                                                                <  1:1370625/143‑1 (MQ=255)
              cGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGAt                                                                                                                                      <  1:625427/149‑1 (MQ=255)
               ggAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTATAACACCTTTCACCGCTTCGAAGCCCGTCGCACCGGCATCGATg                                                                                                                                     >  2:862225/1‑149 (MQ=255)
                taaaTTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGc                                                                                                                                    <  1:1289635/148‑1 (MQ=255)
                gAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGc                                                                                                                                    >  2:1335048/1‑149 (MQ=255)
                gAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGc                                                                                                                                    <  1:836406/149‑1 (MQ=255)
                gAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGc                                                                                                                                    >  2:1452203/1‑149 (MQ=255)
                gAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGGGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGc                                                                                                                                    <  1:410797/149‑1 (MQ=255)
                      ataAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAg                                                                                                                              <  1:1135163/149‑1 (MQ=255)
                      ataAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAg                                                                                                                              >  1:5977/1‑149 (MQ=255)
                        aaaTATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTg                                                                                                                            >  1:1381731/1‑149 (MQ=255)
                          atatGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGcc                                                                                                                          <  1:47225/149‑1 (MQ=255)
                           tatGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCg                                                                                                                         >  2:1198266/1‑149 (MQ=255)
                                  gtgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCAt                                                                                                                  >  1:112661/1‑149 (MQ=255)
                                       tagtaCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGaag                                                                                                                                     >  1:34870/1‑123 (MQ=255)
                                       tagtaCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGaag                                                                                                                                     <  2:34870/125‑3 (MQ=255)
                                       tagtaCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGCTGCCAATAGTGCCGTTTTCATGTAAg                                                                                                             >  2:870427/1‑149 (MQ=255)
                                       tagtaCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAg                                                                                                             >  1:1115514/1‑149 (MQ=255)
                                       tagtaCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAg                                                                                                             >  1:1322136/1‑149 (MQ=255)
                                          taCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATc                                                                                                          <  2:380953/149‑1 (MQ=255)
                                           aCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATTGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATTGAAGATcc                                                                                                         <  2:922845/149‑1 (MQ=255)
                                            cccACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCa                                                                                                        <  1:757023/149‑1 (MQ=255)
                                             ccACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGAAGc                                                                                                                                    <  1:879781/120‑1 (MQ=255)
                                             ccACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGAAGc                                                                                                                                    >  2:879781/1‑120 (MQ=255)
                                             ccACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAg                                                                                                       <  1:916965/149‑1 (MQ=255)
                                             ccACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGAATGGTTTTCGCTTTAACAACTTTCACCGCTTCGATGCACGTAGCACCGGCATCGATGCCAATACTGCCGTTTTCATGGAAGATCCat                                                                                                       >  2:257606/1‑148 (MQ=255)
                                              cACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGt                                                                                                      >  1:1541918/1‑149 (MQ=255)
                                               aCGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTc                                                                                                     >  1:180303/1‑149 (MQ=255)
                                                    atTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCaaaaa                                                                                                <  1:1335048/149‑1 (MQ=255)
                                                     ttGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAAt                                                                                               <  2:674340/149‑1 (MQ=255)
                                                     ttGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAAt                                                                                               <  2:1592783/149‑1 (MQ=255)
                                                       gtgATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAAAAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGt                                                                                             <  1:710661/149‑1 (MQ=255)
                                                          aTCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGt                                                                                          >  2:1479668/1‑149 (MQ=255)
                                                           tCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCa                                                                                                        >  2:910769/1‑134 (MQ=255)
                                                           tCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCa                                                                                                        <  1:910769/134‑1 (MQ=255)
                                                           tCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTc                                                                                         >  2:14115/1‑149 (MQ=255)
                                                           tCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTc                                                                                         <  2:41374/149‑1 (MQ=255)
                                                             gATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCAt                                                                                       <  2:112661/149‑1 (MQ=255)
                                                             gATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCAt                                                                                       >  1:267760/1‑149 (MQ=255)
                                                              aTAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATa                                                                                      <  1:1198266/149‑1 (MQ=255)
                                                               tAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGGATGGTTTTCGCTTTAACACCGTTCACCGCTTCGATGCCCGTAGCCCCGGCATCGATGCCATGAGTGACGTATTCATGGAAGATCCAGTCACAAATGTCGTCaaac                                                                                     >  1:78105/1‑146 (MQ=255)
                                                               tAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAg                                                                                     <  2:40865/149‑1 (MQ=255)
                                                                   cGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTAc                                                                                 <  2:1243633/149‑1 (MQ=255)
                                                                   cGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTAc                                                                                 <  2:1237522/149‑1 (MQ=255)
                                                                       ccGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCAc                                                                             <  2:190254/149‑1 (MQ=255)
                                                                        cGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATTGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACt                                                                            <  2:1166354/149‑1 (MQ=255)
                                                                             cgcgTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACt                                                                       <  2:1115514/149‑1 (MQ=255)
                                                                              gcgTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTg                                                                      >  2:1427308/1‑149 (MQ=255)
                                                                                 tgtCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATc                                                                   >  1:178849/1‑149 (MQ=255)
                                                                                 tgtCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATc                                                                   <  2:697967/149‑1 (MQ=255)
                                                                                     aTCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAgc                                                               <  1:1423345/149‑1 (MQ=255)
                                                                                       cTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATACAGTCAAAAATGTATTCATAGTTACCCACTGTACTGATCCAgcgc                                                             >  1:646961/1‑149 (MQ=255)
                                                                                        tttCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGca                                                            <  1:862225/149‑1 (MQ=255)
                                                                                        tttCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGca                                                            >  1:187874/1‑149 (MQ=255)
                                                                                        tttCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGca                                                            >  2:374280/1‑149 (MQ=255)
                                                                                            gCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGcaccac                                                        <  2:125242/149‑1 (MQ=255)
                                                                                             ccGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATTCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACt                                                       <  1:949923/149‑1 (MQ=255)
                                                                                                   cATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTcc                                                 >  2:906827/1‑149 (MQ=255)
                                                                                                    aTCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTccc                                                >  2:699803/1‑149 (MQ=255)
                                                                                                    aTCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAAAGTAGTCATAGTTACCCACTGTACTGATCCAGCGCCCCACTAACTccc                                                >  2:991755/1‑149 (MQ=255)
                                                                                                        tGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCg                                            <  2:121750/149‑1 (MQ=255)
                                                                                                        tGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGACAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCg                                            >  2:180330/1‑149 (MQ=255)
                                                                                                         gCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGt                                           >  1:294172/1‑149 (MQ=255)
                                                                                                               ttttCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCTTCGATGCCAATAGTGACGTTTTCATGGTAGATCCAGTCAAAAATGTAGTCATAGATCACCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGCTa                                     >  1:319321/1‑149 (MQ=255)
                                                                                                               ttttCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTa                                     <  1:1479668/149‑1 (MQ=255)
                                                                                                                tttCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTAc                                    >  1:313681/1‑149 (MQ=255)
                                                                                                                    gCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACGCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGAc                                >  1:626982/1‑149 (MQ=255)
                                                                                                                    gCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGAc                                <  1:1452203/149‑1 (MQ=255)
                                                                                                                      tttAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGt                              <  2:116144/149‑1 (MQ=255)
                                                                                                                          acacCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTggg                          >  1:1299580/1‑149 (MQ=255)
                                                                                                                           cacCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGc                         >  1:1028989/1‑149 (MQ=255)
                                                                                                                            acCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCt                        >  1:1587836/1‑149 (MQ=255)
                                                                                                                              cTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGTGCTgg                      <  2:5977/149‑1 (MQ=255)
                                                                                                                               tttCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACACACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGc                     >  1:254178/1‑149 (MQ=255)
                                                                                                                                  cACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCa                  <  2:1541918/149‑1 (MQ=255)
                                                                                                                                   ttggcatCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCAt                 >  1:475694/7‑149 (MQ=255)
                                                                                                                                   aCCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCAt                 >  1:505147/1‑149 (MQ=255)
                                                                                                                                   aCCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCAt                 >  2:242567/1‑149 (MQ=255)
                                                                                                                                    ccGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCAtt                >  2:173871/1‑149 (MQ=255)
                                                                                                                                    ccGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCAtt                >  2:144786/1‑149 (MQ=255)
                                                                                                                                    ccGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCAtt                <  2:1333131/149‑1 (MQ=255)
                                                                                                                                        ttCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt            <  2:1262150/149‑1 (MQ=255)
                                                                                                                                          cGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGa          <  2:1572929/149‑1 (MQ=255)
                                                                                                                                          cGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGa          <  1:16640/149‑1 (MQ=255)
                                                                                                                                           gATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGAt         >  2:201871/1‑149 (MQ=255)
                                                                                                                                              gcacggaGCACCGGCATGGATGCCAATAGTGCCGGTTTGATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCCGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCt      <  2:78105/143‑1 (MQ=255)
                                                                                                                                               cccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCTCTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCtt     <  2:1322136/149‑1 (MQ=255)
                                                                                                                                                ccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTACCTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTta    >  2:316701/1‑149 (MQ=255)
                                                                                                                                                ccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTta    <  2:267760/149‑1 (MQ=255)
                                                                                                                                                  gTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTtata  >  1:590887/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTATA  >  NC_000913/1449654‑1449948

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: