Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,432,812 A→T 5.8% L250* (TTG→TAG)  smf ← DNA recombination‑mediator A family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,432,8120AT 5.8% 223.8 / 3.7 69L250* (TTG→TAG) smfDNA recombination‑mediator A family protein
Reads supporting (aligned to +/- strand):  ref base A (29/36);  new base T (2/2);  total (31/38)
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.40e-01

AAATGAATTTTCAGGGGCGTCTGGCAACCAGTGCAATCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAA  >  NC_000913/3432667‑3432959
                                                                                                                                                 |                                                                                                                                                   
aaaTGAATTTTCAGGGGCGTCTGGCAACCAGTGCAATCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTggcaagg                                                                                                                                                  >  3:224610/1‑149 (MQ=255)
aaaTGAATTTTCAGGGGCGTCTGGCAACCAGTGCAATCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTggcaagg                                                                                                                                                  <  1:56452/149‑1 (MQ=255)
aaaTGAATTTTCAGGGGCGTCTGGCAACCAGTGCAATCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTggcaag                                                                                                                                                   <  1:307542/148‑1 (MQ=255)
  aTGAATTTTCAGGGGCGTCTGGCAACCAGTGCAATCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTggcaaggc                                                                                                                                                 >  3:220699/1‑148 (MQ=255)
            aGGGGCGTCTGGCAACCAGTGCAATCCAAATTGCATGTTTTCGAGAATTTCTTCCGGTTCCTTCACAATAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtctc                                                                                                                                      <  2:92089/149‑1 (MQ=255)
            aGGGGCGTCTGGCAACCAGTGCAATCCAAATTGCAAGTTTTCCAGAATTTCTTCGGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGTCAAAAACTtct                                                                                                                                       <  7:184480/148‑1 (MQ=255)
            aGGGGCGTCTGGCAACCAGTGCAATCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtctc                                                                                                                                      <  2:309706/149‑1 (MQ=255)
                gCGTCTGGCAACCAGTGCAATCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGccc                                                                                                                                  <  2:188724/149‑1 (MQ=255)
                      ggCAACCAGTGCAATCCAAATTGCGAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCt                                                                                                                             >  8:276390/1‑148 (MQ=255)
                       gCAACCAGTGCAATCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTc                                                                                                                            >  6:189578/1‑148 (MQ=255)
                       gCAACCAGGGCAATCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCt                                                                                                                             <  6:173258/147‑1 (MQ=255)
                                 cAATCCATTTTGCAAGTTTTCCAGTTTTTCTTCGGGTTGCGTCACTAGTATCGCCCCTTGTTTTATTAACCAGTGAGGGCCTTCGCTTCCCGTATTCCCTATTGGTCCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAAc                                                                                                                 <  4:220699/149‑1 (MQ=255)
                                 cAATCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAAc                                                                                                                 >  5:254570/1‑149 (MQ=255)
                                  aaTCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACg                                                                                                                <  8:277845/149‑1 (MQ=255)
                                   aTCCAAATTGCAAGTTTTCCAGATTTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGATGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGt                                                                                                               <  8:43810/149‑1 (MQ=255)
                                   aTCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtctc                                                                                                                                      >  1:214186/1‑126 (MQ=255)
                                   aTCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTtctc                                                                                                                                      <  2:214186/126‑1 (MQ=255)
                                    tCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTg                                                                                                              >  2:348078/1‑149 (MQ=255)
                                    tCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTg                                                                                                              >  3:130658/1‑149 (MQ=255)
                                      cAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGc                                                                                                             <  4:191948/148‑1 (MQ=255)
                                      cAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCt                                                                                                            <  8:319640/149‑1 (MQ=255)
                                                    tCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAAc                                                                                               <  6:254570/148‑1 (MQ=255)
                                                     ccAGAATTTCTTCAGATTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCa                                                                                             >  2:53646/1‑149 (MQ=255)
                                                       aGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCAc                                                                                            >  3:219849/1‑148 (MQ=255)
                                                       aGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGAAAAAACTTCTCGCCACTGCTCAATCGCACAACGTGCTGTCACCAGCGAccccc                                                                                            >  8:76293/1‑146 (MQ=255)
                                                       aGAATTTCTTCCGGATCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCt                                                                                                            <  6:153951/132‑1 (MQ=255)
                                                       aGAATTTCTTCCGGATCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCt                                                                                                            >  5:153951/1‑132 (MQ=255)
                                                             tcttcCGGTTCCGTCACAAGAATGGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCa                                                                                      <  1:301376/148‑1 (MQ=255)
                                                               ttcCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCACTATTGGACCTGGCAAGGCAAAAAATTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAg                                                                                   >  8:223218/1‑149 (MQ=255)
                                                                  cGGTGCCGTCACAAGAAGCGCACCTTGTTTTATTAACCAGTGAGGCCCTGCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGcc                                                                                 <  6:19678/148‑1 (MQ=255)
                                                                   ggTTCCGTCAGAAGAATCGCACGTTGTTTTATTCAGCAGTGAGGCCCTTCGCTTCCCGGATTCGCTATTGGACCTGGCAAGTCAAATTCTTCTGGCCCCTGCTCAAGCGCACAACTTGCTGTCACCAGCGAACCACTACGCAAAGCCGc                                                                               <  3:296436/149‑1 (MQ=255)
                                                                    gTTCCGTCAAAAGAATCGCACCTTGTTTTATTACCCAGTGAGGCCCTTCGCTTCCCGGCTTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCCCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCt                                                                              <  5:80721/149‑1 (MQ=255)
                                                                     ttCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGa                                                                                                 >  2:313504/1‑129 (MQ=255)
                                                                     ttCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGa                                                                                                 <  1:313504/129‑1 (MQ=255)
                                                                        cGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTcca                                                                          <  4:224610/149‑1 (MQ=255)
                                                                        cGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTcca                                                                          <  5:189578/149‑1 (MQ=255)
                                                                         gTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTccac                                                                         <  7:276390/149‑1 (MQ=255)
                                                                               aGAATCGCACCTTGTTTTATTAACCACTGAGGCCCTTCGCTTCCCGGAGTCCCTCCTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCACGCGCACAACGTACTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAg                                                                   <  4:28344/149‑1 (MQ=255)
                                                                                   tCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAAc                                                                                               <  6:234869/117‑1 (MQ=255)
                                                                                   tCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAAc                                                                                               >  5:234869/1‑117 (MQ=255)
                                                                                         cTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCa                                                                                      <  4:250054/120‑1 (MQ=255)
                                                                                         cTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCa                                                                                      >  3:250054/1‑120 (MQ=255)
                                                                                            gTTTTATTAGCCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACt                                                       >  3:229092/1‑148 (MQ=255)
                                                                                                tattaACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGa                                                   <  7:15557/148‑1 (MQ=255)
                                                                                                      ccAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGa                                            <  4:229092/149‑1 (MQ=255)
                                                                                                         gTGAGGCCCTTCGCTTCCCGGATTCCCGATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATaa                                         >  3:20658/1‑149 (MQ=255)
                                                                                                          tGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACg                                                                                        >  3:101650/1‑101 (MQ=255)
                                                                                                          tGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACg                                                                                        <  4:101650/101‑1 (MQ=255)
                                                                                                          tGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAAt                                        >  5:322855/1‑149 (MQ=255)
                                                                                                            aGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGc                                      <  3:154391/149‑1 (MQ=255)
                                                                                                                ccTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGAtt                                  >  3:253609/1‑149 (MQ=255)
                                                                                                                ccTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGAtt                                  >  4:85140/1‑149 (MQ=255)
                                                                                                                 cTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGAttt                                 >  2:35884/1‑149 (MQ=255)
                                                                                                                    cGCTTCCCGGCTTCCCTAGTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATttc                                <  6:322855/147‑1 (MQ=255)
                                                                                                                         cccGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTGCGTgg                         >  7:274690/1‑149 (MQ=255)
                                                                                                                            ggATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCATTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGa                       >  8:106378/1‑148 (MQ=255)
                                                                                                                             gATTCCCTATTGGACCTGGCTAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGaaa                     >  5:342531/1‑149 (MQ=255)
                                                                                                                             gATTCCCTATTGGACCTGGCTAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGaaa                     >  5:337158/1‑149 (MQ=255)
                                                                                                                                  ccTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGt                 <  1:115570/148‑1 (MQ=255)
                                                                                                                                  ccTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTATACATAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGt                 >  3:50213/1‑148 (MQ=255)
                                                                                                                                  ccTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCGCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTCCACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGt                 <  2:27855/148‑1 (MQ=255)
                                                                                                                                   cTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCAc                                               >  6:47700/1‑117 (MQ=255)
                                                                                                                                   cTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCAc                                               <  5:47700/117‑1 (MQ=255)
                                                                                                                                   cTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTaa               <  7:76293/149‑1 (MQ=255)
                                                                                                                                        ggACCTGGCTAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCaa           <  6:342531/148‑1 (MQ=255)
                                                                                                                                        ggACCTGGCTAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCaa           <  6:337158/148‑1 (MQ=255)
                                                                                                                                          aCCTGGCAAGGCAAAACCTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGCACCACTACGCAAAGCCGCTTCCACCACCAGTCCACCTTTACTTAGCCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAggg        >  6:98064/1‑149 (MQ=255)
                                                                                                                                          aCCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAggg        <  4:130658/149‑1 (MQ=255)
                                                                                                                                            cTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATaa                                         <  4:15027/114‑1 (MQ=255)
                                                                                                                                            cTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGt      <  2:366239/149‑1 (MQ=255)
                                                                                                                                            cTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGt      >  1:54244/1‑149 (MQ=255)
                                                                                                                                            cTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGt      >  1:366239/1‑149 (MQ=255)
                                                                                                                                             tGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTg     >  7:129257/1‑149 (MQ=255)
                                                                                                                                                caaggcaaAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGaa  >  1:341778/1‑149 (MQ=255)
                                                                                                                                                caaggcaaAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGaa  <  2:341778/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AAATGAATTTTCAGGGGCGTCTGGCAACCAGTGCAATCCAAATTGCAAGTTTTCCAGAATTTCTTCCGGTTCCGTCACAAGAATCGCACCTTGTTTTATTAACCAGTGAGGCCCTTCGCTTCCCGGATTCCCTATTGGACCTGGCAAGGCAAAAACTTCTCGCCCCTGCTCAAGCGCACAACGTGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAA  >  NC_000913/3432667‑3432959

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: