Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,143,679 C→A 12.5% S162S (TCG→TCT hybB ← putative hydrogenase 2 cytochrome b type component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,143,6790CA12.5% 410.0 / 38.1 128S162S (TCG→TCThybBputative hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  ref base C (53/59);  new base A (7/9);  total (60/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 129 total aligned reads displayed.

AGCGGCAACATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTAGAAGTACGGCAGGTTCCAG  >  NC_000913/3143533‑3143826
                                                                                                                                                  |                                                                                                                                                   
tgCGTCAGCATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCTGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAgaga                                                                                                                                                   <  1:668747/148‑1 (MQ=255)
   ggCAACATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGacac                                                                                                                                                <  2:160478/149‑1 (MQ=255)
   ggCAACATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGAGacac                                                                                                                                                <  3:538162/149‑1 (MQ=255)
    gCAACATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTACAGCGactaa                                                                                                                                               <  4:364082/149‑4 (MQ=255)
      aaCATTTCATAGCTCTGCCACAACGGATTCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac                                                                                                                                                  >  3:61695/1‑144 (MQ=255)
      aaCATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGac                                                                                                                                                  <  4:61695/144‑1 (MQ=255)
       aCATTGCATGGCTCTGCCACAACGGATGCGCCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGAGACACCTTc                                                                                                                                            <  4:501150/149‑1 (MQ=255)
         aTTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCa                                                                                                                                          >  2:465385/1‑149 (MQ=255)
          tttCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCAc                                                                                                                                         <  3:511842/149‑1 (MQ=255)
              aTAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCCTCCACCCCa                                                                                                                                     <  3:424828/149‑1 (MQ=255)
                 gCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGAc                                                                                                                                  <  4:766516/149‑1 (MQ=255)
                  ctctGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACg                                                                                                                                 >  4:802312/1‑149 (MQ=255)
                   tctGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGt                                                                                                                                >  1:464379/1‑149 (MQ=255)
                    ctGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGtt                                                                                                                               <  3:434222/149‑1 (MQ=255)
                     tGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTc                                                                                                                              <  1:896216/149‑1 (MQ=255)
                       ccACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCaa                                                                                                                            >  4:501847/1‑149 (MQ=255)
                        cacaACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCaaa                                                                                                                           >  2:51652/1‑149 (MQ=255)
                         acaACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAAc                                                                                                                          >  4:594596/1‑149 (MQ=255)
                           aaCGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAg                                                                                                                        <  1:298398/149‑1 (MQ=255)
                           aaCGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAg                                                                                                                        >  3:186179/1‑149 (MQ=255)
                            aCGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGt                                                                                                                       <  3:905204/149‑1 (MQ=255)
                             cGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTg                                                                                                                      <  3:291543/149‑1 (MQ=255)
                                  gCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGa                                                                                                                 <  3:27941/149‑1 (MQ=255)
                                  gCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGa                                                                                                                 <  2:757022/149‑1 (MQ=255)
                                   cACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAg                                                                                                                >  2:39195/1‑149 (MQ=255)
                                    aCCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGc                                                                                                               >  1:139971/1‑149 (MQ=255)
                                    aCCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGc                                                                                                               >  1:937659/1‑149 (MQ=255)
                                    aCCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGc                                                                                                               >  3:123721/1‑149 (MQ=255)
                                     ccTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCa                                                                                                              <  1:658152/149‑1 (MQ=255)
                                          tAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTc                                                                                                         <  2:261503/149‑1 (MQ=255)
                                            gCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCa                                                                                                       <  2:139971/149‑1 (MQ=255)
                                            gCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCa                                                                                                       <  3:420305/149‑1 (MQ=255)
                                            gCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCa                                                                                                       <  4:294711/149‑1 (MQ=255)
                                               cgccgAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTg                                                                                                    >  1:717550/1‑149 (MQ=255)
                                               cgccgAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTg                                                                                                    >  4:175499/1‑149 (MQ=255)
                                                       tcatcaGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACg                                                                                            >  3:284294/1‑149 (MQ=255)
                                                            aGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGat                                                                                       <  1:789084/149‑1 (MQ=255)
                                                                 ccccATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGa                                                                                  >  4:54042/1‑149 (MQ=255)
                                                                  cccATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGAt                                                                                 >  4:401324/1‑149 (MQ=255)
                                                                  cccATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGAt                                                                                 <  4:62440/149‑1 (MQ=255)
                                                                  cccATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGAt                                                                                 >  4:726678/1‑149 (MQ=255)
                                                                  cccATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGAt                                                                                 >  1:856179/1‑149 (MQ=255)
                                                                   ccATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATg                                                                                >  2:194252/1‑149 (MQ=255)
                                                                   ccATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATg                                                                                >  3:846024/1‑149 (MQ=255)
                                                                    cATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATgg                                                                               >  1:471865/1‑149 (MQ=255)
                                                                    cATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATgg                                                                               <  1:302805/149‑1 (MQ=255)
                                                                      ttGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTc                                                                             >  2:721883/1‑149 (MQ=255)
                                                                         aaGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCa                                                                                                       >  2:755791/1‑120 (MQ=255)
                                                                         aaGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCa                                                                                                       <  1:755791/120‑1 (MQ=255)
                                                                         aaGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCGTCACGCCGATATAGATGGTCATa                                                                          <  3:802312/149‑1 (MQ=255)
                                                                         aaGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATa                                                                          <  2:748212/149‑1 (MQ=255)
                                                                         aaGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATa                                                                          <  1:389898/149‑1 (MQ=255)
                                                                         aaGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATa                                                                          >  4:396468/1‑149 (MQ=255)
                                                                             cTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGa                                                                      <  2:412941/149‑1 (MQ=255)
                                                                              tGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGAt                                                                                 <  3:218215/137‑1 (MQ=255)
                                                                              tGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGAt                                                                                 >  4:218215/1‑137 (MQ=255)
                                                                              tGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAc                                                                     <  4:62669/149‑1 (MQ=255)
                                                                                  gCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAGCAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccgcc                                                                 >  3:875052/1‑149 (MQ=255)
                                                                                   cATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGAccgccg                                                                >  1:238290/1‑149 (MQ=255)
                                                                                      ggTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtct                                                             >  1:106706/1‑149 (MQ=255)
                                                                                       gTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGtctc                                                            <  3:892432/149‑1 (MQ=255)
                                                                                          ggCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGaa                                                         <  3:125303/149‑1 (MQ=255)
                                                                                           gcagcaGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAAc                                                        >  3:673036/1‑149 (MQ=255)
                                                                                            cagcagCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACa                                                       <  3:130796/149‑1 (MQ=255)
                                                                                               cagcGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTa                                                    <  2:4982/149‑1 (MQ=255)
                                                                                                 gcgcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGa                                                                                                                 <  1:683292/86‑1 (MQ=255)
                                                                                                 gcgcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGa                                                                                                                 >  2:683292/1‑86 (MQ=255)
                                                                                                   gcACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGa                                                <  1:111846/149‑1 (MQ=255)
                                                                                                     aCCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATAGTCATACAGACCGCCGTCTCGAACAGTACCGAGt                                              >  1:768370/1‑149 (MQ=255)
                                                                                                       cGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTc                                            <  2:464379/149‑1 (MQ=255)
                                                                                                           cgcgATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGt                                        <  3:501847/149‑1 (MQ=255)
                                                                                                              gatgatGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGa                                     >  3:284388/1‑149 (MQ=255)
                                                                                                               atgatgAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa                                    <  1:762096/149‑1 (MQ=255)
                                                                                                               atgatgAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa                                    <  2:51419/149‑1 (MQ=255)
                                                                                                               atgatgAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGaa                                    <  4:19946/149‑1 (MQ=255)
                                                                                                                  atgaAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAATGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGGCTCGAACAGTACCGAGTTCACGTTGAAGTg                                 >  2:292491/1‑149 (MQ=255)
                                                                                                                  atgaAGAACATCACCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTg                                 <  1:455600/149‑1 (MQ=255)
                                                                                                                    gaagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAg                                               >  2:413087/1‑133 (MQ=255)
                                                                                                                    gaagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAg                                               <  1:413087/133‑1 (MQ=255)
                                                                                                                    gaagaaCATCACCTTGTGTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCCAGTTCACGTTGAAGTGAc                               <  2:717550/149‑1 (MQ=255)
                                                                                                                     aagaaCATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCAAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGAcc                              >  1:433628/1‑149 (MQ=255)
                                                                                                                          cATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGaa                         <  1:500434/149‑1 (MQ=255)
                                                                                                                          cATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGaa                         <  4:668619/149‑1 (MQ=255)
                                                                                                                           aTCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGt                                                                                                     <  4:16040/72‑1 (MQ=255)
                                                                                                                           aTCACCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAAt                        <  1:661438/149‑1 (MQ=255)
                                                                                                                             cACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGt                      >  2:500984/1‑149 (MQ=255)
                                                                                                                             cACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGt                      <  1:500984/149‑1 (MQ=255)
                                                                                                                             cACCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGt                      >  2:621288/1‑149 (MQ=255)
                                                                                                                              aCCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTa                     <  3:594596/149‑1 (MQ=255)
                                                                                                                              aCCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCGTACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTa                     >  4:232458/1‑149 (MQ=255)
                                                                                                                              aCCTTGTTTAGTCGCTGTAGAGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTa                     <  2:305457/149‑1 (MQ=255)
                                                                                                                                    tttAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTAGAAGTa               >  2:240054/1‑149 (MQ=255)
                                                                                                                                     ttAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTAGAAGTAc              <  4:748801/149‑1 (MQ=255)
                                                                                                                                     ggtgTCGCTACAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTAGAAGTAc              >  3:666777/4‑149 (MQ=255)
                                                                                                                                       aGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTAGAAGTACgg            <  1:731294/149‑1 (MQ=255)
                                                                                                                                          cgctacaGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTAGAAGTACGGCAg         <  4:666777/143‑1 (MQ=255)
                                                                                                                                          cGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTAGAAGTACGGCAg         >  2:907564/1‑149 (MQ=255)
                                                                                                                                           gCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTAGAAGTACGGCAgg        <  3:223472/149‑1 (MQ=255)
                                                                                                                                           gCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTAGAAGTACGGCAgg        >  1:706374/1‑149 (MQ=255)
                                                                                                                                                aGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTAGAAGTACGGCAGGTTCCa   <  2:471865/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGCGGCAACATTTCATAGCTCTGCCACAACGGATGCACCTTGTAGCCCGCCGAGATCATCAGCGACCCCATTGAAGACTGGTGCATGGTCGGCAGCAGCGCACCGAGCGCGATGATGAAGAACATCACCTTGTTTAGTCGCTGTAGCGACACCTTCCACCCCAGACGTTCAAACAGTGCCGGAGCAAACTCCAGTGCCATCACGCCGATATAGATGGTCATACAGACCGCCGTCTCGAACAGTACCGAGTTCACGTTGAAGTGACCCGGAATGTAGAAGTACGGCAGGTTCCAG  >  NC_000913/3143533‑3143826

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