Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,002,271 T→C 7.7% F282F (TTT→TTC corA → Ni(2(+))/Co(2(+))/Mg(2(+)) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,002,2710TC 7.7% 537.7 / 25.8 181F282F (TTT→TTCcorANi(2(+))/Co(2(+))/Mg(2(+)) transporter
Reads supporting (aligned to +/- strand):  ref base T (93/74);  new base C (6/8);  total (99/82)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 185 total aligned reads displayed.

CTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTC  >  NC_000913/4002122‑4002420
                                                                                                                                                     |                                                                                                                                                     
     ccAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGcc                                                                                                                                                  >  1:1287427/1‑150 (MQ=255)
     ccAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCAAAGTTCATGcc                                                                                                                                                  <  1:1055052/150‑1 (MQ=255)
      cAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCa                                                                                                                                                 >  2:1110646/1‑150 (MQ=255)
      cAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCa                                                                                                                                                 <  1:838974/150‑1 (MQ=255)
      cAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCa                                                                                                                                                 >  2:898711/1‑150 (MQ=255)
       aGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAg                                                                                                                                                >  1:151773/1‑150 (MQ=255)
        taaGGTGAAATTCCTGATGCATGCTGCTATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCAAAGTTCATGCCata                                                                                                                                               <  1:1120444/149‑3 (MQ=255)
         aaGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCAAAGTTCATGCCatag                                                                                                                                              <  1:1094747/150‑4 (MQ=255)
         aaGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCAAAGTTCATGCCatag                                                                                                                                              <  2:1380346/150‑4 (MQ=255)
             tGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGa                                                                                                                                          <  1:161090/150‑1 (MQ=255)
                aCTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATAAACTCAAAGTTCATGCCATAGCTggaag                                                                                                                                       <  2:526048/150‑5 (MQ=255)
                aCTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTACTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGt                                                                                                                                       >  1:412312/1‑150 (MQ=255)
                 cTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTg                                                                                                                                      >  1:293394/1‑150 (MQ=255)
                  ttCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTgg                                                                                                                                     <  1:1007349/150‑1 (MQ=255)
                   tCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTAt                                                                                                                                                     <  1:1186144/133‑1 (MQ=255)
                   tCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGa                                                                                                                                    >  1:577523/1‑150 (MQ=255)
                    ccTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAg                                                                                                                                   <  1:1427147/150‑1 (MQ=255)
                     cTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGc                                                                                                                                  >  2:1201422/1‑150 (MQ=255)
                      tGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCt                                                                                                                                 >  1:186729/1‑150 (MQ=255)
                      tGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCt                                                                                                                                 <  2:499998/150‑1 (MQ=255)
                          gCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCgg                                                                                                                             >  2:36561/1‑150 (MQ=255)
                            aGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCt                                                                                                                           >  1:15353/1‑150 (MQ=255)
                             ggcggcAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTa                                                                                                                          >  1:754520/1‑150 (MQ=255)
                             ggcggcAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATATTCTCGGTGGTATCCGTGGTATTCATTCCGCCGACGCTCGTTGCTTACAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGTAGCTTCGGCTa                                                                                                                          >  2:191264/1‑150 (MQ=255)
                                            tATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTAt                                                                                                           >  2:413521/1‑150 (MQ=255)
                                             aTCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTgg                                                                                                                    <  2:12546/140‑1 (MQ=255)
                                             aTCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTgg                                                                                                                    >  1:12546/1‑140 (MQ=255)
                                              tCAACATCGAGCAGAACCGCATCATCACAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCt                                                                                                         <  2:386531/150‑1 (MQ=255)
                                               cAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGGGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCtt                                                                                                        <  1:39628/150‑1 (MQ=255)
                                               cAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCtt                                                                                                        >  2:133040/1‑150 (MQ=255)
                                               cAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCtt                                                                                                        <  1:944600/150‑1 (MQ=255)
                                                  cATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTAt                                                                                                     >  1:1038813/1‑150 (MQ=255)
                                                  cATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTAt                                                                                                     >  2:878464/1‑150 (MQ=255)
                                                   aTCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTTGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGATTTTATGCCAGAACTGAAGTTGAGCTTCTGCTACGCTGGCGCGATTATCTTTATg                                                                                                    <  1:740410/150‑1 (MQ=255)
                                                    tCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGa                                                                                                   <  2:389728/150‑1 (MQ=255)
                                                        gCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCt                                                                                               <  2:362728/150‑1 (MQ=255)
                                                          aGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTcg                                                                                             >  2:1427923/1‑150 (MQ=255)
                                                           gAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCAAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTATGATCCTcgcg                                                                                           >  2:723253/1‑150 (MQ=255)
                                                              ccGCATCATCAAAATCTTGGCGGTGGTATCCGGGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGATCTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCggg                                                                                         <  1:407744/150‑1 (MQ=255)
                                                                gCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGcc                                                                                       >  1:723359/1‑150 (MQ=255)
                                                                gCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGcc                                                                                       >  2:772900/1‑150 (MQ=255)
                                                                 catcatCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCt                                                                                      <  2:630921/150‑1 (MQ=255)
                                                                  atcatcAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTg                                                                                     <  2:604501/150‑1 (MQ=255)
                                                                    catcaAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGc                                                                                   >  1:1307666/1‑150 (MQ=255)
                                                                      tcaAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCAc                                                                                 >  1:113948/1‑150 (MQ=255)
                                                                        aaaaTCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTa                                                                                                            >  1:611406/1‑121 (MQ=255)
                                                                         aaaTCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGt                                                                              <  2:113948/150‑1 (MQ=255)
                                                                            tcttctCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATc                                                                           <  1:90467/150‑1 (MQ=255)
                                                                              ttctCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTg                                                                         <  2:950293/150‑1 (MQ=255)
                                                                               tctcGGTGGTATCCGTGGTATTACTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGt                                                                        >  2:561820/1‑150 (MQ=255)
                                                                                   ggtggtATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACtt                                                                    >  2:115765/1‑150 (MQ=255)
                                                                                     tggtATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTa                                                                  <  1:491167/150‑1 (MQ=255)
                                                                                      ggtATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTaa                                                                 >  1:1345487/1‑150 (MQ=255)
                                                                                      ggtATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTaa                                                                 <  1:1164748/150‑1 (MQ=255)
                                                                                      ggtATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTaa                                                                 >  1:645231/1‑150 (MQ=255)
                                                                                         aTCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCg                                                              >  2:1112987/1‑150 (MQ=255)
                                                                                         aTCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCg                                                              <  1:36561/150‑1 (MQ=255)
                                                                                         aTCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCg                                                              >  2:184664/1‑150 (MQ=255)
                                                                                          tCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCgg                                                             >  2:741776/1‑150 (MQ=255)
                                                                                           ccGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGga                                                            >  1:1231717/1‑150 (MQ=255)
                                                                                           ccGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGga                                                            >  2:1265900/1‑150 (MQ=255)
                                                                                            cGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaa                                                           >  1:499193/1‑150 (MQ=255)
                                                                                              tGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaaga                                                         <  2:754520/150‑1 (MQ=255)
                                                                                               ggTATTCCTGCCGCCGACGCTCGTTGATTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaagaa                                                        >  2:433458/1‑150 (MQ=255)
                                                                                                 tATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACt                                                      <  2:499193/150‑1 (MQ=255)
                                                                                                  aTTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTg                                                     >  1:1333099/1‑150 (MQ=255)
                                                                                                   ttCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTgg                                                    <  2:916846/150‑1 (MQ=255)
                                                                                                          cgccgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaaa                                             >  1:338770/1‑150 (MQ=255)
                                                                                                           gccgACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTaaaa                                            >  1:821582/1‑150 (MQ=255)
                                                                                                              gACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAAc                                         >  2:216349/1‑150 (MQ=255)
                                                                                                              gACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAAc                                         <  2:577523/150‑1 (MQ=255)
                                                                                                                cGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGCGCTTCGGCTATCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACga                                       >  1:96673/1‑150 (MQ=255)
                                                                                                                  cTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACgaga                                     >  2:153850/1‑150 (MQ=255)
                                                                                                                   tCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGagag                                    >  2:1319840/1‑150 (MQ=255)
                                                                                                                     gTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCg                                  >  2:666426/1‑150 (MQ=255)
                                                                                                                       tGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaag                                                          <  2:240897/124‑1 (MQ=255)
                                                                                                                       tGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGgaag                                                          >  1:240897/1‑124 (MQ=255)
                                                                                                                        gCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTg                               <  1:133040/150‑1 (MQ=255)
                                                                                                                           tCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCt                            >  2:318012/1‑150 (MQ=255)
                                                                                                                            ccAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCtt                           <  1:665249/150‑1 (MQ=255)
                                                                                                                              aGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAg                         >  2:1043487/1‑150 (MQ=255)
                                                                                                                               gttctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGt                        >  1:924024/5‑150 (MQ=255)
                                                                                                                               gttctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGt                        >  1:1365431/5‑150 (MQ=255)
                                                                                                                                ttctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTc                       >  1:20315/4‑150 (MQ=255)
                                                                                                                                 tctGGCATAAACTCAAAGTTCATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCt                      >  1:502383/3‑150 (MQ=255)
                                                                                                                                 tATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCt                      >  2:1007798/1‑150 (MQ=255)
                                                                                                                                      cATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTa                 >  2:697313/1‑150 (MQ=255)
                                                                                                                                          aaCTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGcc             <  2:64935/150‑1 (MQ=255)
                                                                                                                                            cTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCAc           >  2:1077921/1‑150 (MQ=255)
                                                                                                                                             tttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCAcc          >  1:1400406/1‑150 (MQ=255)
                                                                                                                                             tttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCAcc          >  1:920078/1‑150 (MQ=255)
                                                                                                                                             tttGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAcc          <  2:250658/150‑1 (MQ=255)
                                                                                                                                               tGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAccgc        <  2:338770/150‑3 (MQ=255)
                                                                                                                                               tGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGCCAGACTAAGCCAccgc        <  2:293394/150‑3 (MQ=255)
                                                                                                                                                 aGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGtt      <  2:450184/150‑1 (MQ=255)
                                                                                                                                                  gTTTATGCCAGAACTGAAGTGGAGCTTCGGTTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTa     <  1:31356/150‑1 (MQ=255)
                                                                                                                                                  gTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTa     <  2:668074/150‑1 (MQ=255)
                                                                                                                                                  gTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTa     >  1:1144889/1‑150 (MQ=255)
                                                                                                                                                   tttATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTAt    <  1:433458/150‑1 (MQ=255)
                                                                                                                                                     tATGCCAGAACTGAAGGGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGTAAGAACTGGTTGTAAAAAACGACAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTc  <  2:895004/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGTTCCAGAAGGTGAACTTCCTGATGCAGGCGGCAATGGGTTTTATCAACATCGAGCAGAACCGCATCATCAAAATCTTCTCGGTGGTATCCGTGGTATTCCTGCCGCCGACGCTCGTTGCTTCCAGCTATGGCATGAACTTTGAGTTTATGCCAGAACTGAAGTGGAGCTTCGGCTACCCTGGCGCGATTATCTTTATGATCCTCGCGGGCCTGGCACCGTATCTGTACTTTAAGCGGAAGAACTGGTTGTAAAAAACGAGAGCGGTGGCTTAGTCTGGCTAAGCCACCTGTTATTC  >  NC_000913/4002122‑4002420

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: