Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 793,331 T→G 6.0% E153A (GAG→GCG)  modF ← molybdate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913793,3310TG 6.0% 296.9 / 7.7 101E153A (GAG→GCG) modFmolybdate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (52/43);  new base G (4/2);  total (56/45)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 101 total aligned reads displayed.

CAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTC  >  NC_000913/793185‑793479
                                                                                                                                                  |                                                                                                                                                    
  aaCTCCGGGATCGCATCGAAGCGAGTGAGCCCCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCGGCTGACGGGAGGCAACATCCAGGGCCGCGACCGGCTCATCAAGAATCAACACGTCAGGCTCCgcc                                                                                                                                                 >  2:405935/1‑148 (MQ=255)
    atCCGGTATCTCTTCGAATTGATTTAGCACCAGTAGCAGATTAATACAGGAATGATGTATCGAGGCGTGTCGCTCAGCCAGATGCTGACGTGTGGGAACATCCAGTCCATCGAACGGCTCATCAAGAATAATCTAGTCAGGATCCGACAt                                                                                                                                               <  2:1544061/149‑1 (MQ=255)
          gATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAgcgc                                                                                                                                         >  1:706585/1‑150 (MQ=255)
             ctcATCGAAGCGATTGAGCACCAGTCCCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTGAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTg                                                                                                                                      <  2:78982/150‑1 (MQ=255)
             ctcATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTg                                                                                                                                      <  1:334741/150‑1 (MQ=255)
               cATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCGCTGATGTCGGAGCCTGac                                                                                                                                    <  2:1620659/150‑1 (MQ=255)
               cATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCGCTGATGTCGGAGCCTGac                                                                                                                                    <  2:1712765/150‑1 (MQ=255)
                 tCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGacac                                                                                                                                  <  1:1466984/150‑1 (MQ=255)
                 acGAAGCGATTCAGCACCAGTACCAGAGTAATACCGGACTGATGTACCGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGCCTCCGACATCAGCGCCTGacac                                                                                                                                  <  1:1066985/149‑1 (MQ=255)
                  cGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCATCAAGTCCGGCTCCGACATCAGCGCCTGacacc                                                                                                                                 >  1:23533/1‑149 (MQ=255)
                  cGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGAcaca                                                                                                                                 >  2:1213858/1‑150 (MQ=255)
                  cGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGAGCCTGAcaca                                                                                                                                 >  2:751425/1‑150 (MQ=255)
                       cGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAgg                                                                                                                            >  2:540886/1‑150 (MQ=255)
                       cGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAgg                                                                                                                            <  2:1278611/150‑1 (MQ=255)
                         tttGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGt                                                                                                                          <  2:1604763/149‑1 (MQ=255)
                         aTTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGt                                                                                                                          <  2:1189104/150‑1 (MQ=255)
                             tgCACCAGTACCAGATTAATACCGGACTGATGTAACGAGGCGAGGCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGttttt                                                                                                                      <  2:328031/149‑1 (MQ=255)
                             aGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACATAGCAGGGttttt                                                                                                                      >  2:917000/1‑150 (MQ=255)
                             aGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGttttt                                                                                                                      <  1:1548302/150‑1 (MQ=255)
                                 ccAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACCTCCAGGCCATCGAACGGCTCATCAAGGATCAACAAGTCAGGCTCCGCCATCAGCGCCTGACACAGCCGGGTTTTTcgcg                                                                                                                  >  2:314981/1‑150 (MQ=255)
                                   aGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGtc                                                                                                                >  1:1210168/1‑150 (MQ=255)
                                   aGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGGTTTTCGCGtc                                                                                                                >  1:183637/1‑150 (MQ=255)
                                   aGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCGCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGtc                                                                                                                <  1:1461243/150‑1 (MQ=255)
                                      aCCAGCGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCg                                                                                                             <  2:1531095/150‑1 (MQ=255)
                                       ccAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGc                                                                                                            >  1:513031/1‑150 (MQ=255)
                                       ccAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGc                                                                                                            >  2:735164/1‑150 (MQ=255)
                                            gTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTg                                                                                                       <  1:588042/150‑1 (MQ=255)
                                               aTACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGaa                                                                                                    >  1:427268/1‑150 (MQ=255)
                                               aTACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCATGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGaa                                                                                                    >  1:1124172/1‑150 (MQ=255)
                                                  ccGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAgg                                                                                                 >  2:52017/1‑150 (MQ=255)
                                                     gACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAt                                                                                              <  1:1666229/150‑1 (MQ=255)
                                                     gACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAt                                                                                              >  1:1794/1‑150 (MQ=255)
                                                     gACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAt                                                                                              <  1:1296874/150‑1 (MQ=255)
                                                       cTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAttt                                                                                            <  2:369136/150‑1 (MQ=255)
                                                          aTGTAACGAGTCGAGTAGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCACCAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTaaa                                                                                         <  2:1259827/150‑1 (MQ=255)
                                                            gTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGc                                                                                       <  1:494073/150‑1 (MQ=255)
                                                            gTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGc                                                                                       <  1:657765/150‑1 (MQ=255)
                                                               aCGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCgtc                                                                                    <  2:632163/150‑1 (MQ=255)
                                                                 gAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCgtcgg                                                                                  >  1:1581523/1‑150 (MQ=255)
                                                                 gAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGGCTCGCCAGTGGAAAGGTATTTAAAGCgtcgg                                                                                  >  1:1561200/1‑150 (MQ=255)
                                                                    gCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCgtcggtcg                                                                               <  1:540886/150‑1 (MQ=255)
                                                                     cGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGGATTTAAAGCgtcggtcgg                                                                              >  1:267771/1‑149 (MQ=255)
                                                                       aGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCTACAAGGCAGGCTCCGACATCAGCGCCTGACACAGCAGGGGTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGGCgagg                                                                            <  1:426813/150‑1 (MQ=255)
                                                                          cGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAATAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAATGGTATTTAAAGCGTCGGTCGaggagg                                                                         >  2:1448953/1‑150 (MQ=255)
                                                                          cGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGGAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGaggcgg                                                                         >  2:1079373/1‑150 (MQ=255)
                                                                          cGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAAAAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTGTTTCGCGTCTAGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGaggagg                                                                         >  2:1442155/1‑150 (MQ=255)
                                                                              cagccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCCGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAATGCGTCGGTCGAGGAGGGCgg                                                                     >  2:406716/1‑150 (MQ=255)
                                                                              cagccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCgg                                                                     <  2:529223/150‑1 (MQ=255)
                                                                                gccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGCCATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTAGTTAAAGCGTCGGTCGAGGAGGGCGGTa                                                                   >  1:667548/1‑150 (MQ=255)
                                                                                gccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTa                                                                   <  1:514742/150‑1 (MQ=255)
                                                                                 ccagcTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTaa                                                                  <  2:685436/150‑1 (MQ=255)
                                                                                     ctgctgACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATAAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATAcc                                                              <  1:419585/150‑1 (MQ=255)
                                                                                          gACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTTAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAAct                                                         >  2:213048/1‑150 (MQ=255)
                                                                                          gACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAAct                                                         >  2:1458932/1‑150 (MQ=255)
                                                                                            cGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgc                                                       >  1:1574081/1‑150 (MQ=255)
                                                                                            cGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgc                                                       >  1:822026/1‑150 (MQ=255)
                                                                                            cGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATAAGCGCCTGACACAGCAGGGGTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgc                                                       >  1:827595/1‑150 (MQ=255)
                                                                                             gTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgct                                                      <  2:1558530/150‑1 (MQ=255)
                                                                                             gTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAActgct                                                      <  2:1272865/150‑1 (MQ=255)
                                                                                                   cAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCa                                                >  1:200942/1‑150 (MQ=255)
                                                                                                     aCATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGt                                              >  1:1263348/1‑150 (MQ=255)
                                                                                                         ccAGGCCATCGAACGGCTCATCAAGAATCAACCAGTCAGGCTCCGCCATCAGCGCCTGACCCAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGCGGGCGGTAATACCGAACTGCTGCGCCCGTtgca                                          >  2:173395/1‑150 (MQ=255)
                                                                                                          cAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgcat                                         <  2:1124172/150‑1 (MQ=255)
                                                                                                          cAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgcat                                         <  1:651891/150‑1 (MQ=255)
                                                                                                           aGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTtgcatg                                        <  1:1552383/150‑1 (MQ=255)
                                                                                                              ccATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCaa                                     >  1:112386/1‑150 (MQ=255)
                                                                                                               cATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAAc                                    >  1:145981/1‑150 (MQ=255)
                                                                                                               cATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAAc                                    >  2:914785/1‑150 (MQ=255)
                                                                                                                    atcGGCTCATCACGACTCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTGTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGGCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGt                               <  1:202454/148‑1 (MQ=255)
                                                                                                                    aaCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGt                               <  2:548559/150‑1 (MQ=255)
                                                                                                                      cGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGGTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGc                             <  1:1213858/150‑1 (MQ=255)
                                                                                                                       ggCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCa                            >  1:1174894/1‑150 (MQ=255)
                                                                                                                          tcatcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATcc                         >  2:1541777/1‑150 (MQ=255)
                                                                                                                            atctAGAATCAACAAGTCAGGCTCCGATATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCtt                       <  1:961334/150‑1 (MQ=255)
                                                                                                                            atcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCtt                       >  2:1016236/1‑150 (MQ=255)
                                                                                                                             tcaagaatcaagaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCttt                      <  1:1669186/150‑1 (MQ=255)
                                                                                                                             tcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCACCGCGTTGCAtgcaac                      <  2:706585/150‑6 (MQ=255)
                                                                                                                             tcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCttt                      >  1:1449808/1‑150 (MQ=255)
                                                                                                                             tcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCttt                      >  1:741624/1‑150 (MQ=255)
                                                                                                                               aaGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTAc                    >  2:622300/1‑150 (MQ=255)
                                                                                                                               aaGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGGATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGAGGTGCATCCTTTag                    >  1:251739/1‑149 (MQ=255)
                                                                                                                                aGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACt                   >  1:1163415/1‑150 (MQ=255)
                                                                                                                                aGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACt                   <  2:1558416/150‑1 (MQ=255)
                                                                                                                                    ctgctgtgTCAGGCGCTGATGTCGGAGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCAt               >  2:597886/8‑150 (MQ=255)
                                                                                                                                     caacaaGTCCGGCTCCGACAGCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGGCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATc              <  2:183637/150‑1 (MQ=255)
                                                                                                                                     caacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTTAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATc              >  2:164613/1‑150 (MQ=255)
                                                                                                                                     caacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATc              >  1:414981/1‑150 (MQ=255)
                                                                                                                                     caacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATc              <  2:1210168/150‑1 (MQ=255)
                                                                                                                                       acaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCt            >  2:542976/1‑150 (MQ=255)
                                                                                                                                        tgtgTCAGGCGCTGATGTCGGAGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTg           >  1:514519/4‑150 (MQ=255)
                                                                                                                                         aaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGa          >  2:955916/1‑150 (MQ=255)
                                                                                                                                          aGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGaa         >  1:74895/1‑150 (MQ=255)
                                                                                                                                           gTCAGGCTCCGACATCAGCGCCTGACATAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAAt        <  1:917000/150‑1 (MQ=255)
                                                                                                                                           gTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAAt        <  1:1458932/150‑1 (MQ=255)
                                                                                                                                           gTCAGGCGCTGATGTCGGAGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAAt        >  1:1215600/1‑150 (MQ=255)
                                                                                                                                           gTCAGGCGCTGATGTCGGAGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAAt        >  1:624839/1‑150 (MQ=255)
                                                                                                                                             cAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCCTCCGGAATGa      >  2:450666/1‑150 (MQ=255)
                                                                                                                                             cAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGa      >  1:1343677/1‑150 (MQ=255)
                                                                                                                                                gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGAtct   >  1:470830/1‑150 (MQ=255)
                                                                                                                                                 cTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTAAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGAtctc  <  1:673023/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAAACTCCGGGATCTCATCGAAGCGATTGAGCACCAGTACCAGAGTAATACCGGACTGATGTAACGAGGCGAGTCGCTCAGCCAGCTGCTGACGTGAGGCAACATCCAGGCCATCGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTC  >  NC_000913/793185‑793479

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: