Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,007,466 T→C 7.7% intergenic (+14/‑44) ygeW → / → ygeX putative carbamoyltransferase/2,3‑diaminopropionate ammonia lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,007,4660TC 7.7% 282.8 / 10.5 91intergenic (+14/‑44)ygeW/ygeXputative carbamoyltransferase/2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (42/42);  new base C (3/4);  total (45/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

TACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTC  >  NC_000913/3007318‑3007614
                                                                                                                                                    |                                                                                                                                                    
tACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGAt                                                                                                                                                      >  3:368851/1‑149 (MQ=255)
 aCCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATc                                                                                                                                                     <  1:284511/149‑1 (MQ=255)
          ctacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCg                                                                                                                                            <  3:730302/149‑1 (MQ=255)
              ctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGccc                                                                                                                                        <  3:221588/149‑1 (MQ=255)
               tacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCct                                                                                                                                       >  4:452679/1‑149 (MQ=255)
                acAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGcccga                                                                                                                                      <  2:279690/149‑3 (MQ=255)
                  aaaGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTa                                                                                                                                    >  1:534016/1‑149 (MQ=255)
                       agCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTcc                                                                                                                               <  4:705997/149‑1 (MQ=255)
                           aGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCagag                                                                                                                           >  2:727575/1‑149 (MQ=255)
                           aGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCagag                                                                                                                           >  1:44766/1‑149 (MQ=255)
                                gAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGCGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaa                                                                                                                      >  4:464283/1‑149 (MQ=255)
                                  aGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaaaa                                                                                                                    <  3:452679/149‑1 (MQ=255)
                                      ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggat                                                                                                                >  1:763668/1‑149 (MQ=255)
                                      ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggat                                                                                                                <  1:670988/149‑1 (MQ=255)
                                       tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTTGGCCCTCTATTTCCAGAGGCCAAAaggata                                                                                                               <  1:749190/149‑1 (MQ=255)
                                       tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggata                                                                                                               >  2:600799/1‑149 (MQ=255)
                                        atatCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggatag                                                                                                              <  2:79015/149‑1 (MQ=255)
                                             atcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGatat                                                                                                         >  4:105605/1‑149 (MQ=255)
                                              tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATg                                                                                                        >  4:250482/1‑149 (MQ=255)
                                                gccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAagga                                                                                                                 <  4:487/138‑1 (MQ=255)
                                                gccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAagga                                                                                                                 >  3:487/1‑138 (MQ=255)
                                                     gatgatCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGttt                                                                                                 >  1:80850/1‑149 (MQ=255)
                                                      atgatCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGtttt                                                                                                <  1:445923/149‑1 (MQ=255)
                                                       tgatCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTtc                                                                                               >  1:173778/1‑149 (MQ=255)
                                                        gatCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTtct                                                                                              >  1:523683/1‑149 (MQ=255)
                                                         atCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTtctc                                                                                             >  4:266042/1‑149 (MQ=255)
                                                               tCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGaa                                                                                       >  2:739801/1‑149 (MQ=255)
                                                                cccGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAg                                                                                      <  2:44766/149‑1 (MQ=255)
                                                                 ccGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGa                                                                                     >  1:29921/1‑149 (MQ=255)
                                                                 ccGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGa                                                                                     >  1:298945/1‑149 (MQ=255)
                                                                  cGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGAt                                                                                    >  4:231466/1‑149 (MQ=255)
                                                                   gTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAGAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCTGCACTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGAtt                                                                                   <  4:368851/149‑1 (MQ=255)
                                                                    tAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAATGGATAGGATATGTCCGTTTTCTCATTGAAGATTg                                                                                  <  2:68514/149‑1 (MQ=255)
                                                                        tACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGatat                                                                              <  2:523683/149‑1 (MQ=255)
                                                                           gCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGc                                                                           >  3:317355/1‑149 (MQ=255)
                                                                            ccAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGc                                                                                                                                          >  4:87534/1‑85 (MQ=255)
                                                                            ccAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGc                                                                                                                                          <  3:87534/85‑1 (MQ=255)
                                                                             cAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCg                                                                         <  3:250482/149‑1 (MQ=255)
                                                                             cAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCg                                                                         <  3:103577/149‑1 (MQ=255)
                                                                              aaaCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGa                                                                        <  3:748830/149‑1 (MQ=255)
                                                                              aaaCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGa                                                                        >  4:339822/1‑149 (MQ=255)
                                                                              aaaCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGa                                                                        <  2:140824/149‑1 (MQ=255)
                                                                              aaaCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGa                                                                        <  1:764430/149‑1 (MQ=255)
                                                                                aCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATa                                                                      <  1:697519/149‑1 (MQ=255)
                                                                                   aGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATc                                                                             >  4:13295/1‑139 (MQ=255)
                                                                                   aGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATc                                                                             <  3:13295/139‑1 (MQ=255)
                                                                                        cACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAAttt                                                              <  2:482426/149‑1 (MQ=255)
                                                                                         aCTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAAtttt                                                             >  1:189419/1‑149 (MQ=255)
                                                                                          cTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAAttttt                                                            <  2:695283/149‑1 (MQ=255)
                                                                                           tGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGTGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTCCTCATTGAAGATTGTTATCGGCGATATCaatcttcaa                                                           <  1:634207/149‑8 (MQ=255)
                                                                                             cTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCaa                                                         <  3:104190/149‑1 (MQ=255)
                                                                                             cTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCaa                                                         <  4:141820/149‑1 (MQ=255)
                                                                                                gAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACgg                                                      >  4:698207/1‑149 (MQ=255)
                                                                                                 aGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGc                                                     <  3:231466/149‑1 (MQ=255)
                                                                                                  gCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCg                                                    <  1:581408/149‑1 (MQ=255)
                                                                                                         ctgAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATc                                             <  2:534016/149‑1 (MQ=255)
                                                                                                           gAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCAc                                           <  3:698207/149‑1 (MQ=255)
                                                                                                            aaaGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGTGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCAcc                                          >  2:651744/1‑149 (MQ=255)
                                                                                                            aaaGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCAcc                                          >  1:723889/1‑149 (MQ=255)
                                                                                                            aaaGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCAcc                                          >  2:102709/1‑149 (MQ=255)
                                                                                                              agaagCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGc                                        <  3:28942/149‑1 (MQ=255)
                                                                                                              agaagCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGc                                        >  3:715412/1‑149 (MQ=255)
                                                                                                               gaagCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGct                                       >  1:208737/1‑149 (MQ=255)
                                                                                                                     cAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTc                                 <  4:445575/149‑1 (MQ=255)
                                                                                                                     cAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTc                                 >  3:693021/1‑149 (MQ=255)
                                                                                                                     cAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTc                                 <  2:763668/149‑1 (MQ=255)
                                                                                                                        gAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCa                              >  2:737016/1‑149 (MQ=255)
                                                                                                                        gAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCa                              <  3:339822/149‑1 (MQ=255)
                                                                                                                          aCGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCaaa                            <  2:208737/149‑1 (MQ=255)
                                                                                                                               tGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAg                       <  1:600799/149‑1 (MQ=255)
                                                                                                                                 aaaTAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATTTCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGc                     <  1:92775/149‑1 (MQ=255)
                                                                                                                                 aaaTAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCGCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGc                     <  2:90592/149‑1 (MQ=255)
                                                                                                                                     aaCCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCaaa                 >  1:33465/1‑149 (MQ=255)
                                                                                                                                      aCCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCa                                            <  1:140477/121‑1 (MQ=255)
                                                                                                                                      aCCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCa                                            >  2:140477/1‑121 (MQ=255)
                                                                                                                                       ccccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACt               >  4:372237/1‑149 (MQ=255)
                                                                                                                                        cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGAt                                                                                    <  4:371559/79‑1 (MQ=255)
                                                                                                                                        cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGAt                                                                                    >  3:371559/1‑79 (MQ=255)
                                                                                                                                        cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTg              >  2:711653/1‑149 (MQ=255)
                                                                                                                                        cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTg              >  2:731666/1‑149 (MQ=255)
                                                                                                                                        cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTg              >  3:13669/1‑149 (MQ=255)
                                                                                                                                        cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTg              >  3:545863/1‑149 (MQ=255)
                                                                                                                                        cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTg              <  3:601202/149‑1 (MQ=255)
                                                                                                                                          cTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgc            >  3:173785/1‑149 (MQ=255)
                                                                                                                                            cGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgcgc          <  1:739801/149‑1 (MQ=255)
                                                                                                                                             gggCCGACCAGCTGATCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgcgcg         >  1:766543/1‑149 (MQ=255)
                                                                                                                                              ggCCGACCAGCTGATCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTc                                 <  3:204039/124‑1 (MQ=255)
                                                                                                                                              ggCCGACCAGCTGATCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTc                                 >  4:204039/1‑124 (MQ=255)
                                                                                                                                              ggCCGACCAGCTGATCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgcgcgc        >  4:660320/1‑149 (MQ=255)
                                                                                                                                                   aTCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGtt   >  1:221032/1‑149 (MQ=255)
                                                                                                                                                   aTCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGtt   >  4:742466/1‑149 (MQ=255)
                                                                                                                                                   aTCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGtt   <  1:711653/149‑1 (MQ=255)
                                                                                                                                                    tCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAACTTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGttt  <  2:173778/149‑2 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTC  >  NC_000913/3007318‑3007614

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: