Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 891,425 C→T A81V (GCC→GTC)  nfsA → nitroreductase A, NADPH‑dependent, FMN‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913891,4250CT98.8% 268.2 / NA 86A81V (GCC→GTC) nfsAnitroreductase A, NADPH‑dependent, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (40/45);  total (40/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

AACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTGGGCGGGGTATATATC  >  NC_000913/891277‑891570
                                                                                                                                                    |                                                                                                                                                 
aaCAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCg                                                                                                                                                 <  2:1096467/151‑1 (MQ=255)
aaCAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGAGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCg                                                                                                                                                 >  2:596258/1‑151 (MQ=255)
    gcgcCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACtt                                                                                                                                             <  1:1114530/151‑1 (MQ=255)
           tGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGc                                                                                                                                      >  2:280228/1‑151 (MQ=255)
           tGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGGGTTCTGTGTCGACTTTaaac                                                                                                                                        >  2:894363/1‑147 (MQ=255)
            gCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGc                                                                                                                                      <  2:678946/150‑1 (MQ=255)
             cGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGGGTTCTGTGTCGACTTTAACCGcc                                                                                                                                     >  1:450151/1‑150 (MQ=255)
              gACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGCAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCAt                                                                                                                                   >  2:855416/1‑151 (MQ=255)
              gACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAACATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCa                                                                                                                                    <  2:639333/150‑1 (MQ=255)
               aCGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGGGTTCTGTGTCGACTTTAACCGCCAtt                                                                                                                                  >  1:129774/1‑151 (MQ=255)
                   ccagttccagttTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTa                                                                                                                                <  1:443960/149‑1 (MQ=255)
                    cagttccagttTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACGGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAg                                                                                                                             <  1:695237/151‑1 (MQ=255)
                      gttccagttTTTTGCAGTGCAGTAGCATTATTCGCATTCCCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGAt                                                                                                                           <  2:660810/151‑1 (MQ=255)
                        tccagttTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCt                                                                                                                         >  1:227469/1‑151 (MQ=255)
                         ccagttTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATc                                                                                                                          >  1:687497/1‑149 (MQ=255)
                         ccagttTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTg                                                                                                                        <  2:243777/151‑1 (MQ=255)
                         ccagttTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTg                                                                                                                        <  1:119386/151‑1 (MQ=255)
                          cagttTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGt                                                                                                                       <  2:1182269/151‑1 (MQ=255)
                           agttTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTc                                                                                                                      >  1:862878/1‑151 (MQ=255)
                            gttTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTc                                                                                                                      <  2:600130/150‑1 (MQ=255)
                             ttttttGCAGGCGAGTGGCATTATCCCCATCGCAGACAATGGGTTACGGGAAGAAATGTGGACCGTGACCGGCGGGCAAAACCAGGTAGCGAAAGCGGCGGAGTTTTGGGTTTTCTGTGCGTACTTTAACCGCCATTTACAGAt                                                                                                                           <  1:539472/144‑1 (MQ=255)
                              tttttGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCgg                                                                                                                   >  1:513697/1‑151 (MQ=255)
                                 ttGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGAt                                                                                                                 <  2:61354/150‑1 (MQ=255)
                                 ttGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGAt                                                                                                                 <  2:617320/150‑1 (MQ=255)
                                    cagtgcagtAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTc                                                                                                             >  2:3357/1‑151 (MQ=255)
                                    cagtgcagtAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTc                                                                                                             >  2:1109642/1‑151 (MQ=255)
                                     agtgcagtAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTc                                                                                                             <  2:28068/150‑1 (MQ=255)
                                     agagaagtAGCTTTATTCGCATTACCGACAAAGGGTTGCGGAAAGAACTGGTGACGCGGACCGGCGGACAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTTTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCa                                                                                                            <  2:856824/146‑1 (MQ=255)
                                      gtgcagtAGCATTATTCGCATTACCGACAAAGCGTTCCGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCGGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCa                                                                                                            <  1:757541/150‑1 (MQ=255)
                                           gtAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCgg                                                                                                      <  1:327603/151‑1 (MQ=255)
                                               cATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCt                                                                                                   <  2:539054/150‑1 (MQ=255)
                                               cATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTg                                                                                                  >  1:475287/1‑151 (MQ=255)
                                                attattCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTgg                                                                                                 >  2:1128952/1‑151 (MQ=255)
                                                attattCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTg                                                                                                  <  1:1138021/150‑1 (MQ=255)
                                                 ttattCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGc                                                                                                >  2:508738/1‑151 (MQ=255)
                                                   attCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCgg                                                                                              >  2:1125296/1‑151 (MQ=255)
                                                      cGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGaac                                                                                           <  2:879200/151‑1 (MQ=255)
                                                      cGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGa                                                                                             <  1:135121/149‑1 (MQ=255)
                                                        cATTACCGACAAAGCGATACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAACGCAAGCGGCGGAGTTCTGGGAGTTCTGTGTTGACTTTAACCGCCAATTACAGAACTGTCCGGATCATCAGCTCGGCCtctccgaac                                                                                           >  1:820650/1‑141 (MQ=255)
                                                          ttACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGCCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGaacaa                                                                                         <  1:767895/149‑1 (MQ=255)
                                                          ttACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACt                                                                                       >  1:478995/1‑151 (MQ=255)
                                                           tACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTg                                                                                                  <  1:645065/139‑1 (MQ=255)
                                                           tACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTg                                                                                                  >  2:645215/1‑139 (MQ=255)
                                                           tACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTg                                                                                      >  1:957854/1‑151 (MQ=255)
                                                           tACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTg                                                                                      >  1:667697/1‑151 (MQ=255)
                                                           tACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGAGGCGGAGATCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGACTGGCGGAACAACt                                                                                       >  1:642652/1‑150 (MQ=255)
                                                             ccGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGtt                                                                                    >  1:1117432/1‑151 (MQ=255)
                                                             ccGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGt                                                                                     <  1:1060844/150‑1 (MQ=255)
                                                              cGACAAAGCGTTACGTGAAGAACTGGTGCCGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTg                                                                                   <  1:399255/151‑1 (MQ=255)
                                                              cGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGaacaa                                                                                         >  1:922445/1‑145 (MQ=255)
                                                              cGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACAAGTAGCGCAAGCGGCGGAGCTCTGGGGGGTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGAAGCTCAGCCCGGCCTGGCGGGACAACTggt                                                                                    >  1:659745/1‑148 (MQ=255)
                                                               gACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGGTCAGCTCGGCCTGGCGGAACAACTGTTGc                                                                                  >  1:331855/1‑151 (MQ=255)
                                                                aCAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCt                                                                                 <  2:536067/151‑1 (MQ=255)
                                                                aCAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGGTCTGGGGGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGc                                                                                  >  1:159770/1‑150 (MQ=255)
                                                                    aGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGgt                                                                             <  1:514869/151‑1 (MQ=255)
                                                                     gCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGgtg                                                                            <  1:865721/151‑1 (MQ=255)
                                                                     gCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGgt                                                                             <  1:63162/150‑1 (MQ=255)
                                                                      cGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGgtgt                                                                           >  2:291330/1‑151 (MQ=255)
                                                                      cGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGgtg                                                                            >  2:29473/1‑150 (MQ=255)
                                                                            gTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTg                                                                      <  1:553993/150‑1 (MQ=255)
                                                                              gaagaaCTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGAt                                                                    <  2:129805/150‑1 (MQ=255)
                                                                                     tGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAat                                                            >  1:351484/1‑151 (MQ=255)
                                                                                      ggTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAat                                                            <  1:993050/150‑1 (MQ=255)
                                                                                          aCGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATgg                                                       <  1:814981/151‑1 (MQ=255)
                                                                                               gACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGgcgc                                                    >  2:378550/1‑149 (MQ=255)
                                                                                                       ggCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCAtt                                          <  1:410042/151‑1 (MQ=255)
                                                                                                         cAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTa                                         >  2:588945/1‑150 (MQ=255)
                                                                                                           aaaaCACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATc                                      <  1:978229/151‑1 (MQ=255)
                                                                                                             aaCACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGc                                    <  2:475403/151‑1 (MQ=255)
                                                                                                               cacGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGc                                    >  1:268923/1‑149 (MQ=255)
                                                                                                                 cGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGCCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAg                                  <  1:385502/149‑1 (MQ=255)
                                                                                                                     gcgcAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAAt                            >  1:683732/1‑151 (MQ=255)
                                                                                                                      cgcAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAAt                            <  1:274741/150‑1 (MQ=255)
                                                                                                                        cAAGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGc                         >  2:548434/1‑151 (MQ=255)
                                                                                                                         aaGCGGCGGAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCt                        <  1:1125612/151‑1 (MQ=255)
                                                                                                                              gcggAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGAt                   >  2:113302/1‑151 (MQ=255)
                                                                                                                                ggAGTTCTGGGTGTGCTGTGTCGACTTAAACCGCCATTTACAGATCTGTGCGGATGCTCAGCTCGGCCTTGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGAtt                  <  2:368658/150‑1 (MQ=255)
                                                                                                                                 gAGTTCTGGGTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTgg                >  2:309365/1‑151 (MQ=255)
                                                                                                                                    ttCTGGGTGTTCTGTGTCGCCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTggg               <  1:643279/149‑1 (MQ=255)
                                                                                                                                     tCTGGGAGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTGGGCgg            >  2:1123637/1‑151 (MQ=255)
                                                                                                                                        gggTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTGGGCGGGGt         <  2:919569/151‑1 (MQ=255)
                                                                                                                                        gggTGTTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTGGGCGGGGt         <  1:433898/151‑1 (MQ=255)
                                                                                                                                          gtgtTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTGGGCGGGGtat       <  2:180500/151‑1 (MQ=255)
                                                                                                                                          gtgtTCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTGGGCGGGGtat       <  2:277487/151‑1 (MQ=255)
                                                                                                                                             ttCTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTGGGCGGGGtata      >  2:803026/1‑149 (MQ=255)
                                                                                                                                               cTGTGTCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTGGGCGGGGTATATATc  >  2:661313/1‑151 (MQ=255)
                                                                                                                                                 gtgtCGCCTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTGGGCGGGGTATATATc  <  2:682197/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AACAGCGCCCGTGCGACGTCCAGTTCCAGTTTTTTGCAGTGCAGTAGCATTATTCGCATTACCGACAAAGCGTTACGTGAAGAACTGGTGACGCTGACCGGCGGGCAAAAACACGTAGCGCAAGCGGCGGAGTTCTGGGTGTTCTGTGCCGACTTTAACCGCCATTTACAGATCTGTCCGGATGCTCAGCTCGGCCTGGCGGAACAACTGTTGCTCGGTGTCGTTGATACGGCAATGATGGCGCAGAATGCATTAATCGCAGCGGAATCGCTGGGATTGGGCGGGGTATATATC  >  NC_000913/891277‑891570

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