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

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

ACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTAT  >  NC_000913/1449653‑1449947
                                                                                                                                                     |                                                                                                                                                 
aCATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCCTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTa                                                                                                                                                   <  1:130300/150‑1 (MQ=255)
aCATCCAGATCGAGGCGGAAATTATAAATATGTTGGGGTGTAGTACCCACGATATTGTGATCGATAAGCGGGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTa                                                                                                                                                   <  2:982300/150‑1 (MQ=255)
   tCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCa                                                                                                                                                >  1:1032439/1‑150 (MQ=255)
   tCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCa                                                                                                                                                >  2:1117970/1‑150 (MQ=255)
    ccAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCACGTAGCAc                                                                                                                                               >  1:306897/1‑150 (MQ=255)
      aGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTTATCGATAAGCGTGCCTTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCg                                                                                                                                             >  2:1330440/1‑150 (MQ=255)
      aGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCg                                                                                                                                             >  2:1630153/1‑150 (MQ=255)
               cGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGaag                                                                                                                                    <  2:216748/150‑3 (MQ=255)
               cGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATg                                                                                                                                    >  2:1266014/1‑150 (MQ=255)
               cGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATg                                                                                                                                    >  1:1526843/1‑150 (MQ=255)
                ggAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAATCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGAAGc                                                                                                                                   <  2:964933/150‑1 (MQ=255)
                 gAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGcc                                                                                                                                  <  1:279263/150‑1 (MQ=255)
                  aaaTTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCa                                                                                                                                 >  1:325858/1‑150 (MQ=255)
                  aaaTTATAAATATGTTGATGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCGGTGCTACGGGCATCGAAGcgg                                                                                                                                 <  1:590598/150‑3 (MQ=255)
                   aaTTATAAATCTGTTGGGGTGTAGTACCCACGATATTGTGATCGATAAGCGGGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCaa                                                                                                                                <  1:682562/150‑1 (MQ=255)
                        taAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTg                                                                                                                           >  1:411538/1‑150 (MQ=255)
                          aaTATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGcc                                                                                                                         >  1:404858/1‑150 (MQ=255)
                           atatGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCg                                                                                                                        >  2:1464740/1‑150 (MQ=255)
                            tatGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGt                                                                                                                       <  2:985831/150‑1 (MQ=255)
                            tatGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGt                                                                                                                       <  2:272168/150‑1 (MQ=255)
                              tGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGttt                                                                                                                     >  1:269721/1‑150 (MQ=255)
                              tGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGttt                                                                                                                     <  2:845441/150‑1 (MQ=255)
                                 tGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCa                                                                                                                  <  1:1524589/150‑1 (MQ=255)
                                  ggTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCAt                                                                                                                 >  2:1132277/1‑150 (MQ=255)
                                   gtgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATg                                                                                                                <  2:1242479/150‑1 (MQ=255)
                                   gtgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATg                                                                                                                <  1:747633/150‑1 (MQ=255)
                                    tgtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATgg                                                                                                               <  1:1425729/150‑1 (MQ=255)
                                    ggtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATgg                                                                                                               <  2:1260631/149‑1 (MQ=255)
                                     gtgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGa                                                                                                              <  2:306897/150‑1 (MQ=255)
                                     gtgtAGTACCCACGATATTGTGATAGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGa                                                                                                              >  2:83061/1‑150 (MQ=255)
                                      tgtAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGaa                                                                                                             <  1:72610/150‑1 (MQ=255)
                                       gtagtaCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAg                                                                                                            <  1:805071/150‑1 (MQ=255)
                                           taCCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATcc                                                                                                        <  2:259516/150‑1 (MQ=255)
                                           taCCCACGATATTGTGAGCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTGTCGCTTTACCACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATcc                                                                                                        <  2:159657/150‑1 (MQ=255)
                                                  gATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCaaa                                                                                                 >  1:918637/1‑150 (MQ=255)
                                                      ttGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATg                                                                                             <  1:678891/150‑1 (MQ=255)
                                                         tgATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAg                                                                                          <  2:404858/150‑1 (MQ=255)
                                                         tgATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAg                                                                                          <  2:1072757/150‑1 (MQ=255)
                                                           aTCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTc                                                                                        <  2:336448/150‑1 (MQ=255)
                                                           aTCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTc                                                                                        <  2:41164/150‑1 (MQ=255)
                                                              gATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATa                                                                                     >  2:1271406/1‑150 (MQ=255)
                                                              gATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATa                                                                                     <  1:250676/150‑1 (MQ=255)
                                                                tAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGt                                                                                   <  1:83061/150‑1 (MQ=255)
                                                                  aGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTa                                                                                 >  2:730823/1‑150 (MQ=255)
                                                                       gCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCAc                                                                            >  1:1431050/1‑150 (MQ=255)
                                                                       gCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCAc                                                                            <  2:325858/150‑1 (MQ=255)
                                                                       gCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCAc                                                                            >  1:55008/1‑150 (MQ=255)
                                                                            tgcgcGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTAc                                                                       <  1:95738/149‑1 (MQ=255)
                                                                               gcgTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGa                                                                    <  1:230135/150‑1 (MQ=255)
                                                                                 gtgtCATCTTTCGCCGTCTCATCGTGCATGTTTTTCGCTTTAACAGCTTTCACCGCTTCGATGGCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGTTCCAGTCAAATATGTAGTCATAGTTACCCAGTGTACTGAtg                                                                  <  2:968250/150‑2 (MQ=255)
                                                                                 gtgtCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATc                                                                  <  2:1076253/150‑1 (MQ=255)
                                                                                        cTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGca                                                           >  1:1133087/1‑150 (MQ=255)
                                                                                        cTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGca                                                           >  1:367478/1‑150 (MQ=255)
                                                                                          ttCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGcacc                                                         <  2:1517871/150‑1 (MQ=255)
                                                                                              ccGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTa                                                     >  2:1051929/1‑150 (MQ=255)
                                                                                               cGTCTCATCGTTCATGGTTTTAGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCTTCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCATCGCACCACTaa                                                    <  1:766606/150‑1 (MQ=255)
                                                                                                   tcATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACt                                                                      >  2:96461/1‑128 (MQ=255)
                                                                                                   tcATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACt                                                                      <  1:102287/128‑1 (MQ=255)
                                                                                                   tcATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACt                                                                      <  1:96461/128‑1 (MQ=255)
                                                                                                   tcATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACt                                                                      >  2:102287/1‑128 (MQ=255)
                                                                                                   tcATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTcc                                                >  1:764058/1‑150 (MQ=255)
                                                                                                    cATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTccc                                               >  2:1362238/1‑150 (MQ=255)
                                                                                                    cATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTccc                                               >  2:1098812/1‑150 (MQ=255)
                                                                                                    cATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTccc                                               >  2:1606150/1‑150 (MQ=255)
                                                                                                               gTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTa                                    >  2:1310781/1‑150 (MQ=255)
                                                                                                               gTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTa                                    <  1:730823/150‑1 (MQ=255)
                                                                                                                  ttCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTg                                 <  2:918637/150‑1 (MQ=255)
                                                                                                                    cGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGAc                               >  1:643902/1‑150 (MQ=255)
                                                                                                                     gCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACg                              >  2:1330344/1‑150 (MQ=255)
                                                                                                                     gCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACg                              >  1:1232161/1‑150 (MQ=255)
                                                                                                                        ttAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTg                           <  2:55008/150‑1 (MQ=255)
                                                                                                                        ttAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTg                           >  2:36814/1‑150 (MQ=255)
                                                                                                                        ttAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTg                           <  2:1108309/150‑1 (MQ=255)
                                                                                                                         tAACAGCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATTGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTgg                          <  2:83587/150‑1 (MQ=255)
                                                                                                                         tAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGcac                                                          <  1:1256059/118‑1 (MQ=255)
                                                                                                                         tAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGcac                                                          >  2:1256059/1‑118 (MQ=255)
                                                                                                                         tAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTgg                          >  2:307506/1‑150 (MQ=255)
                                                                                                                          aaCACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATATTTACCCACTGTACTGTTCCAGCGCACCACTATCTCCCGGCGTTCGGTACTGAAGTTggg                         <  2:1032439/150‑1 (MQ=255)
                                                                                                                            cacCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCt                       >  2:858518/1‑150 (MQ=255)
                                                                                                                               cTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGc                    <  1:1051929/150‑1 (MQ=255)
                                                                                                                               cTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCATTCGGTACTGACGTTGGGCTGGc                    <  2:414430/150‑1 (MQ=255)
                                                                                                                                tttCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGcc                   >  1:444736/1‑150 (MQ=255)
                                                                                                                                      ggcatCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTc             >  1:930178/5‑150 (MQ=255)
                                                                                                                                       gcatCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTcc            >  1:1244793/4‑150 (MQ=255)
                                                                                                                                       gCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTcc            >  2:360109/1‑150 (MQ=255)
                                                                                                                                        catCGATGCCGGTGCTACGGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt           >  1:179487/3‑150 (MQ=255)
                                                                                                                                        cTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCt           <  1:455967/150‑1 (MQ=255)
                                                                                                                                           cGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCGACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGAt        <  1:1330440/150‑1 (MQ=255)
                                                                                                                                           cGATGCCCGTAGCACCGGCATCGATGCCAAAAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGAt        >  2:342350/1‑150 (MQ=255)
                                                                                                                                            gATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATg       >  1:132682/1‑150 (MQ=255)
                                                                                                                                            gATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATg       <  2:19346/150‑1 (MQ=255)
                                                                                                                                              tGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCt     <  1:1271406/150‑1 (MQ=255)
                                                                                                                                                 ccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTtat  <  2:1431050/150‑1 (MQ=255)
                                                                                                                                                 ccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTtat  >  1:1260234/1‑150 (MQ=255)
                                                                                                                                                 ccGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTtat  >  1:444259/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ACATCCAGATCGAGGCGGAAATTATAAATATGTTGGTGTGTAGTACCCACGATATTGTGATCGATAAGCGTGCCGTAGCGCGTGTCATCTTTCGCCGTCTCATCGTGCATGGTTTTCGCTTTAACACCTTTCACCGCTTCGATGCCCGTAGCACCGGCATCGATGCCAATAGTGCCGTTTTCATGGAAGATCCAGTCAAAAATGTAGTCATAGTTACCCACTGTACTGATCCAGCGCACCACTAACTCCCGGCGTTCGGTACTGACGTTGGGCTGGCCCATTTCCTGATGCTTAT  >  NC_000913/1449653‑1449947

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