Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 792,560 G→A 5.2% A410V (GCG→GTG)  modF ← molybdate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913792,5600GA 5.2% 255.6 / 2.6 77A410V (GCG→GTG) modFmolybdate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (40/33);  new base A (2/2);  total (42/35)
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.57e-01

CAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAATGCCAATCGAATCAAAATAGCC  >  NC_000913/792417‑792705
                                                                                                                                               |                                                                                                                                                 
caaacaaCAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGa                                                                                                                                              <  2:402268/149‑1 (MQ=255)
     aacaaTTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACgctg                                                                                                                                         >  3:151540/1‑149 (MQ=255)
     aacaaTTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACgctg                                                                                                                                         >  4:164542/1‑149 (MQ=255)
     aacaaTTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACgctg                                                                                                                                         >  1:501023/1‑149 (MQ=255)
           tGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTcc                                                                                                                                   <  4:591798/149‑1 (MQ=255)
              gTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCa                                                                                                                                >  3:289093/1‑149 (MQ=255)
              gTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCa                                                                                                                                <  1:100132/149‑1 (MQ=255)
               tttCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAg                                                                                                                               <  1:58765/149‑1 (MQ=255)
                   aCCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGaaa                                                                                                                           >  4:6002/1‑149 (MQ=255)
                   aCCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGaaa                                                                                                                           <  3:80236/149‑1 (MQ=255)
                   aCCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGaaa                                                                                                                           >  2:337695/1‑149 (MQ=255)
                   aCCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGATTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGaaa                                                                                                                           >  2:701832/1‑149 (MQ=255)
                    ccTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAg                                                                                                                          >  1:32027/1‑149 (MQ=255)
                    ccTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAg                                                                                                                          >  4:80745/1‑149 (MQ=255)
                     cTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGa                                                                                                                         >  3:61579/1‑149 (MQ=255)
                      ttCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGAc                                                                                                                        >  3:303771/1‑149 (MQ=255)
                      ttCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGAc                                                                                                                        >  4:39495/1‑149 (MQ=255)
                             aTCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGaa                                                                                                                 <  4:532113/149‑1 (MQ=255)
                                aGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCACCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACgg                                                                                                              >  3:359911/1‑149 (MQ=255)
                                 gCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGa                                                                                                             <  2:582132/149‑1 (MQ=255)
                                 gCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGa                                                                                                             <  1:337695/149‑1 (MQ=255)
                                 gCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGa                                                                                                             <  3:39495/149‑1 (MQ=255)
                                     aTCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGt                                                                                                         >  1:662568/1‑149 (MQ=255)
                                      tCAACTAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGATCATAAGCAACGTCGGATGTCTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTc                                                                                                        >  4:582156/1‑149 (MQ=255)
                                      tCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGTCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTc                                                                                                        >  1:376943/1‑149 (MQ=255)
                                        aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGc                                                                                                           >  1:325745/1‑144 (MQ=255)
                                        aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGc                                                                                                           <  2:325745/144‑1 (MQ=255)
                                        aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAg                                                                                                      >  1:140366/1‑149 (MQ=255)
                                        aacaaAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAg                                                                                                      >  2:57022/1‑149 (MQ=255)
                                              aCGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTAt                                                                                                                     <  4:61100/128‑1 (MQ=255)
                                              aCGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTAt                                                                                                                     >  3:61100/1‑128 (MQ=255)
                                                   ggATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCACCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTa                                                                                           <  4:448742/149‑1 (MQ=255)
                                                     aTAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATc                                                                                         >  4:26669/1‑149 (MQ=255)
                                                       aaGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCaa                                                                                       <  1:38372/149‑1 (MQ=255)
                                                        aGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAAt                                                                                      <  2:501023/149‑1 (MQ=255)
                                                        aGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAAt                                                                                      <  4:151540/149‑1 (MQ=255)
                                                        aGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAAt                                                                                      <  4:476045/149‑1 (MQ=255)
                                                           tGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGcc                                                                                   >  3:625080/1‑149 (MQ=255)
                                                                     gCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTc                                                                                                        >  3:500419/1‑118 (MQ=255)
                                                                     gCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTc                                                                                                        <  4:500419/118‑1 (MQ=255)
                                                                        gATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGc                                                                      >  3:66259/1‑149 (MQ=255)
                                                                         aTCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTAt                                                                                          >  3:143391/1‑128 (MQ=255)
                                                                         aTCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTAt                                                                                          <  4:143391/128‑1 (MQ=255)
                                                                          tCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCa                                                                    >  2:1379/1‑149 (MQ=255)
                                                                          tCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCa                                                                    >  4:666374/1‑149 (MQ=255)
                                                                           cAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCAc                                                                   >  1:462953/1‑149 (MQ=255)
                                                                             aGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActg                                                                 >  1:347819/1‑149 (MQ=255)
                                                                              gCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActgc                                                                >  4:448792/1‑149 (MQ=255)
                                                                               ccccTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActgct                                                               <  3:26669/149‑1 (MQ=255)
                                                                                 ccTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACtgctgc                                                             <  3:103475/149‑1 (MQ=255)
                                                                                 ccTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACtgctgc                                                             <  3:1796/149‑1 (MQ=255)
                                                                                  cTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCa                                                            >  2:315005/1‑149 (MQ=255)
                                                                                       gTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGt                                                       >  3:230375/1‑149 (MQ=255)
                                                                                           ttCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGttttt                                                   <  2:451675/149‑1 (MQ=255)
                                                                                             cATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGc                                                 <  4:513891/149‑1 (MQ=255)
                                                                                             cATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGc                                                 <  3:7658/149‑1 (MQ=255)
                                                                                               tCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTg                                               >  2:659814/1‑149 (MQ=255)
                                                                                                cGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTgg                                              <  1:649836/149‑1 (MQ=255)
                                                                                                 gagaATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGc                                             >  3:174292/1‑149 (MQ=255)
                                                                                                   gaATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGa                                           >  3:166099/1‑149 (MQ=255)
                                                                                                    aaTAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGAt                                          <  3:656133/149‑1 (MQ=255)
                                                                                                       aaGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCg                                       >  2:485225/1‑149 (MQ=255)
                                                                                                           aaCGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAAc                                   <  2:347819/149‑1 (MQ=255)
                                                                                                           aaCGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAAc                                   >  2:565934/1‑149 (MQ=255)
                                                                                                               tCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGcc                               <  4:66259/149‑1 (MQ=255)
                                                                                                                      tttCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATaaa                        <  2:258041/149‑1 (MQ=255)
                                                                                                                       ttCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAAt                       <  4:289093/149‑1 (MQ=255)
                                                                                                                         cACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAATGc                     <  1:485225/149‑1 (MQ=255)
                                                                                                                              gTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAATGCCAATc                <  2:626311/149‑1 (MQ=255)
                                                                                                                               tGCGCGGACAATCAGCACCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAATGCCAATCg               <  4:359911/149‑1 (MQ=255)
                                                                                                                                        aaTCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAATGCCAATCGAATCAAAAt      >  1:87362/1‑149 (MQ=255)
                                                                                                                                        aaTCAGCACCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAATGCCAATCGAATCAAAAt      >  1:98041/1‑149 (MQ=255)
                                                                                                                                         aTCAGCGTCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAATGCCAATCGAATCAAAATa     <  2:376943/149‑1 (MQ=255)
                                                                                                                                          tCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAATGCCAATCGAATCAAAATAg    >  2:223741/1‑149 (MQ=255)
                                                                                                                                           cAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTGAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGGTTTTGCTGGCGATCCGAAACCGCCTGATAAATGCCAATCGAATCAAAATAGc   >  2:591331/1‑149 (MQ=255)
                                                                                                                                            aGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAATGCCAATCGAATCAAAATAGcc  >  1:131316/1‑149 (MQ=255)
                                                                                                                                            aGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTACTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAATGCCAATCGAATCAAAATAGcc  <  1:701832/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CAAACAACAATTGCGTTTCACCTTCGCTAATCAGCACATCAACAAAACGGCGGATAAGCTGGCGATTCAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCGATCCGAAACGGCCTGATAAATGCCAATCGAATCAAAATAGCC  >  NC_000913/792417‑792705

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: ‑