Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,143,128 A→G 5.1% F346S (TTC→TCC)  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,1280AG 5.1% 260.9 / 3.8 78F346S (TTC→TCC) hybBputative hydrogenase 2 cytochrome b type component
Reads supporting (aligned to +/- strand):  ref base A (51/23);  new base G (2/2);  total (53/25)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01

CTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGAACGACCAGCGGGAA  >  NC_000913/3142985‑3143276
                                                                                                                                               |                                                                                                                                                    
cTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAg                                                                                                                                                 <  3:264182/149‑1 (MQ=255)
 tCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAgg                                                                                                                                                >  4:417524/1‑149 (MQ=255)
 tCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAg                                                                                                                                                 >  1:413532/1‑148 (MQ=255)
 tCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAg                                                                                                                                                 <  2:413532/148‑1 (MQ=255)
   aTGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCg                                                                                                                                              >  2:207341/1‑149 (MQ=255)
   aTGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCg                                                                                                                                              >  1:32354/1‑149 (MQ=255)
    tGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGt                                                                                                                                             >  4:537557/1‑149 (MQ=255)
       tttGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAc                                                                                                                                          >  4:485905/1‑149 (MQ=255)
          gCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACACAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAAccgc                                                                                                                                       <  3:307361/149‑1 (MQ=255)
                cTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCg                                                                                                                                 >  1:94511/1‑149 (MQ=255)
                 tCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGGAGTAGGCGTAACCGCCGCCCgg                                                                                                                                >  4:336704/1‑149 (MQ=255)
                 tCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCgg                                                                                                                                >  2:216779/1‑149 (MQ=255)
                   aTGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGt                                                                                                                              >  2:556593/1‑149 (MQ=255)
                    tGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGtt                                                                                                                             >  2:383650/1‑149 (MQ=255)
                      aCGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGa                                                                                                                           >  4:456767/1‑149 (MQ=255)
                       cGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGaa                                                                                                                          <  1:216779/149‑1 (MQ=255)
                        gATTATGATCGTTTTGTTTTAAAGGAGGACGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAAt                                                                                                                         >  3:395916/1‑149 (MQ=255)
                        gATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAAt                                                                                                                         <  3:456767/149‑1 (MQ=255)
                        gATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAAt                                                                                                                         <  1:321026/149‑1 (MQ=255)
                              gATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGaa                                                                                                                                                     <  3:350547/115‑1 (MQ=255)
                              gATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGaa                                                                                                                                                     >  4:350547/1‑115 (MQ=255)
                              gATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGccacc                                                                                                                   >  4:142705/1‑149 (MQ=255)
                               aTCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGaa                                                                                                                                                     <  2:139871/114‑1 (MQ=255)
                               aTCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGaa                                                                                                                                                     >  1:139871/1‑114 (MQ=255)
                                 cGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGc                                                                                                                <  2:578317/149‑1 (MQ=255)
                                   ttttgttttAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGa                                                                                                              >  2:109233/1‑149 (MQ=255)
                                   ttttgttttAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGa                                                                                                              >  4:289507/1‑149 (MQ=255)
                                       gttttAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAg                                                                                                          <  3:369288/149‑1 (MQ=255)
                                       gttttAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAg                                                                                                          <  4:572248/149‑1 (MQ=255)
                                       gttttAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAg                                                                                                          <  2:382353/149‑1 (MQ=255)
                                        ttttAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAgg                                                                                                         >  3:8566/1‑149 (MQ=255)
                                         tttAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGt                                                                                                        <  3:550093/149‑1 (MQ=255)
                                           tAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCa                                                                                                      <  3:417524/149‑1 (MQ=255)
                                                gaggaAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGc                                                                                                 >  2:285475/1‑149 (MQ=255)
                                                  ggaAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCa                                                                                               >  4:131410/1‑149 (MQ=255)
                                                    aaGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGGAGTAGGCGTAACCGCCGCCCGGGTTGAATg                                                                                                                        <  1:580520/122‑1 (MQ=255)
                                                    aaGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGATCTTCCCAGGTCGGGGAGTAGGCGTAACCGCCGCCCGGGTTGAATg                                                                                                                        >  2:580520/1‑122 (MQ=255)
                                                     aGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGt                                                                                            <  3:499402/149‑1 (MQ=255)
                                                       tATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGGAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTg                                                                                          <  3:336704/149‑1 (MQ=255)
                                                            gCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAAt                                                                                                                         >  2:307157/1‑113 (MQ=255)
                                                            gCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAAt                                                                                                                         <  1:307157/113‑1 (MQ=255)
                                                            gCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAAt                                                                                                                         >  1:398849/1‑113 (MQ=255)
                                                            gCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAAt                                                                                                                         <  2:398849/113‑1 (MQ=255)
                                                             cAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAAc                                                                                    <  1:383650/149‑1 (MQ=255)
                                                             cAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAAc                                                                                    <  2:8523/149‑1 (MQ=255)
                                                              aGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAAcc                                                                                   >  3:14789/1‑149 (MQ=255)
                                                                  gACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAAc                                                                               >  2:459097/1‑149 (MQ=255)
                                                                   aCGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACa                                                                              >  3:121191/1‑149 (MQ=255)
                                                                     gAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGt                                                                            >  3:206943/1‑149 (MQ=255)
                                                                            tACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCa                                                                     >  3:33255/1‑149 (MQ=255)
                                                                                   tAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGa                                                              <  3:525654/149‑1 (MQ=255)
                                                                                    aaGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGAc                                                             >  2:493574/1‑149 (MQ=255)
                                                                                        gcAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACa                                                            <  1:325296/146‑1 (MQ=255)
                                                                                        gcAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACa                                                            >  2:325296/1‑146 (MQ=255)
                                                                                         cAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGaa                                                        <  2:94511/149‑1 (MQ=255)
                                                                                               tcAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCAt                                                  >  3:260678/1‑149 (MQ=255)
                                                                                                cAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATg                                                 >  2:601572/1‑149 (MQ=255)
                                                                                                 aaTAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGc                                                >  4:345111/1‑149 (MQ=255)
                                                                                                   tAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCgg                                              >  3:102286/1‑149 (MQ=255)
                                                                                                      ccACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGTAACAGCATGCGGGaa                                           >  3:204299/1‑149 (MQ=255)
                                                                                                        acaAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCGCGGGTTGAATGCCACCAGCGAATAGGTCTGACGCCAGGTTGCACAACCTCACAGTGCACTCAGTGCTGACAGGAACAGCATGCGGGAATc                                         >  3:386128/1‑149 (MQ=255)
                                                                                                                aaTAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCa                                 >  4:107032/1‑149 (MQ=255)
                                                                                                                 aTAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTAc                                    <  2:282911/145‑1 (MQ=255)
                                                                                                                 aTAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTAc                                    >  1:282911/1‑145 (MQ=255)
                                                                                                                  tAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGc                               >  4:413540/1‑149 (MQ=255)
                                                                                                                  tAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGc                               >  3:94781/1‑149 (MQ=255)
                                                                                                                        tCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGcc                         >  2:440886/1‑149 (MQ=255)
                                                                                                                         cAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCca                        <  1:285475/149‑1 (MQ=255)
                                                                                                                          aaCAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCcac                       >  2:596469/1‑149 (MQ=255)
                                                                                                                          aaCAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCcac                       >  2:484628/1‑149 (MQ=255)
                                                                                                                           aCAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCcaca                      >  1:176858/1‑149 (MQ=255)
                                                                                                                             aGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACg                    >  1:102703/1‑149 (MQ=255)
                                                                                                                                  ttcCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGaa               >  2:323920/1‑149 (MQ=255)
                                                                                                                                  ttcCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGaa               <  3:485905/149‑1 (MQ=255)
                                                                                                                                  ttcCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGaa               >  1:377362/1‑149 (MQ=255)
                                                                                                                                   tcCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGAAc              >  2:591599/1‑149 (MQ=255)
                                                                                                                                     ccAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGAACGa            >  1:219219/1‑149 (MQ=255)
                                                                                                                                     ccAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGAACGa            >  3:240040/1‑149 (MQ=255)
                                                                                                                                      cAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGAACGAc           >  1:360245/1‑149 (MQ=255)
                                                                                                                                       aGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGAACGAcc          >  2:485963/1‑149 (MQ=255)
                                                                                                                                       aGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGAACGAcc          >  2:77601/1‑149 (MQ=255)
                                                                                                                                       aGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATAATTACGCAGCTTCGCCACACGCAGAACGAcc          >  1:287697/1‑149 (MQ=255)
                                                                                                                                               aaGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGAACGACCAGCGGGaa  <  2:294731/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CTCATGCTTTGCTCGCCTCATGACGATTATGATCGTTTTGTTTTAAAGGAGGAAGTATCGGCAGTAGACGAATGAGTACGATGTAAGCGCAAATCTCAATAGCCACAAAACCAATAGAAATCAACAGTTCTTCCCAGGTCGGGAAGTAGGCGTAACCGCCGCCCGGGTTGAATGCCACCAGCGAATAGGTCAGACGCCAGGTTGCACAACCTAACAGTGCGCTCAGTGCTGACAGGAACAGCATGCGGGAATCATTACGCAGCTTCGCCACACGCAGAACGACCAGCGGGAA  >  NC_000913/3142985‑3143276

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