Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,363,075 A→G K111K (AAA→AAG puuC → gamma‑glutamyl‑gamma‑aminobutyraldehyde dehydrogenase; succinate semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,363,0750AG95.3% 296.1 / 2.9 86K111K (AAA→AAGpuuCgamma‑glutamyl‑gamma‑aminobutyraldehyde dehydrogenase; succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (42/40);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-01

GATGAGCGCAGCACGCGGCGTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAACCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTG  >  NC_000913/1362928‑1363222
                                                                                                                                                   |                                                                                                                                                   
gATGAGCGCAGCACGCGGCGTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCg                                                                                                                                                  >  1:280117/1‑151 (MQ=255)
    aGCGCAGCACGCGGCGTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTc                                                                                                                                              <  1:880716/151‑1 (MQ=255)
      cgcAGCACGCGGCGTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGt                                                                                                                                            >  2:1076001/1‑151 (MQ=255)
        cagcaCGCGGCGTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTc                                                                                                                                           >  2:326360/1‑150 (MQ=255)
           caCGCGGCGTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTcaca                                                                                                                                        <  1:849863/150‑1 (MQ=255)
             cgcgGCGTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAACCGATTCGTCACAGTc                                                                                                                                     >  1:237907/1‑151 (MQ=255)
              gcggcgTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCt                                                                                                                                    >  2:433861/1‑151 (MQ=255)
              gcggcgTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCt                                                                                                                                    <  1:1029769/151‑1 (MQ=255)
               cggcgTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTg                                                                                                                                   >  2:238291/1‑151 (MQ=255)
                ggcgTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGc                                                                                                                                  <  2:806467/151‑1 (MQ=255)
                ggcgTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGc                                                                                                                                  >  1:1086363/1‑151 (MQ=255)
                ggcgTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGc                                                                                                                                  >  2:522902/1‑151 (MQ=255)
                        taaaCGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGtgatga                                                                                                                           <  2:333800/148‑1 (MQ=255)
                               ggcgACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAACCGATTCGTCACAGTCTGCGTGATGATATTCCcg                                                                                                                   >  2:1248217/1‑151 (MQ=255)
                                gcgACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCcgg                                                                                                                  <  1:336878/151‑1 (MQ=255)
                                 cgACTGGTCACTCGCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGCAAACTCTCGACACCGGCAAGCCGATTCGTCTCAGTCTGCGTGATGATATTCCcggc                                                                                                                 <  2:930311/151‑1 (MQ=255)
                                  gACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCcggcg                                                                                                                >  1:59968/1‑151 (MQ=255)
                                   aCTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCcggcgc                                                                                                               <  1:168789/151‑1 (MQ=255)
                                      ggTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCcggcgcgg                                                                                                             <  1:117598/150‑1 (MQ=255)
                                      ggTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCcggcgcg                                                                                                              <  1:812238/149‑1 (MQ=255)
                                         cACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGgcgcggcgcg                                                                                                         >  2:520855/1‑151 (MQ=255)
                                         cACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGgcgcggcgcg                                                                                                         <  2:225470/151‑1 (MQ=255)
                                         cACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGgcgcggcgcg                                                                                                         >  2:969091/1‑151 (MQ=255)
                                             ctctTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGgcgcgcg                                                                                                       <  1:1079683/149‑1 (MQ=255)
                                               cttctCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGgcgcgcgc                                                                                                      <  2:920204/148‑1 (MQ=255)
                                                 tctcCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTc                                                                                                 >  2:389633/1‑151 (MQ=255)
                                                    ccGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGc                                                                                               >  1:707275/1‑150 (MQ=255)
                                                     cGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTg                                                                                             <  1:1087020/151‑1 (MQ=255)
                                                      ggCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTgg                                                                                            <  2:995294/151‑1 (MQ=255)
                                                       gCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGt                                                                                           <  2:853897/151‑1 (MQ=255)
                                                          aaaCGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTAc                                                                                         <  2:1218636/150‑1 (MQ=255)
                                                          aaaCGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACg                                                                                        <  1:101063/151‑1 (MQ=255)
                                                          aaaCGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACg                                                                                        >  2:1020251/1‑151 (MQ=255)
                                                          aaaCGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACg                                                                                        <  2:1150643/151‑1 (MQ=255)
                                                           aaCGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCGCTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGc                                                                                       <  2:1058298/151‑1 (MQ=255)
                                                            aCGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGcc                                                                                      <  2:332820/151‑1 (MQ=255)
                                                            aCGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGcc                                                                                      >  1:913964/1‑151 (MQ=255)
                                                             cGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGcc                                                                                      <  1:692770/150‑1 (MQ=255)
                                                                aaaGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGa                                                                                    >  1:459685/1‑149 (MQ=255)
                                                                      gTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCg                                                                            <  1:855046/151‑1 (MQ=255)
                                                                       tACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCg                                                                            >  1:206441/1‑150 (MQ=255)
                                                                           gAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACaa                                                                        <  1:101317/150‑1 (MQ=255)
                                                                           gAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACaa                                                                        >  2:712483/1‑150 (MQ=255)
                                                                             aTAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAgt                                                                     >  1:837860/1‑151 (MQ=255)
                                                                             aTAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAACCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAg                                                                      >  2:436638/1‑150 (MQ=255)
                                                                               aaaCTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAgtgt                                                                   <  1:799571/151‑1 (MQ=255)
                                                                               aaaCTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAgtgt                                                                   >  1:451500/1‑151 (MQ=255)
                                                                               aaaCTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAgtgt                                                                   >  1:910713/1‑151 (MQ=255)
                                                                               aaaCTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAgtg                                                                    <  2:515988/150‑1 (MQ=255)
                                                                                aaCTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAgtgt                                                                   >  1:339179/1‑150 (MQ=255)
                                                                                 aCTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTAt                                                                 >  1:8071/1‑151 (MQ=255)
                                                                                 aCTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTAt                                                                 <  2:1045593/151‑1 (MQ=255)
                                                                                  cTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATg                                                                >  1:761981/1‑151 (MQ=255)
                                                                                     gCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCg                                                             >  1:472872/1‑151 (MQ=255)
                                                                                     gCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCg                                                             >  2:508431/1‑151 (MQ=255)
                                                                                           ttAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTgg                                                       >  1:366518/1‑151 (MQ=255)
                                                                                            tAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTgg                                                       <  1:452446/150‑1 (MQ=255)
                                                                                                 gAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGacca                                                 <  1:15794/151‑1 (MQ=255)
                                                                                                    gCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCTCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGGGGCGaccac                                                >  1:220808/1‑149 (MQ=255)
                                                                                                    gCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGaccacc                                               >  2:345292/1‑150 (MQ=255)
                                                                                                      ccACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGt                                            >  2:129001/1‑151 (MQ=255)
                                                                                                      ccACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGt                                            >  1:619047/1‑151 (MQ=255)
                                                                                                      ccACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGt                                            >  2:1023387/1‑151 (MQ=255)
                                                                                                         cGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGc                                         >  2:125409/1‑151 (MQ=255)
                                                                                                             gAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGc                                         >  1:1083648/1‑147 (MQ=255)
                                                                                                               agagCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATg                                     >  1:1193016/1‑149 (MQ=255)
                                                                                                               agagCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAg                                   <  1:78720/151‑1 (MQ=255)
                                                                                                                gagCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCTTTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGc                                  <  2:942566/151‑1 (MQ=255)
                                                                                                                  gCTGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTg                                <  1:745610/151‑1 (MQ=255)
                                                                                                                    tGGCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGc                              >  2:1064153/1‑151 (MQ=255)
                                                                                                                     ggCACTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGc                              <  2:1161353/150‑1 (MQ=255)
                                                                                                                        aCTGCTGGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCgatg                          <  2:976539/151‑1 (MQ=255)
                                                                                                                          tgctgGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCgatgat                        >  1:751298/1‑151 (MQ=255)
                                                                                                                           gctgGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCgatgat                        <  2:585872/150‑1 (MQ=255)
                                                                                                                            ctgGAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGCGGCGACCACCAGCAGCCATGAGCTGGCGATGATc                       >  1:761021/1‑150 (MQ=255)
                                                                                                                              ggAACCTCTCAACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgt                     <  1:214806/150‑1 (MQ=255)
                                                                                                                              ggAAACTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcg                      <  1:1165607/149‑1 (MQ=255)
                                                                                                                              ggAAACTCTCGACACCGGCAAACCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgtg                    <  2:283330/151‑1 (MQ=255)
                                                                                                                                aaaCTCTCGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATcgtgcg                  >  1:961262/1‑151 (MQ=255)
                                                                                                                                      tcGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAAcc            >  2:271037/1‑151 (MQ=255)
                                                                                                                                      tcGACACCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAAcc            >  2:291497/1‑151 (MQ=255)
                                                                                                                                           acCGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCg       >  1:14382/1‑151 (MQ=255)
                                                                                                                                            ccGGCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCg       <  1:548629/150‑1 (MQ=255)
                                                                                                                                              ggCAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCg    <  2:885753/151‑1 (MQ=255)
                                                                                                                                                cAAGCCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTg  <  2:853193/151‑1 (MQ=255)
                                                                                                                                                cAAGCCGATTCGACACAGTCTGCGTGATGATATTCCCGCCGCGGCGCGCCCCATTCGCTGGTACGCCGAAGCGATCGACAAACAGGATGGCGAAGTGGCGACTACCACTAGTCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGc     >  1:397622/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATGAGCGCAGCACGCGGCGTATTTGAACGCGGCGACTGGTCACTCTCTTCTCCGGCTAAACGTAAAGCGGTACTGAATAAACTCGCCGATTTAATGGAAGCCCACGCCGAAGAGCTGGCACTGCTGGAAACTCTCGACACCGGCAAACCGATTCGTCACAGTCTGCGTGATGATATTCCCGGCGCGGCGCGCGCCATTCGCTGGTACGCCGAAGCGATCGACAAAGTGTATGGCGAAGTGGCGACCACCAGTAGCCATGAGCTGGCGATGATCGTGCGTGAACCGGTCGGCGTG  >  NC_000913/1362928‑1363222

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